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This article describes the various experimental bounds on the 
variation of the fundamental constants of nature. After a discussion 
on the role of fundamental constants, of their definition and link 
with metrology, the various constraints on the variation of the fine 
structure constant, the gravitational, weak and strong interactions 
couplings and the electron to proton mass ratio are reviewed. This 
review aims (1) to provide the basics of each measurement, (2) to 
show as clearly as possible why it constrains a given constant and 
(3) to point out the underlying hypotheses. Such an investigation is 
of importance to compare the different results, particularly in view 
of understanding the recent claims of the detections of a variation of 
the fine structure constant and of the electron to proton mass ratio 
in quasar absorption spectra. The theoretical models leading to the 
prediction of such variation are also reviewed, including Kaluza- 
Klcin theories, string theories and other alternative theories and 
cosmological implications of these results are discussed. The links 
with the tests of general relativity arc emphasized. 
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I. INTRODUCTION 

The development of physics rehed considerably on the 
Copernician principle, which states that we are not living 
in a particular place in the universe and stating that the 
laws of physics do not differ from one point in spacetime 
to another. This contrasts with the Aristotelian point of 
view in which the laws on Earth and in Heavens differ. It 
is however natural to question this assumption. Indeed, 
it is difficult to imagine a change of the form of physical 
laws (e.g. a Newtonian gravitation force behaving as the 
inverse of the square of the distance on Earth and as 
another power somewhere else) but a smooth change in 
the physical constants is much easier to conceive. 

Comparing and reproducing experiments is also a root 
of the scientific approach which makes sense only if the 
laws of nature does not depend on time and space. This 
hypothesis of constancy of the constants plays an impor- 
tant role in particular in astronomy and cosmology where 
the redshift measures the look-back time. Ignoring the 
possibility of varying constants could lead to a distorted 
view of our universe and if such a variation is estab- 
lished corrections would have to be applied. It is thus 
of importance to investigate this possibility especially as 
the measurements become more and more precise. Obvi- 
ously, the constants have not undergone huge variations 
on Solar system scales and geological time scales and one 
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is looking for tiny effects. Besides, the question of the 
values of the constants is central to physics and one can 
hope to explain them dynamically as predicted by some 
high-energy theories. Testing for the constancy of the 
constants is thus part of the tests of general relativity. 
Let us emphasize that this latter step is analogous to the 
transition from the Newtonian description of mechanics 
in which space and time were just a static background in 
which matter was evolving to the relativistic description 
where spacctime becomes a dynamical quantity deter- 
mined by the Einstein equations (Damour, 2001). 

Before discussing the properties of the constants of na- 
ture, we must have an idea of which constants to consider. 
First, all constants of physics do not play the same role, 
and some have a much deeper one than others. Follow- 
ing Levy-Lcblond (1979), we can define three classes of 
fundamental constants, class A being the class of the con- 
stants characteristic of particular objects, class B being 
the class of constants characteristic of a class of physi- 
cal phenomena, and class C being the class of universal 
constants. Indeed, the status of a constant can change 
with time. For instance, the velocity of light was a ini- 
tially a type A constant (describing a property of light) 
which then became a type B constant when it was real- 
ized that it was related to electro-magnetic phenomena 
and, to finish, it ended as a type C constant (it enters 
many laws of physics from electromagnetism to relativity 
including the notion of causality...). It has even become 
a much more fundamental constant since it has been cho- 
sen as the definition of the meter (Petley, 1983). A more 
conservative definition of a fundamental constant would 
thus be to state that it is any parameter that can not be 
calculated with our present knowledge of physics, i.e. a 
free parameter of our theory at hand. Each free param- 
eter of a theory is in fact a challenge for future theories 
to explain its value. 

How many fundamental constants should wc consider? 
The set of constants which are conventionally consid- 
ered as fundamental (Flowers and Petley, 2001) con- 
sists of the electron charge e, the electron mass mc, the 
proton mass mp, the reduced Planck constant ?i, the 
velocity of light in vacuum c, the Avogadro constant 
A^^, the Boltzmann constant fcg, the Newton constant 
G, the permeability and permittivity of space, eo and 
liQ. The latter has a fixed value in the SI system of 
unit (/xo = 47r X 10~^Hm^^) which is implicit in the 
definition of the Ampere; eo is then fixed by the rela- 
tion ffoMo = . The inclusion of Nj^ in the former 
list has been debated a lot (see e.g. Birge, 1929). To 
compare with, the minimal standard model of particle 
physics plus gravitation that describes the four known 
interactions depends on 20 free parameters (Cahn, 1996; 
Hogan, 2000): the Yukawa coefficients determining the 
masses of the six quark (u, d, c, s, t, b) and three lepton 
(e,/U, r) flavors, the Higgs mass and vacuum expecta- 



tion value, three angles and a phase of the Cabibbo- 
Kobayashi-Maskawa matrix, a phase for the QCD vac- 
uum and three coupling constants gs,g^, gi for the gauge 
group SU{3) X SU{2) x f/(l) respectively Below the Z 
mass, gi and combine to form the electro-magnetic 
coupling constant 

_o 5 _2 —2 /-I \ 

^BM = 3-91 + .9w • (1) 

The number of free parameters indeed depends on the 
physical model at hand (see Weinberg, 1983). This is- 
sue has to be disconnected from the number of required 
fundamental dimensionful constants. Duff, Okun and 
Veneziano (2002) recently debated this question, respec- 
tively arguing for none, three and two (see also Wignall, 
2000). Arguing for no fundamental constant leads to 
consider them simply as conversion parameters. Some 
of them are, like the Boltzmann constant, but some oth- 
ers play a deeper role in the sense that when a physi- 
cal quantity becomes of the same order of this constant 
new phenomena appear, this is the case e.g. of h and c 
which are associated respectively to quantum and rela- 
tivistic effects. Okun (1991) considered that only three 
fundamental constants are necessary, the underlying rea- 
son being that in the international system of units which 
has 7 base units an 17 derived units, four of the seven 
base units are in fact derived (Ampere, Kelvin, mole and 
candela). The three remaining base units (meter, second 
and kilogram) are then associated to three fundamental 
constants (c, h and G). They can be seen as limiting 
quantities: c is associated to the maximum velocity and 
U to the unit quantum of angular momentum and sets 
a minimum of uncertainty whereas G is not directly as- 
sociated to any physical quantity (see Martins 2002 who 
argues that G is the limiting potential for a mass that 
does not form a black hole). In the framework of quan- 
tum field theory -|- general relativity, it seems that this 
set of three constants has to be considered and it al- 
lows to classify the physical theories (see figure 1). How- 
ever, Veneziano (1986) argued that in the framework of 
string theory one requires only two dimensionful funda- 
mental constants, c and the string length Ag. The use 
of h seems unnecessary since it combines with the string 
tension to give Ag. In the case of the Goto-Nambu ac- 
tion S/h = (T/n) Jd{Area) = Xj'^ J d{Area) and the 
Planck constant is just given by A^"^. In this view, h has 
not disappeared but has been promoted to the role of a 
UV cut-off that removes both the infinities of quantum 
field theory and singularities of general relativity. This 
situation is analogous to pure quantum gravity (Novikov 
and Zel'dovich. 1982) where h and G never appear sepa- 
rately but only in the combination i^^ = Ghjc? so that 
only c and are needed. Volovik (2002) considered the 
analogy with quantum liquids. There, an observer knows 
both the effective and microscopic physics so that he can 
judge whether the fundamental constants of the effective 
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theory remain fundamental constants of the microscopic 
theory. The status of a constant depends on the consid- 
ered theory (effective or microscopic) and, more interest- 
ingly, on the observer measuring them, i.e. on whether 
this observer belongs to the world of low-energy quasi- 
particles or to the microscopic world. 

Resolving this issue is indeed far beyond the scope 
of this article and can probably be considered more as 
an epistemological question than a physical one. But, 
as the discussion above tends to show, the answer de- 
pends on the theoretical framework considered [see also 
Cohen- Tannoudji (1985) for arguments to consider the 
Boltzmann constant as a fundamental constant] . A more 
pragmatic approach is then to choose a theoretical frame- 
work, so that the set of undetermined fixed parameters 
is fully known and then to wonder why they have the 
values they have and if they are constant. 

We review in this article both the status of the ex- 
perimental constraints on the variation of fundamental 
constants and the theoretical motivations for considering 
such variations. In section II, we recall Dirac's argument 
that initiated the consideration of time varying constants 
and we briefly discuss how it is linked to anthropic ar- 
guments. Then, since the fundamental constants are en- 
tangled with the theory of measurement, we make some 
very general comments on the consequences of metrology. 
In Sections III and IV, we review the observational con- 
straints respectively on the variation of the fine structure 
and of gravitational constants. Indeed, we have to keep 
in mind that the obtained constraints depend on under- 
lying assumptions on a certain set of other constants. We 
summarize more briefly in Section V, the constraints on 
other constants and we give, in Section VI, some hints 
of the theoretical motivations arising mainly from grand 
unified theories, Kaluza-Klein and string theories. We 
also discuss a number of cosmological models taking these 
variations into account. For recent shorter reviews, see 
Varshalovich et al. (2000a), Chiba (2001), Uzan (2002) 
and Martins (2002). 

Notations: In this work, we use SI units and the follow- 
ing values of the fundamental constants today^ 

c = 299,792,458m-s-^ (2) 

h = 1.054571596(82) x 10"^^ J • s (3) 

G = 6.673(10) X IQ-^^ m^ • kg"^ • s"^ (4) 

me = 9.10938188(72) x lO'^i kg (5) 

mp = 1.67262158(13) x lO^^Tj^g (g) 

mn = 1.67492716(13) x IO^^t (y-, 

e = 1.602176462(63) x 10"^^ C (8) 



^soo http://physics.nist.gov/cuu/Constaiits/ for an up 
to date list of the recommended values of the constants of 
nature. 



for the velocity of light, the reduced Planck constant, the 
Newton constant, the masses of the electron, proton and 
neutron, and the charge of the electron. We also define 





2 ^ 


(9) 


and the following 


; dimensionless ratios 






— ~ 1/137.03599976(50) 


(10) 


"w = 


G„rnf.c ^ 
~ 1.03 X 10"^ 


(11) 


o/ (E\ = 
"sl-^J — 


he 




= 


Grrii ,„ 
P ~ 5 X 10-39 

nc 


(13) 


n = 


— ~ 5.44617 X 10"^ 

TOp 


(14) 


X = 


5P"L-" ~ 1-62 X 10-^ 


(15) 


y = 


g^al^ ~ 2.977 x 10"^ 


(16) 



which characterize respectively the; strength of the 
electro-magnetic, weak, strong and gravitational forces 
and the electron-proton mass ratio, gp ~ 5.585 is the 
proton gyro-magnetic factor. Note that the relation (12) 
between two quantities that depend strongly on energy; 
this will be discussed in more details in Section V. We 
introduce the notations 

ao = ^-^^ — = 0.5291771 A (17) 



-Ei = Iruec^al,, = 13.60580 eV 



(18) 



E 

R^ = --l = 1.0973731568549(83) x 10^ m'^ (19) 
he 

respectively for the Bohr radius, the hydrogen ionization 
energy and the Rydbcrg constant. 

While working in cosmology, we assume that the uni- 
verse is described by a Friedmann-Lemaitre spacetime 



(20) 



where t is the cosmic time, a the scale factor and 7^ the 
metric of the spatial sections. We define the redshift as 



l + z = 



a 



(21) 



where oq is the value of the scale factor today while v^, 
and uq are respectively the frequencies at emission and 
today. We decompose the Hubble constant today as 



if-i = 9.7776 X 10^ /i-i yr 



(22) 



where h = 0.68 ±0.15 is a dimensionless number, and the 
density of the universe today is given by 

po = 1-879 X IQ-^'^nh'^ kg • m-3. (23) 
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II. GENERALITIES 

A. From Dirac nutnerologlcal principle to anthropic 
arguments 

The question of the constancy of the constants of 
physics was probably first addressed by Dirac (1937, 
1938, 1979) who expressed, in his "Large Numbers hy- 
pothesis" , the opinion that very large (or small) dimen- 
sionless universal constants cannot be pure mathemat- 
ical numbers and must not occur in the basic laws of 
physics. He suggested, on the basis of this numerological 
principle, that these large numbers should rather be con- 
sidered as variable parameters characterizing the state 
of the universe Dirac formed the five dimensionless ra- 
tios a^„, a^, a^, 5 = Hoh/mpC^ ~ 2/i x 10"'*^ and 
e = Gpo/Hq ~ 5/z~^ X 10~* and asked the question of 
which of these ratio is constant as the universe evolves. 
Usually, only 6 and e vary as the inverse of the cosmic 
time (note that with the value of the density chosen by 
Dirac, the imiverse is not flat so that a (x t and p oc t^^). 
Dirac then noticed that a^^n/a^^j^, representing the rel- 
ative magnitude of electrostatic and gravitational forces 
between a proton and an electron, was of the same or- 
der as Hoe-^/meC^ = Sa^^/p, representing the age of the 
universe in atomic time so that the five previous num- 
bers can be "harmonized" if one assumes that and 
S vary with time and scale as the inverse of the cosmic 
time^ . This implies that the intensity of all gravitational 
effects decrease with a rate of about 10""'^'' yr"'^ and that 
p oc t~'^ (since e is constant) which corresponds to a 
flat universe. Kothari (1938) and Chandrasekhar (1939) 
were the first to point out that some astronomical con- 
sequences of this statement may be detectable. Similar 
ideas were expressed by Milne (1937). 

Dicke (1961) pointed out that in fact the density of the 
universe is determined by its age, this age being related 
to the time needed to form galaxies, stars, heavy nuclei... 
This led him to formulate that the presence of an ob- 
server in the universe places constraints on the physical 
laws that can be observed. In fact, what is meant by 
observer is the existence of (highly?) organized systems 
and the anthropic principle can be seen as a rephrasing of 
the question "why is the universe the way it is?" (Hogan, 
2000). Carter (1974, 1976, 1983), who actually coined 
the term "anthropic principle" , showed that the numero- 
logical coincidence found by Dirac can be derived from 
physical models of stars and the competition between the 
weakness of gravity with respect to nuclear fusion. Carr 



and Rees (1979) then showed how one can scale up from 
atomic to cosmological scales only by using combinations 
of "em> "g and me/mp. 

Dicke (1961, 1962b) brought Mach's principle into the 
discussion and proposed (Brans and Dicke, 1961) a the- 
ory of gravitation based on this principle. In this the- 
ory the gravitational constant is replaced by a scalar 
field which can vary both in space and time. It follows 
that, for cosmological solutions, G oc H oc and 
p oc where n is expressible in terms of an arbitrary 
parameter oj-g-^ as = 2 + 3oj-g-^/2. Einstein gravity is 
recovered when w^^ oc. This predicts that oc 
and S oc whereas a^Mi ^ kept constant. 

This kind of theory was further generalized to obtain var- 
ious functional dependences for G in the formalization of 
scalar-tensor theories of gravitation (see e.g. Damour and 
Esposito-Farese, 1992). 

The first extension of Dirac's idea to non-gravitational 
forces was proposed by Jordan (1937, 1939) who still 
considered that the weak interaction and the proton to 
electron mass ratio were constant. He realized that the 
constants has to become dynamical fields and used the 
action 



s = l^/^dV ^"Kt) 



(24) 



^The ratio Sa^y^/p, represents roughly the inverse of the 
number of times an electron orbits around a proton during 
the age of the universe. Already, this suggested a link be- 
tween micro-physics and cosmological scales. 



r] and ^ being two parameters. Fierz (1956) realized that 
with such a Lagrangian, atomic spectra will be space- 
dependent. But, Dirac's idea was revived after Teller 
(1948) argued that the decrease of G contradicts paleon- 
tological evidences [see also Pochoda and Schwarzschild 
(1964) and Gamow (1967c) for evidences based on the 
nuclear resources of the Sun]. Gamow (1967a, 1967b) 
proposed that a^^ might vary as t in order to save 
the, according to him, "elegant" idea of Dirac (see also 
Stanyukovich, 1962). In both Gamow (1967a, 1967b) 
and Dirac (1937) theories the ratio a^/a^^ decreases as 
t-^ Teller (1948) remarked that o'^^ -InFotp, so 
that would become the logarithm of a large num- 

EM ^ ^ 

ber. Landau (1955), de Witt (1964) and Isham et al. 

(1971) advocated that such a dependence may arise if 
the Planck length provides a cut-off to the logarithmic 
divergences of quantum electrodynamics. In this latter 
class of models a^-^^ oc 1/lnt, a^. oc t~^, S CK and 
and e remain constant. Dyson (1967), Peres (1967) and 
then Davies (1972) showed, using respectively geological 
data of the abundance of rhenium and osmium and the 
stability of heavy nuclei, that these two hypothesis were 
ruled out observationally (see Section III for details on 
the experimental results). Modern theories of high en- 
ergy physics offer new arguments to reconsider the varia- 
tion of the fundamental constants (see Section VI). The 
most important OTitcome of Dirac's proposal and of the 
following assimilated theories [among which a later ver- 
sion of Dirac (1974) theory in which there is matter cre- 
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ation either where old matter was present or uniformly 
throughout the universe] is that the hypothesis of the 
constancy of the fundamental constants can and must be 
checked experimentally. 

A way to reconcile some of the large numbers is to 
consider the energy dependence of the couplings as de- 
termined by the renormalization group (see e.g. Itzkyson 
and Zuber, 1980). For instance, concerning the fine struc- 
ture constant, the energy-dependence arises from vacuum 
polarization that tends to screen the charge. This screen- 
ing is less important at small distance and the charge ap- 
pears bigger so that the effective coupling constant grows 
with energy. It follows from this approach that the three 
gauge groups get unified into a larger grand unification 
group so that the three couplings ; ctw ^"^^ stem 
from the same dimensionless number a^.-^,^. This might 
explain some large numbers and answer some of Dirac 
concerns (Hogan, 2000) but indeed, it does not explain 
the weakness of gravity which has become known as the 
hierarchy problem. 

Let us come back briefly to the anthropic considera- 
tions and show that they allow to set an interval of ad- 
missible values for some constants. Indeed, the anthropic 
principle does not tell whether the constants are varying 
or not but it gives an insight on how special our uni- 
verse is. In such an approach, one studies the effect of 
small variations of a constant around its observed value 
and tries to find a phenomenon highly dependent on this 
constant. This does not ensure that there is no other set 
of constants (very different of the one observed today) 
for which an organized universe may exist. It just tells 
about the stability in a neighborhood of the location of 
our universe in the parameter space of physical constants. 
Rozental (1988) argued that requiring that the lifetime 
of the proton Tp ~ a~^{n / m^c^) c-x:p(l / a^„) ~ lO^^ yr 
is larger than the age of the universe ~ c/Hq ~ 10^^ s 
implies that a-^^^ < 1/80. On the other side, if we be- 
lieve in a grand unified theory, this unification has to take 
place below the Planck scale implying that a^,^ > 1/170, 
this bound depending on assumptions on the particle con- 
tent. Similarly requiring that the electromagnetic repul- 
sion is much smaller than the attraction by strong in- 
teraction in nuclei (which is necessary to have nuclei) 
leads to a-^^ < 1/20. The thermonuclear reactions in 
stars are efficient if fc^T ^ a^jj^mpC^ and the temper- 
ature of a star of radius and mass Mg can roughly 
be estimated as k^^T ~ GMsfUp/ R,^, which leads to the 
estimate a^,^ ~ 10~^. One can indeed think of many 
other examples to put such bounds. From the previous 
considerations, we retain that the most stringent is 

1/170 < a,^ < 1/80. (25) 

It is difficult to believe that these arguments can lead to 
much sharper constraints. They are illustrative and give 
a hint that the constants may not be "random" parame- 



ters without giving any explanation for their values. 

Rozental (1988) also argued that the existence of hy- 
drogen and the formation of complex elements in stars 
(mainly the possibility of the reaction 3a ^^C) set 
constraints on the values of the strong coupling constant. 
The production of ^^C in stars requires a triple tuning: 
(i) the decay lifetime of *Be, of order 10~^ s, is four orders 
of magnitude longer than the time for two a particles to 
scatter, (ii) an excited state of the carbon lies just above 
the energy of *Be + a and finally (iii) the energy level of 
^^O at 7.1197 MeV is non resonant and below the energy 
of ^^C-I-q:, of order 7.1616 MeV, which ensures that most 
of the carbon synthetized is not destroyed by the capture 
of an a-particle (see Livio et al, 2000). Ober hummer et 
al. (2000) showed that outside a window of respectively 
0.5% and 4% of the values of the strong and electromag- 
netic forces, the stellar production of carbon or oxygen 
will be reduced by a factor 30 to 1000 (see also Pochet et 
al, 1991). Concerning the gravitational constant, galaxy 
formation require < 10^. Other such constraints on 
the other parameters listed in the previous section can 
be obtained. 

B. Metrology 

The introduction of constants in physical law is closely 
related to the existence of systems of units. For instance, 
Newton's law states that the gravitational force between 
two masses is proportional to each mass and inversely 
proportional to their separation. To transform the pro- 
portionality to an equality one requires the use of a quan- 
tity with dimension of m^ • kg""'^ • independent of the 
separation between the two bodies, of their mass, of their 
composition (equivalence principle) and on the position 
(local position invariance) . With an other system of units 
this constant could have simply been anything. 

The determination of the laboratory value of constants 
relies mainly on the measurements of lengths, frequen- 
cies, times,... (see Petley, 1985 for a treatise on the mea- 
surement of constants and Flowers and Petley, 2001, for 
a recent review) . Hence, any question on the variation of 
constants is linked to the definition of the system of units 
and to the theory of measurement. The choice of a base 
units affects the possible time variation of constants. 

The behavior of atomic matter is mainly determined 
by the value of the electron mass and of the fine structure 
constant. The Rydberg energy sets the (non-relativistic) 
atomic levels, the hyperfine structure involves higher 
powers of the fine structure constant, and molecular 
modes (including vibrational, rotational... modes) depend 
on the ratio rriQ/mp. As a consequence, if the fine struc- 
ture constant is spacetime dependent, the comparison be- 
tween several devices such as clocks and rulers will also 
be spacetime dependent. This dependence will also differ 
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from one clock to another so that metrology becomes both 
device and spacetime dependent. 

Besides this first metrologic problem, the choice of 
units has implications on the permissible variations of 
certain dimensionful constant. As an illustration, we fol- 
low Petley (1983) who discusses the implication of the 
definition of the meter. The definition of the meter via 
a prototype platinum-iridium bar depends on the inter- 
atomic spacing in the material used in the construction 
of the bar. Atkinson (1968) argued that, at first order, it 
mainly depends on the Bohr radius of the atom so that 
this definition of the meter fixes the combination (17) 
as constant. Another definition was based on the wave- 
length of the orange radiation from krypton-86 atoms. 
It is likely that this wavelength depends on the Rydberg 
constant and on the reduced mass of the atom so that it 
ensures that rUec'^a-^ /2h is constant. The more recent 

^ EM / 

definition of the meter as the length of the path traveled 

by light in vacuum during a time of 1/299792458 of a 
second imposes the constancy of the speed of light^ c. 
Identically, the definitions of the second as the duration 
of 9,192,631,770 periods of the transition between two 
hyperfine levels of the ground state of cesium-133 or of 
the kilogram via an international prototype respectively 
impose that rri^c^ ot^^^ j fi and mp are fixed. 

Since the definition of a system of units and the value 
of the fundamental constants (and thus the status of their 
constancy) are entangled, and since the measurement of 
any dimensionful quantity is in fact the measurements of 
a ratio to standards chosen as units, it only makes sense 
to consider the variation of dimensionless ratios. 

In theoretical physics, we often use the fundamental 
constants as units (see McWeeny, 1973 for the relation 
between natural units and SI units). The international 
system of units (SI) is more appropriate to human size 
measurements whereas natural systems of units are more 
appropriate to the physical systems they refer to. For 
instance h, c and G allows to construct the Planck mass, 
time and length which are of great use as units while 
studying high-energy physics and the same can be done 
from h, e, and Eq to construct a unit mass (too)j 
length (47reo/i^/rr/ce^) and time (2eo/i^/7rmee^). A phys- 
ical quantity can always be decomposed as the product of 
a label representing a standard quantity of reference and 
a numerical value representing the number of times the 
standard has to be taken to build the required quantity. 
It follows that a given quantity X that can be expressed 
as X = fcii^i(m, kg, s, . . .) with fci a dimensionless quan- 
tity and Fi a function of the base units (here SI) to some 
power. Let us decompose X as X = k2F2{h,e,c, . . .) 



^Note that the velocity of light is not assigned a fixed value 
directly, but rather the value is fixed as a consequence of the 
definition of the meter. 



where k2 is another dimensionless constant and F2 a func- 
tion of a sufficient number of fundamental constants to be 
consistent with the initial base units. The time variation 
of X is given by 

dlnX _ dlnfci dlnFi _ dlnfca dlnFa 
dt ~ ~dt ^ dr~ ~ ~dt ^ dt~' 
Since only dk\/dt or dk2/dt can be measured, it is nec- 
essary to have chosen a system of units, the constancy 
of which is assumed (i.e. that either dFi/dt = or 
dF2/dt = 0) to draw any conclusion concerning the time 
variation oi X, in the same way as the description of a 
motion needs to specify a reference frame. 

To illustrate the importance of the choice of units and 
the entanglement between experiment and theory while 
measuring a fundamental constant, let us sketch how one 
determines rUe in the SI system (following Mohr and Tay- 
lor, 2001), i.e. in kilogram (see figure 2). The kilo- 
gram is defined from a platinum-iridium bar to which 
we have to compare the mass of the electron. The key 
to this measurement is to express the electron mass as 
me = 2hRoa/a'^^c. From the definition of the second, 
i?oo is determined by precision laser-spectroscopy on hy- 
drogen and deuterium and the theoretical expression for 
the ls-2s hydrogen transition as v = (3/4)i?ooc[l — H + 
lla|„/48+(56a3^)/(97r) Ina,^ + . . .] arising from QED. 
The fine structure constant is determined by compar- 
ing theory and experiment for the anomalous magnetic 
moment of the electron (involving again QED). Finally, 
the Planck constant is determined by a Watt balance 
comparing a Watt electrical power to a Watt mechanical 
power (involving classical mechanics and classical elec- 
tromagnetism only: h enters through the current and 
voltage calibration based on two condensed-matter phe- 
nomena: Josephson and quantum Hall effects so that it 
involves the theories of these two effects). 

As a conclusion, let us recall that (i) in general, the val- 
ues of the constants are not determined by a direct mea- 
surement but by a chain involving both theoretical and 
experimental steps, (ii) they depend on our theoretical 
understanding, (iii) the determination of a self-consistent 
set of values of the fundamental constants results from an 
adjustment to achieve the best match between theory and 
a defined set of experiments (see e.g., Birge, 1929) (iv) 
that the system of units plays a crucial role in the mea- 
surement chain, since for instance in atomic units, the 
mass of the electron could have been obtained directly 
from a mass ratio measurement (even more precise!) and 
(v) fortunately the test of the variability of the constants 
does not require a priori to have a high-precision value 
of the considered constant. 

In the following, we will thus focus on the variation of 
dimensionless ratios which, for instance, characterize the 
relative magnitude of two forces, and are independent of 
the choice of the system of units and of the choice of stan- 
dard rulers or clocks. Let us note that some (hopeless) 



6 



attempts to constraint the time variation of dimcnsionful 
constants have been tried and will be briefly discussed in 
Section V.F. This does not however mean that a physical 
theory cannot have dimcnsionful varying constants. For 
instance, a theory of varying fine structure constant can 
be implemented either as a theory with varying electric 
charge or varying speed of light. 

C. Overview of the methods 

Before going into the details of the constraints, it is 
worth taking some time to discuss the kind of experi- 
ments or observations that we need to consider and what 
we can hope to infer from them. 

As emphasized in the previous section, we can only 
measure the variation of dimensionless quantities (such 
as the ratio of two wavelengths, two decay rates, two 
cross sections ...) and the idea is to pick up a physical 
system which depends strongly on the value of a set of 
constants so that a small variation will have dramatic 
effects. The general strategy is thus to constrain the 
spacetime variation of an observable quantity as precisely 
as possible and then to relate it to a set of fundamental 
constants. 

Basically, we can split all the methods into three 
classes: (i) atomic methods including atomic clocks, 
quasar absorption spectra and observation of the cos- 
mic microwave background radiation (CMBR) where one 
compares ratios of atomic transition frequencies. The 
CMB observation depends on the dependence of the re- 
combination process on Q!j,„; (ii) nuclear methods in- 
cluding nucleosynthesis, a- and /3-decay, Oklo reactor for 
which the observables are respectively abundances, life- 
times and cross sections; and (iii) gravitational methods 
including the test of the violation of the universality of 
free fall where one constrains the relative acceleration of 
two bodies, stellar evolution. . . 

These methods are either experimental (e.g. atomic 
clocks) for which one can have a better control of the sys- 
tematics, observational (e.g. geochemical, astrophysical 
and cosmological observations) or mixed (a- and /3-decay, 
universality of free fall). This sets the time scales on 
which a possible variation can be measured. For instance, 
in the case of the fine structure constant (see Section III) , 
one expects to be able to constrain a relative variation 
of a^jM of order 10~^ [geochemical (Oklo)], 10~^ [astro- 
physical (quasars)], 10"'^ — 10~^ [cosmological methods] , 
10~^^ — 10~^^ [laboratory methods] respectively on time 
scales of order 10^ yr, 10^ - 10^° yr, 10^° yr and 1-12 
months. This brings up the question of the comparison 
and of the compatibility of the different measurements 
since one will have to take into account e.g. the rate of 
change of 0:^,^; which is often assumed to be constant. In 
general, this requires to specify a model both to deter- 



mine the law of evolution and the links between the con- 
stants. Long time scale experiments allow to test a slow 
drift evolution while short time scale experiments enable 
to test the possibility of a rapidly varying constant. 

The next step is to convert the bound on the varia- 
tion of some measured physical quantities (decay rate, 
cross section,...) into a bound on some constants. It is 
clear that in general (for atomic and nuclear methods 
at least) it is impossible to consider the electromagnetic, 
weak and strong effects independently so that this latter 
step involves some assumptions. 

Atomic methods are mainly based on the comparison 
of the wavelengths of different transitions. The non rel- 
ativistic spectrum depends mainly on Roo and jj, the 
fine structure on R^oa'^^ and the hyperfine structure on 
9pRooCf^^^- Extending to molecular spectra to include 
rotational and vibrational transitions allows to have ac- 
cess to /U. It follows that we can hope to disentangle the 
observations of the comparisons of different transitions 
to constrain on the variation of {a^^, jj,, gp). The excep- 
tion is CMB which involves a dependence on a^^ and 
mainly due to the Thomson scattering cross section and 
the ionization fraction. Unfortunately the effect of these 
parameters have to be disentangled from the dependence 
on the usual cosmological parameters which render the 
interpretation more difficult. 

The internal structure and mass of the proton and neu- 
tron are completely determined by strong gauge fields 
and quarks interacting together. Provided we can ignore 
the quark masses and electromagnetic effects, the whole 
structure is only dependent on an energy scale A^^j, . It 
follows that the stability of the proton greatly depends 
on the electromagnetic effects and the masses and 
of the up and down quarks. In nuclei, the interaction of 
hadrons can be thought to be mediated by pions of mass 
~ mp(mu -I- rud)- Since the stability of the nucleus 
mainly results from the balance between this attractive 
nuclear force, the nucleon degeneracy pressure and the 
Coulomb repulsion, it will mainly involve m^, mj, Q^em- 

Big bang nucleosynthesis depends on G (expansion 
rate), Gp (weak interaction rates), ctg (binding of light 
elements), a^^ (via the electromagnetic contribution to 
rrin — m-p but one will also have to take into account the 
contribution of a possible variation of the mass of the 
quarks, and Wd). Besides, if to„ — rrip falls below 
rric the /?-decay of the neutron is no longer energetically 
possible. The abundance of helium is mainly sensitive to 
the freeze-out temperature and the neutron lifetime and 
heavier element abundances to the nuclear rates. 

All nuclear methods involve a dependence on the mass 
of the nuclei of charge Z and atomic number A 

m{A, Z) = Zmp + (A - Z)m„ + + E^^ , 

where and E^^ are respectively the strong and elec- 
tromagnetic contributions to the binding energy. The 
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Bethe-Weizacker formula gives that 

^BM = 98.25^^^^a,,, MeV. 



(26) 



If we decompose nip and nin as (see Gasser and 
Leutwyler, 1982) m(p „) = U3 + &(u,d)"i-u + ^(d,u)TOd + 
i3(p jj)Q:^„ where U3 is the pure QCD approximation of 
the nucleon mass {b^, and pj/ws being pure num- 
bers), it reduces to 

m{A,Z) = {Au3 + E^) (27) 
+ {Zbu + iV6d)TOu + {Zbd + Nbn)md 

ZBp + NBn + 98.25 MeV^ , 

with N = A — Z , the neutron number. This depends 
on our understanding of the description of the nucleus 
and can be more sophisticated. For an atom, one would 
have to add the contribution of the electrons. Zrric. The 
form (27) depends on strong, weak and electromagnetic 
quantities. The numerical coefficients B(^n,p) are given 
explicitly by (Gasser and Leutwiller, 1982) 

Bpa^^ =0.63 MeV B^a^^ = -OASMeV. (28) 

It follows that it is in general difficult to disentangle the 

effect of each parameter and compare the different meth- 
ods. For instance comparing the constraint on fi obtained 
from electromagnetic methods to the constraints on 
and Gp from nuclear methods requires to have some the- 
oretical input such as a theory to explain the fermion 
masses. Moreover, most of the theoretical models pre- 
dict a variation of the coupling constants from which one 
has to infer the variation of /x etc... 

For macroscopic bodies, the mass has also a negative 
contribution 



AmiG) = 



G fpif}p{f) 



d¥dV 



(29) 



from the gravitational binding energy. As a conclusion, 
from (27) and (29), we expect the mass to depend on all 
the coupling constant, 777,(0:^,^;, a.^^,, a.^, a.,.. ...). 

This has a profound consequence concerning the mo- 
tion of any body. Let a be any fundamental constant, 
assumed to be a scalar function and having a time vari- 
ation of cosmological origin so that in the privileged cos- 
mological rest-frame it is given by a{t). A body of mass 
m moving at velocity v will experience an anomalous ac- 
celeration 



6a = 



1 dmv dv 
m dt dt 



9 In TO 
da 



-av. 



(30) 



Now, in the rest-frame the body, a has a spatial depen- 
dence a[{t' + v.f' /c^)/ \/\ — f^/c^] so that, as long as 
<C c, Va = {a/c^)v. The anomalous acceleration can 
thus be rewritten as 



5a 



a 5mc? 
TO Sa 



Vina. 



(31) 



In the most general case, for non-relativistically moving 
body. 



, a Smc^ \ ( Va d v 
(5a = - ( — 1 -h — ^ 



(32) 



It reduces to Eq. (30) in the appropriate limit and the 
additional gradient term will be produced by local matter 
sources. This anomalous acceleration is generated by the 
change in the (electromagnetic, gravitational,...) binding 
energy (Dicke, 1964; Dicke, 1969; Haugan, 1979; Eardley, 
1979; Nordtvedt, 1990). Besides, the a-dependence is a 
priori composition-dependent (see e.g. Eq. 27). As a con- 
sequence, any variation of the fundamental constants will 
entail a violation of the universality of free fall: the total 
mass of the body being space dependent, an anomalous 
force appears if energy is to be conserved. The variation 
of the constants, deviation from general relativity and 
violation of the weak equivalence principle are in general 
expected together, e.g. if there exists a new interaction 
mediated by a massless scalar field. 

Gravitational methods include the constraints that can 
be derived from the test of the theory of gravity such as 
the test of the universality of free fall, the motion of the 
planets in the Solar system, stellar and galactic evolu- 
tions. They arc based on the comparison of two time 
scales, the first (gravitational time) dictated by grav- 
ity (ephemeris, stellar ages,. . . ) and the second (atomic 
time) is determined by any system not determined by 
gravity (e.g. atomic clocks,. . . ) (Canuto and Goldman, 
1982). For instance planet ranging, neutron star bina- 
ries observations, primordial nucleosynthesis and palcon- 
tological data allow to constraint the relative variation of 
G respectively to a level of IQ-^^ _ jq-ii^ ^Q-x^ _ 

Attacking the full general problem is a hazardous and 
dangerous task so that we will first describe the con- 
straints obtained in the literature by focusing on the fine 
structure constant and the gravitational constant and 
we will then extend to some other (less studied) com- 
binations of the constants. Another and complementary 
approach is to predict the mutual variations of different 
constants in a given theoretical model (see Section VI). 



III. FINE STRUCTURE CONSTANT 
A. Geological constraints 

1. The Oklo phenomenon 

Oklo is a prehistoric natural fission reactor that oper- 
ated about 2 X 10^ yr ago during (2.3 ± 0.7) x 10^ yr 
in the Oklo uranium mine in Gabon. This phenomenon 
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was discovered by the French Commissariat a I'Encrgie 
Atomique in 1972 (see Naudet, 1974, Maurette, 1976 and 
Petrov, 1977 for early studies and Naudet, 2000 for a 
general review) while monitoring for uranium ores. Two 
biUion years ago, uranium was naturally enriched (due 
to the difference of decay rate between ^35^ and ^^^U) 
and 235^ represented about 3.68% of the total uranium 
(compared with 0.72% today). Besides, in Oklo the con- 
centration of neutron absorbers which prevent the neu- 
trons from being available for the chain fission was low; 
water played the role of moderator and slowed down fast 
neutrons so that they can interact with other ^35^ and 
the reactor was large enough so that the neutrons did not 
escape faster than they were produced. 

Prom isotopic abundances of the yields, one can ex- 
tract informations about the nuclear reactions at the time 
the reactor was operational and reconstruct the reaction 
rates at that time. One of the key quantity measured is 
the ratio gf^Sm/g|^Sm of two light isotopes of samarium 
which are not fission products. This ratio of order of 0.9 
in normal samarium, is about 0.02 in Oklo ores. This low 
value is interpreted by the depletion of Sm by thermal 
neutrons to which it was exposed while the reactor was 
active. 

Shlyakhter (1976) pointed out that the capture cross 
section of thermal neutron by g^^Sm 



n + Sm 



6rSm + 7 



(33) 



is dominated by a capture resonance of a neutron of en- 
ergy of about 0.1 eV. The existence of this resonance is a 
consequence of an almost cancellation between the elec- 
tromagnetic repulsive force and the strong interaction. 

To obtain a constraint, one first needs to measure the 
neutron capture cross section of gt^Sm at the time of the 
reaction and to relate it to the energy of the resonance. 
One has finally to translate the constraint on the varia- 
tion of this energy on a constraint on the time variation 
of the considered constant. 

The cross section of the neutron capture (33) is 
strongly dependent on the energy of a resonance at 
Er = 97.3 meV and is well described by the Breit-Wigner 
formula 



(34) 



where go = (2J + l)(2s + 1)-^{2I + 1)-^ is a statistical 
factor which depends on the spin of the incident neutron 
,s = 1/2, of the target nucleus / and of the compound 
nucleus J; for the reaction (33), we have go = 9/16. The 
total width F = Fn + F-y is the sum of the neutron par- 
tial width Fn = 0.533 meV (at Er) and of the radiative 
partial width F^ = 60.5 meV. 

The effective absorption cross section is defined by 



1 2 f PTf c^^/^B^ 

a{Er,T) = / a^^^^){E)J — ——-VEdE 

(35) 

where the velocity vq = 2200 m-s~^ corresponds to an 
energy Eo — 25.3 meV and the effective neutron flux is 
similarly given by 



[2E e-^/'^B^ 



\/EdE. 



(36) 



The samples of the Oklo reactors were exposed 

(Naudet, 1974) to an integrated effective fluence J (j)dt 
of about 10^^ neutron-cm~^ = 1 kb~^. It implies that 
any process with a cross section smaller than 1 kb can 
be neglected in the computation of the abundances; this 
includes neutron capture by gl^Sm and gl^Sm. On the 
other hand, the fission of H^U, the capture of neutron 
by go^Nd and by g2^Sm with respective cross sections 
(T5 ~ 0.6 kb, ai43 ~ 0.3 kb and (7149 > 70 kb are the domi- 
nant processes. It follows that the equations of evolution 
for the number densities Nnj, Ni4^, N14C, and A'235 of 
62^Sm, 62*Sm, gl^Sm and gl^U are (Damour and Dyson, 
1996; Fujn et al, 2000) 



djVl47 

dt 

dNus 
dt 

dNi49 

dt 

diV235 

dt 



-ai4,7(f>Ni4,f + a-/235^-^235 (37) 

(T147M47 (38) 

-ai4g4>Nl49 + 17/2350^235 (39) 

-^5*^235 (40) 



where the system has to be closed by using a modified 

absorption cross section = (75(1 — C) (see references 
in Damour and Dyson, 1996). This system can be inte- 
grated under the assumption that the cross sections are 
constant and the result compared with the natural abun- 
dances of the samarium to extract the value of (T149 at 
the time of the reaction. Shlyakhter (1976) first claimed 
that (7149 = 55 ±8 kb (at cited by Dyson, 1978). Damour 
and Dyson (1996) re-analized this result and found that 
57 kb < 0-149 < 93 kb. Fujn et al. (2000) found that 
0-149 = 91 ± 6 kb. 

By comparing this measurements to the current value 
of the cross section and using (35) one can transform it 
into a constraint on the variation of the resonance energy. 
This step requires to estimate the neutron temperature. 
It can be obtained by using informations from the abun- 
dances of other isotopes such as lutetium and gadolinium. 
Shlyakhter (1976) deduced that \AEr\ < 20meV but as- 
sumed the much too low temperature of T = 20° C. 
Dyson and Damour (1996) allowed the temperature to 
vary between 180° C and 700° C and deduced the con- 
servative bound — 120meV < AEr < 90meV and Fujii 
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et al. (2000) obtained two branches, the first compatible 
with a null variation /S.Er = 9 ± 11 meV and the second 
indicating a non-zero effect AE^. = — 97 ± 8 meV both 
for T = 200 - 400° C and argued that the first branch 
was favored. 

Damour and Dyson (1996) related the variation of 

to the fine structure constant by taking into account 
that the radiative capture of the neutron by gl^Sm cor- 
responds to the existence of an excited quantum state 
g2°Sm (so that Er = E*r^Q — Eug — nin) and by assuming 
that the nuclear energy is independent of . It follows 
that the variation of a^.^, can be related to the difference 
of the Coulomb energy of these two states. The com- 
putation of this latter quantity is difficult and requires 
to be related to the mean-square radii of the protons in 
the isotopes of samarium and Damour and Dyson (1996) 
showed that the Bethe-Weizacker formula (26) overesti- 
mates by about a factor the 2 the qiej^ -sensitivity to the 
resonance energy. It follows from this analysis that 



AEr 

"EMX-^^-l-lMeV, 



(41) 



which, once combined with the constraint on AE^, im- 
plies 

-0.9 X 10-^ < Aa^^Ja^^^ < 1.2 x 10"^ (42) 

corresponding to the range —6.7 x 10~^^yr~^ < 
•^em/q^em < 5-0 ^ 10~^^yr~-^ if a^^ is assumed constant. 
This tight constraint arises from the large amplification 
between the resonance energy (~ 0.1 eV) and the sensi- 
tivity (~ 1 MeV). Fujii et al. (2000) re-analyzed the data 
and included data concerning gadolinium and found the 
favored result Q:em/q:em = (— 0.2±0.8)xl0~^^yr~^ which 
corresponds to 



Aa^/ttEM = (-0.36 ±1.44) X 10- 



(43) 



and another branch ciem/'^em = (4.9 ±0.4) x 10~^'^yr~-^. 
The first bound is favored given the constraint on the 
temperature of the reactor. Nevertheless, the non-zero 
result cannot be eliminated, even using results from 
gadolinium abundances (Fujii, 2002). Note however that 
spliting the analysis in two branches seems to be at odd 
with the aim of obtaining a constraint. Olive et al. (2002) 
refined the analysis and confirmed the previous results. 

Earlier studies include the original work by Shlyakhter 
(1976) who found that |q;em/q^e 
sponding to 



< 10 ^^yr ^ corre- 



|AaEM/aEMl <l-8xl0-^ 



(44) 



In fact he stated that the variation of the strong interac- 
tion coupling constant was given by Ag^/g^ ~ AEr/Vo 
where Vq ~ 50 MeV is the depth of a square potential 
well. Arguing that the Coulomb force increases the av- 
erage inter-nuclear distance by about 2.5% for A ~ 150, 



EM 



20Ag^/g^, leading to 



he concluded that Aa^^j^/a^ 

IcKem/q^emI < 10~"^'^yr~^. Irvine (1983a, b) quoted the 
bound Id 



/q^emI < 5 X 10~-^'^yr~^. The analysis of Sis- 
terna and Vucetich (1990) used, according to Damour 
and Dyson (1996) an ill-motivated finite-temperature de- 
scription of the excited state of the compound nucleus. 
Most of the studies focus on the effect of the fine structure 
constant mainly because the effects of its variation can 
be well controlled but, one would also have to take the ef- 
fect of the variation of the Fermi constant, or identically 
a^, (see Section V.A). Horvath and Vucetich (1988) in- 
terpreted the results from Oklo in terms of null-redshift 
experiments. 



2. a- decay 

The fact that a-decay can be used to put constraints 
on the time variation of the fine structure constant was 
pointed out by Wilkinson (1958) and then revived by 
Dyson (1972, 1973). The main idea is to extract the 
-dependence of the decay rate and to use geological 
samples to bound its time variation. 

The decay rate. A, of the a-decay of a nucleus of 
charge Z and atomic number A 



He 



(45) 



is governed by the penetration of the Coulomb barrier 
described by the Gamow theory and well approximated 
by 



,c/v 



(46) 



where v is the escape velocity of the a particle and where 

A is a function that depends slowly on a^^^ and v. It 
follows that the variation of the decay rate with respect 
to the fine structure constant is well approximated by 



da„,. V 



^ 1 dlnAg \ 
2dlnaEM/ 



(47) 



where AE = 2mv^ is the decay energy. Consider- 
ing that the total energy is the sum of the nuclear en- 
ergy E^^^ and of the Coulomb energy Ej^^/80 MeV ~ 
Z{Z — l)A~^/^ajj,j, and that the former does not depend 
on a„ 



one deduces that 
dlnAS / AE 



dln( 



0.6 MeV 



fiAZ) 



(48) 



with f{A,Z) = [(Z + 2)(Z + l)(A + 4)-V3 
-Z{Z - l)A-i/3] . It follows that the sensitivity of the 
decay rate on the fine structure constant is given by 

dlnA 



s = 



din, 



~ 47r.^-a^ 



10 



This result can be qualitatively understood since an in- 
crease of Qfj,^ induces an increase in the height of the 
Coulomb barrier at the nuclear surface while the depth 
of the nuclear potential below the top remains the same. 
It follows that the a particle escapes with greater en- 
ergy but at the same energy below the top of the barrier. 
Since the barrier becomes thinner at a given energy be- 
low its top, the penetrability increases. This computation 
indeed neglects the effect of a variation of on the nu- 
cleus that can be estimated to be dilated by about 1% if 
a^^ increases by 1%. 

Wilkinson (1958) considered the most favorable a- 
decay reaction which is the decay of 92^ U 

^fU^^^^Th + ^He (50) 

for which /\E ^ 4.27 McV (s ~ 540). By comparing the 
geological dating of the Earth by different methods, he 
concluded that the decay constant A of ^^^U, ^35^ ^nd 
232rjij^ have not changed by more than a factor 3 or 4 
during the last 3 — 4 x 10^ years from which it follows 
that laBivi/o^EMl < 2 X 10~^^yr~^ and thus 

|Aa™/aEMl <8x 10-3. (51) 

This bound is very rough but it agrees with Oklo on 
comparable time scale. This constraint was revised by 
Dyson (1972) who claimed that the decay rate has not 
changed by more than 20%, during the past 2 x 10^ years, 
which implies 

lAa^M/^KMl <4xl0-^ (52) 

These data were recently revisited by Olive et al. (2002). 
Using laboratory and meteoric data for -'^^^Sm (AE ~ 
2.31 MeV, s ~ 770) for which AA/A was estimated to be 
of order 7.5 x 10"^ they concluded that 

|Aa™/aKMl<10-^ (53) 



3. Spontaneous fission 

a-emitting nuclei are classified into four gcnerically in- 
dependent decay series (the thorium, neptunium, ura- 
nium and actinium series). The uranium series is the 
longest known series. It begins with gf^U, passes a sec- 
ond time through Z = 92 Q^'^V) as a consequence of an 
a-0-decay and then passes by five a-decays and finishes 
by an a-/3-(3-decay to end with |2^Pb. The longest lived 
member is gl^U with a half-life of 4.47 x 10^ yr, which 
four orders of magnitude larger than the second longest 
lived elements. g2*U thus determines the time scale of 
the whole series. 

The expression of the lifetime in the case of sponta- 
neous fission can be obtained from Gamow theory of a- 
decay by replacing the charge Z by the product of the 
charges of the two fission products. 



Gold (1968) studied the fission of 92^!! with a decay 
time of 7 X 10"^^ yr"^. He obtained a sensitivity (49) 
of s = 120. Ancient rock samples allow to conclude, 
after comparison of rock samples dated by potassium- 
argon and rubidium-strontium, that the decay time of 
g2^U has not varied by more than 10% in the last 2 x 
lO'' yr. Indeed, the main imccrtainty comes from the 
dating of the rock. Gold (1968) concluded on that basis 
that 

|Aa™/"BMl <4.66x 10-^ (54) 

which corresponds to |aj,^,/aj,^, | < 2.3 x lO-^^yr"^ if 
one assumes that a^^ is constant. This bound is indeed 
comparable, in order of magnitude, to the one obtained 
by a-decay data. 

Chitre and Pal (1968) compared the uranium-lead and 
potassium-argon dating methods respectively governed 
by a- and (3- decay to date stony meteoric samples. Both 
methods have different -dependence (see below) and 
they concluded that 

|Aa,M/aBMl<(l-5)xlO-^ (55) 

Dyson (1972) argued on similar basis that the decay rate 
of gf^U has not varied by more than 10% in the past 
2 X lO'-^ yr so that 

|Aa,,/a,,| <10-3. (56) 



4. /3-decay 

Dicke (1959) stressed that the comparison of the 

rubidium-strontium and potassium-argon dating meth- 
ods to uranium and thorium rates constrains the varia- 
tion of a^^ . He concluded that there was no evidence to 
rule out a time variation of the /3-decay rate. 

Peres (1968) discussed qualitatively the effect of a fine 
structure constant increasing with time arguing that the 
nuclei chart would have then been very different in the 
past since the stable heavy element would have had N/Z 
ratios much closer to unity (because the deviation from 
unity is mainly due to the electrostatic repulsion between 
protons). For instance 238 y would be unstable against 
double /?-decay to ^^^Pu. One of its arguments to claim 
that a^.^, has almost not varied lies in the fact that ^''^Pb 
existed in the past as ^''^Rn, which is a gas, so that the 
lead ores on Earth would be uniformly distributed. 

As long as long-lived isotopes are concerned for which 
the decay energy AE is small, we can use a non- 
relativistic approximation for the decay rate 

A = A± (AS)P± (57) 

respectively for /3~-decay and electron capture. A± are 
functions that depend smoothly on and which can 
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thus be considered constant, p+ = £ + 3 and p_ = 2£ + 2 
are the degrees of forbidckunic^ss of the transition. For 
high-Z nuclei with small dc^icu- energy AE, the exponent 
p becomes p = 2 + i^l — a^^^Z^ and is independent of i. 
It follows that the sensitivity (49) becomes 



lAa^M/aEMl <9x 10- 



(63) 



dlnAE 



din 



(58) 



The second factor can be estimated exactly as in 
Eq. (48) for a-decay but with f{A,Z) = ±{2Z + 
l)A~^/^[0.6MeV/Ai<^], the — , + signs corresponding re- 
spectively to /3-decay and electron capture. 

The laboratory determined decay rates of rubidium to 
strontium by /3-decay 



lyRb 



||Sr + De + e- 



and to potassium to argon by electron capture 



(59) 



(60) 



are respectively 1.41 X 10 ^^yr ^ and 4.72 x 10 yr ^. 
The decay energies are respectively AE = 0.275 MeV 
and AE = 1.31 MeV so that .s ~ -180 and s ~ -30. 
Peebles and Dicke (1962) compared these laboratories de- 
termined values with their abundances in rock samples 
after dating by uranium-lead method and with meteorite 
data (dated by uranium-lead and lead- lead). They con- 
cluded that the variation of a^y^ with cannot be ruled 
out by comparison to meteorite data. Later, Yahil (1975) 
used the concordance of the K-Ar and Rb-Sr geochemical 
ages to put the limit 



< 1.2 



(61) 



over the past 10"'^° yr. 

The case of the decay of osmium to rhenium by electron 
emission 



Re 



Os + i/e + e" 



(62) 



was first considered by Peebles and Dicke (1962). They 

noted that the very small value of its decay energy 
AE ~ 2.5 keV makes it a very sensitive indicator of 
the variation of Q!em- ^'^^ case p ~ 2.8 so that 
s ~ -18000. It follows that a change of about 10"^% 
of will induce a change in the decay energy of order 
of the keV, that is of the order of the decay energy itself. 
With a time decreasing a^^^^ , the decay rate of rhenium 
will have slowed down and then osmium will have become 
unstable. Peebles and Dicke (1962) did not have reliable 
laboratory determination of the decay rate to put any 
constraint. Dyson (1967) compared the isotopic analysis 
of molybdenite ores, the isotopic analysis of 14 iron me- 
teorites and laboratory measurements of the decay rate. 
Assuming that the variation of the decay energy comes 
entirely from the variation of a^y^ , he concluded that 



during the past 3 x 10^ years. In a re-analysis (Dyson, 
1972) he concluded that the rhenium decay-rate did not 
change by more than 10% in the past 10® years so that 



< 5 X 10^^ 



(64) 



Using a better determination of the decay rate of yf^Re 
based on the growth of ^^^Os over a 4-year period into 
a large source of osmium free rhenium, Lindner et al. 
(1986) deduced that 



Aa™/aEM = (-4.5±9)xlO- 



(65) 



over a 4.5 x 10® yr period. This was rcccnlty updated 
(Olive et al., 2002) to take into account the improvements 
in the analysis of the meteorite data which now show that 
the half-life has not varied by more than 0.5% in the past 
4.6 Gyr (i.e. a redshift of about 0.45). This implies that 



lAaEM/^EMl < 3 x 10 ^. 



(66) 



We just reported the values of the decay rates as used 
at the time of the studies. One could want to update 
these constraints by using new results on the measure- 
ments on the decay rate,. . . . Even though, they will not, 
in general, be competitive with the bounds obtained by 
other methods. These results can also be altered if the 
neutrinos are massive. 



5. Conclusion 

All the geological studic;s are on time scales of order of 
the age of the Earth (typically z ^ 0.1 — 0.15 depending 
on the values of the cosmological parameters) . 

The Oklo data are probably the most powerful geo- 
chemical data to study the variation of the fine structure 
constant but one has to understand and to model care- 
fully the correlations of the variation of q:^^, and as 
well as the effect of /i (see the recent study by Olive et 
al., 2002). This difiicult but necessary task remains to 
be done. 

The /3-decay results depend on the combination 
Q:f,.,aL and have the advantage not to depend on G. 

EM W or' 

They may be considered more as historical investigations 
than as competitive methods to constraint the variation 
of the fine structure constant, especially in view of the 
Oklo results. The dependence and use of this method on 
ttg was studied by Broulik and Trefil (1971) and Davies 
(1972) (see section V.B). 



B. Atomic spectra 

The previous bounds on the fine structure constant as- 
sume that other constants like the Fermi constant do not 
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vary. The use of atomic spectra may offer cleaner tests 
since we expect them to depend mainly on combinations 
of Q!em) M and Qp. 

We start by recalling some basics concerning atomic 
spectra in order to desribe the modelling of the spec- 
tra of many-electron systems which is of great use while 
studying quasar absorption spectra. We then focus on 
laboratory experiments and the results from quasar ab- 
sorption spectra. 



1. aj5^^ -dependence of atomic spectra 

As an example, let us briefly recall the spectrum of 
the hydrogen atom (see e.g. Cohen- Tannoudji et al, 
1977). As long as we neglect the effect of the spins and we 
work in the non-relativistic approximation, the spectrum 
is simply obtained by solving the Schrodinger equation 
with Hamiltonian 



Ho 



2mo 



(67) 



the eigenfunctions of which is of the form ipnim = 
Rn{r)Yim{0 , (t>) where n is the principal quantum num- 
ber. This solution has an energy 



Ei 



Er, = -^ll- — 



(68) 



independently of the quantum numbers I and m satisfy- 
ing < I < n, \m\ < I. It follows that there are states 
with the same energy. The spectroscopic nomenclature 
refers to a given energy level by the principal quantum 
number and a letter designing the quantum number I 
{s,p, d,f,g... respectively for / = 0, 1, 2, 3, 4 . . .). 

This analysis neglects rclativistic effects which arc ex- 
pected to be typically of order a'^^ (since in the Bohr 
model, v/c = a^^ for the orbit n = 1), to give the fine 
structure of the spectrum. The derivation of this fine 
structure spectrum requires to solve the Dirac equation 
for a particle in a potential —q'^/r and then to develop 
the solution in the non-relativist limit. Here, we simply 
use a perturbative approach in which the Hamiltonian of 
the system is expanded in v/c as 



H = Ho + W 



(69) 



where the corrective term W has different contributions. 
The spin-orbit interaction is described by 



W = 



he 

EM L S 



(70) 



and 2pi/2j where we have added in indices the total elec- 
tron angular moment quantum number J. The second 
correction arises from the (t;/c)^-relativistic terms and is 
of the form 



W 



(71) 



and it is easy to see that its amplitude is also of order 
W , ~ Hn. The third and last correction, known as 

rel EM ' 

the Darwin term, arises from the fact that in the Dirac 
equation the interaction between the electron and the 
Coulomb field is local. But, the non-relativist approx- 
imation leads to a non-local equation for the electron 
spinor that is sensitive to the field on a zone of order of 
the Compton wavelength centered in r. It follows that 



ttU q 



2^2 



■d{v). 



(72) 



The average in an atomic state is of order {W^) = 

nnV/i2mymO)\^ ~ me<^<u ~ "L^o- In con- 
elusion all the rclativistic corrections are of order q:^„ ^ 

EM 

(u/c)^. The energy of a fine structure level is 



EnU = irieC - 



Er 



rricC 



2n4 \J+l/2 



«EM + • • 



and is independent of the quantum number I. 

A much finer effect, referred to as hyperfine strueture, 
arises from the interaction between the spins of the elec- 
tron, S, and the proton, I. They are respectively associ- 
ated to the magnetic moments 



Ms 



qU S 



2me Ji 



Mi 



-9p 



qh I 



2mp h 



(73) 



Note that at this stage, the spectrum becomes depen- 
dent on the strong interaction via g-p (and via g/ in more 
general cases). This effect can be taken into account by 
adding the Hamiltonian 



+ 1 [3(Ms.n)(Mz.n) - M/.Ms] 



(74) 



where n is the unit vector pointing from the proton to the 
electron. The order of magnitude of this effect is typically 
e'^h^ / {menipC^r^) hence roughly 2000 times smaller than 
the effect of the spin-orbit coupling. It splits each fine 



Since r is of order of the Bohr radius, it follows that 

Q ^ '^EM'^^' ^^'^ splitting is indeed small: for in- 
stance, it is of order 4 x 10~^ eV between the levels 2^3/2 



*This 



valid 



to 



all 



order in and the Dirac equation directly gives E^u = 



1 + «EM (n-J-l/2 + ^(J + l/2)2-a|„)- 



-1/2 
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level in a series of hyperfine levels labelled by <E [| J — 
I\,I + J]. For instance for the level 2si/2 and 2pi/2, we 
have J =1/2 and F can take the two values and 1, for 
the level 2^3/2, J = 3/2 and F = 1 or = 2 etc. . . (see 
figure 4 for an example). This description neglects the 
quantum aspect of the electromagnetic field; one effects 
of the coupling of the atom to this field is to lift the 
degeneracy between the levels 2si/2 and 2pi/2- This is 
called the Lamb effect. 

In more complex situations, the computation of the 
spectrum of a given atom has to take all these effects 
into account but the solution of the Schrodinger equa- 
tion depends on the charge distribution and has to be 
performed numerically. 

The easiest generalization concerns hydrogen-like 
atoms of charge Z for which the spectrum can be ob- 
tained by replacing by (Ze)^ and rup by Amp. For 
an external electron in a many-electron atoms, the elec- 
tron density near the nucleus is given (see e.g. Dzuba 
et al., 1999a) by Z^Z/(n*ao)'^ where Za is the effective 
charge felt by the external electron outside the atom, 
n* an effective principal quantum number defined by 
En, = —EjZl/nl. It follows that the relativistic cor- 
rections to the energy level are given by 



u) = wo+qix + q2y 



(76) 



EM 



J +1/2 



'Y [ 4Z 



Such a formula does not take into account many-body 
effects and one expects in general a formula of the form 
AEn„i,j = En^Z^alJl/ J +1/2- C{Z,J,l)]/n,. Dzuba 
et al. (1999b) developed a method to compute the atomic 
spectra of many-clcctrons atoms including relativistic ef- 
fects. It is based on many-body perturbation theory 
(Dzuba et al., 1996) including electron-electron correla- 
tions and use a correlation-potential method for the atom 
(Dzuba et al, 1983). 

Laboratory measurements can provide these spectra 
but only for = a^^^ . In order to detect a variation 
of Qfj,^ , one needs to compute them for different values 
of a]5M- Dzuba et al. (1999a) describe the energy levels 
within one fine-structure multiplet as 



E = Eo + Qi 



(0) 



Q2 



(0) 



- 1 



(0) 



(0) 

. ^EA4 



(75) 



where Eq, Qi and Q2 describe the configuration center. 
The terms in L.S induce the spin-orbit coupling, sec- 
ond order spin-orbit interaction and the first order of the 
Breit interaction. Experimental data can be fitted to get 
Ki and K2 and then numerical simulations determine Qi 
and Q2. The result is conveniently written as 



with X = [a™/aW]' - 1 and y = [a.^/a^]" - 1- As 
an example, let us cite the result of Dzuba et al. (1999b) 
for Fe II 



6d 


J 


= 9/2 


to 


= 38458.9871 + 1394a; + 38y 




J 


= 7/2 


w 


= 38660.0494 -M632a;-h0y 


6/ 


J- 


= 11/2 


OJ 


= 41968.0642 + 1622a; + 3y 




J 


= 9/2 


OJ 


= 42114.8329+ 1772a; + 0y 




J 


= 7/2 


OJ 


= 42237.0500 + 1894a; + Oy 


6p 


J 


= 7/2 


OJ 


= 42658.2404 + 1398a; - Uy (77) 



with the frequency in cm~^ for transitions from the 
ground-state. An interesting case is Ni II (Dzuba et al., 
2001) which has largo relativistic effects of opposite signs 



2/ J = 7/2 w = 57080.373 - 300a; 
6d J = 5/2 uj = 57420.013 - 700a; 
6/ J = 5/2 oj = 58493.071 + 800x. 



(78) 



Such results are particularly useful to compare with spec- 
tra obtained from quasar absorption systems as e.g. in 
the analysis by Murphy et al. (2001c). 

In conclusion, the key point is that the spectra of atoms 
depend mainly on ji, aj^.^ and gp and contain terms both 
in Q^i^ and and that typically 

H = al^Ho + al^W,^^ + gpfJ,' al^W„^^ , (79) 

so that by comparing different kind of transitions in dif- 
ferent atoms there is hope to measure those constants 
despite the fact that plays a role via the nuclear mag- 
netic moment. We describe in the next section the lab- 
oratory oxporimonts and then turn to the measurement 
of quasar absorption spectra. 



2. Laboratory experiments 

Laboratory experiments are based on the comparison 
either of different atomic clocks or of atomic clock with 
ultra-stable oscillators. They are thus based only on the 
quantum mechanical theory of the atomic spectra. They 
also have the advantage to be more reliable and repro- 
ducible, thus allowing a better control of the systematics 
and a better statistics. Their evident drawback is their 
short time scales, fixed by the fractional stability of the 
least precise standards. This time scale is of order of a 
month to a year so that the obtained constraints are re- 
stricted to the instantaneous variation today, but it can 
be compensated by the extreme sensibility. They involve 
the comparison of either ultra-stable oscillators to differ- 
ent composition or of atomic clocks with different species. 
Solid resonators, electronic, fine structure and hyperfine 
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structure transitions respectively give access to -Roo/qJemj 
Roo, Rooal^ and gp/iRaoal^. 

Turneaure and Stein (1974) compared cesium atomic 
clocks with superconducting microwave cavities oscilla- 
tor. The frequency of the cavity-controlled oscillators 
was compared during 10 days that one of a cesium beam. 
The relative drift rate was (-0.4 ± 3.4) x 10""day"^ 
The dimensions of the cavity depends on the Bohr ra- 
dius of the atom while the cesium clock frequency de- 
pends on gpua'^^ (hyperfine transition). It follows that 



Iv t oc q-oLLC? so that 

Ce / cavity a^r^ 



- In (5pMaL)< 4.1x10-1 Vr- 



(80) 



Godone et al. (1993) compared the frequencies of ce- 
sium and magnesium atomic beams. The cesium clock, 
used to define the second in the SI system of units, is 
based on the hyperfine transition F — 3, mp = ^ F = 
4, mp = in the ground-state 6^Si/2 of ^'^^Ce with fre- 
quency given, at lowest order and neglecting relativistic 
and quantum electrodynamic corrections, by 



32cR 



3n3 



gill ~ 9.2 GHz, 



(81) 



where Zg the effective nuclear charge and gi the cesium 
nucleus gyromagnetic ratio. The magnesium clock is 
based on the frequency of the fine structure transition 
3pi 3po, Awj = in the meta-stable triplet of ^^Mg 



Hg 



It follows that 



6n^ 



601 GHz. 



dt 



X (l±10-2). 



The experiment led to the bound 



dt 



ln(5pM) 



< 5.4 X 10-13 



yr 



(82) 



(83) 



(84) 



after using the constraint dlii{gp/ gi)/dt < 5.5 x 
IQ-i^yr-^ (Demidov et ai, 1992). When combined with 
the astrophysical result by Wolfe et al. (1976) on the con- 
straint of gpiia^^ (see Section V.D) it is deduced that 



,/a,„| <2.7x lO-i'^yr 



(85) 



We note that relativistic corrections were neglected. 

Prestage et al. (1995) compared the rates of differ- 
ent atomic clocks based on hyperfine transitions in alkali 
atoms with different atomic numbers. The frequency of 
the hyperfine transition between / ± 1/2 states is given 
by (see e.g. Vanier and Audoin, 1989) 



8 



1 \ 2 ry^ 

2 / «em5/^- 



3 \ -i J n: 

(l-(5)(l-e)/ii?ooC. 



dn 



F,^i{a^^Z) 
(86) 



where z is the charge of the remaining ion once the va- 
lence electron has been removed and A„ = n — n*. The 
term {1 — 5) is the correction to the potential with respect 
to the Coulomb potential and (1 — e) a correction for the 
finite size of the nuclear magnetic dipole moment. It is 
estimated that 5 ~ 4% - 12% and e ~ 0.5%. F,ci(a^^^Z) 
is the Casimir relativistic contribution to the hyperfine 
structure and one takes advantage of the increasing im- 
portance of Fj-oi as the atomic number increases (see fig- 
ure 5). It follows that 



dt v„ 



V d In O' 



(87) 



where v-^ is the frequency of a H maser and when com- 
paring two alkali atoms 



d 2^ ^alkali 1 

dt 



d In Frel 



dlna„ 



dlni^, 



rel 



din 



(88) 



The comparison of different alkali clocks was performed 

and the comparison of Hg"*" ions with a cavity tuned H 
maser over a period of 140 days led to the conclusion that 



./a,„|<3.7xl0-"yr-i. 



(89) 



This method constrains in fact the variation of the quan- 
tity a^y^gp/gi. One delicate point is the evaluation of 
the correction function and the form used by Prestage et 
al. (1995) [i^rci ~ 1 + ll(^aj,„)V6 + ..] differs with 
the Is [Frei ~ 1 + 3{Za^^f/2 + ..] and 2s [frei ~ 
1 -|- n{Za^„Y /% + . . ] results for hydrogen like atoms 
(Breit, 1930). 

Sortais et al. (2001) compared a rubidium to a ce- 
sium clock over a period of 24 months and deduced that 

d\n{y^^/v^^)/dt = (1.9±3.1) x IQ-^^ yr-^, hence improv- 
ing the Tincertainty by a factor 20 relatively to Prestage 
et al. (1995). Assuming g^ constant, they deduced 



(4.2 ±6.9) X 10"'^ yr 



(90) 



if all the drift can be attributed to the Casimir relativistic 
correction F-^cx- 

All the results and characteristics of these experiments 
are summed up in table III. Recently, Braxmaier et al. 
(2001) proposed a new method to test the variability of 
ttj.^ and /i using electromagnetic resonators filled with 
a dielectric. The index of the dielectric depending on 
both a^,^ and /x, the comparison of two oscillators could 
lead to an accuracy of 4 x IQ-^^yr-^. Torgerson (2000) 
proposed to compare atom-stabilized optical frequency 
using an optical resonator. On an explicit example us- 
ing indium and thalium, it is argued that a precision of 
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•^em/q^em 10~^^/t, t being the time of the experiment, 
can be reached. 

Finally, let us note that similar techniques were used 
to test local Lorcntz invariancc (Lamoreaux et at, 1986, 
Chupp et al, 1989) and CPT symmetry (Bluhm et ai, 
2002). In the former case, the breakdown of local Lorentz 
invariancc would cause shifts in the energy levels of atoms 
and nuclei that depend on the orientation of the quantiza- 
tion axis of the state with respect to a universal velocity 
vector, and thus on the quantum numbers of the state. 



3. Astrophysical observations 

The observation of spectra of distant astrophysical ob- 
jects encodes information about the atomic energy levels 
at the position and time of emission. As long as one 
sticks to the non-relativistic approximation, the atomic 
transition energies are proportional to the Rydberg en- 
ergy and all transitions have the same aj.„-dependence, 
so that the variation will affect all the wavelengths by the 
same factor. Such a uniform shift of the spectra can not 
be distinguished from a Dopplcr effect due to the motion 
of the source or to the gravitational field where it sits. 

The idea is to compare different absorption lines from 
different species or equivalently the redshift associated 
with them. According to the lines compared one; can 
extract information about different combinations of the 
constants at the time of emission (see table I). 

While performing this kind of observations a number 
of problems and systematic effects have to be taken into 
account and controlled: 

1. Errors in the determination of laboratory wave- 
lengths to which the observations are compared, 

2. while comparing wavelengths from different atoms 

one has to take into account that they may be lo- 
cated in different regions of the cloud with different 
velocities and hence with different Doppler redshift. 

3. One has to ensure that there is no light blending. 

4. The differential isotopic saturation has to be con- 
trolled. Usually quasars absorption systems are 
expected to have lower heavy element abundances 
(Prochoska and Wolfe, 1996, 1997, 2000). The spa- 
tial inhomogeneity of these abundances may also 
play a role. 

5. Hyperfine splitting can induce a saturation similar 

to isotopic abundances. 

6. The variation of the velocity of the Earth during 
the integration of a quasar spectrum can induce 
differential Doppler shift. 



7. Atmospheric dispersion across the spectral direc- 
tion of the spectrograph slit can stretch the spec- 
trum. It was shown that this can only mimic a 
negative Aa^^^/a^^ (Murphy et al, 2001b). 

8. The prc;sc!ncc! of a magnetic field will shift the en- 
ergy levels by Zeeman effect. 

9. Temperature variations during the observation will 
change the air refractive index in the spectrograph. 

10. Instrumental effects such as variations of the intrin- 
sic instrument profile have to be controlled. 

The effect of these possible systematic errors are dis- 
cussed by Murphy et al. (2001b). In the particular case 
of the comparison of hydrogen and molecular lines, Wik- 
lind and Combes (1997) argued that the detection of the 
variation of /i was limited to A/Lt//i ~ 10~^. A possibility 
to reduce the systematics is to look at atoms having rel- 
ativistic corrections of different signs (see Section III.B) 
since the systematics are not expected, a priori, to simu- 
late the correlation of the shift of different lines of a mul- 
tiplet (see e.g. the example of Ni II Dzuba et al., 2001). 
Besides the systematics, statistical errors were important 
in early studies but have now enormously decreased. 

An efficient method is to observe fine-structure dou- 
blets for which 



2n3 



■ cm 



(91) 



Av being the frequency splitting between the two lines 
of the doublet and 9 the mean frequency (Bethe and 
Salpeter, 1977). It follows that Avjv oc c? and thus 



AlnAU/AlnAlo 
to give Aa^^ja 

E 



[1 -I- Aa^^ja^^ . It can be inverted 



as a function of AA and A as 



Aa, 



V A 



V A 



(92) 



As an example, it takes the following form for Si IV (Var- 
shalovich et al, 1996a) 



Ac 



{z) = 77.55 ( ^ ) - 0-5- (93) 



Since the observed wavelengths are redshifted as Aobs = 
Aem(l + z) it reduces to 



Aa„ 



{z) = 77.55 



Az 



(94) 



As a conclusion, by measuring the two wavelengths of 
the doublet and comparing to laboratory values, one can 
measure the time variation of the fine structure constant. 
This method has been applied to different systems and 
is the only one that gives a direct measurement of a^M • 
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SavcdofF (1956) was the first to realize that the fine and 
hyperfine structures can help to disentangle the redshift 
effect from a possible variation of a^^^ and Wilkinson 
(1958) pointed out that "the interpretation of redshift 
of spectral lines probably implies that atomic constants 
have not changed by less than 10~^ parts per year" . 

Savedoff (1956) used the data by Minkowski and Wol- 
son (1956) of the spectral lines of H, N II, O I, O II, 
Ne III and N V for the radio source Cygnus A of redshift 
z ^ 0.057. Using the data for the fine-structure doublet 
of N II and Ne III and assuming that the splitting was 
proportional to Q;^j^(l + z) led to 

Aa,„/a,„ = (1.8±1.6)xl0-^ (95) 

Bahcall and Salpeter (1965) used the fine structure split- 
ting of the O III and Ne III emission lines in the spectra 
of the quasi-stellar radio sources 3C 47 and 3C 147. Bah- 
call et al. (1967) used the observed fine structure of Si II 
and Si IV in the quasi-stellar radio sources 3C 191 to 
deduce that 

Aa.M/«EM = (-2±5)xlO-2 (96) 

at a redshift z = 1.95. Gamow (1967) criticized this mea- 
surements and suggested that the observed absorption 
lines were not associated with the quasi-stellar source but 
were instead produced in the intervening galaxies. But 
Bahcall et al. (1967) showed on the particular example 
of 3C 191 that the excited fine structure states of Si II 
were seen to be populated in the spectrum of this object 
and that the photon fluxes required to populate these 
states were orders of magnitude too high to be obtained 
in intervening galaxies. 

Bahcall and Schmidt (1967) then used the absorption 
lines of the O III multiplet of the spectra of five radio 
galaxies with redshift of order z ^ 0.2 to improve the 
former bound to 

Aa,„/a,„ = (l±2)xl0-3, (97) 

considering only statistical errors. 

Wolfe et al. (1976) studied the spectrum of 
AO 0235-M64, a BL Lac object with redshift z ~ 0.5. 
From the comparison of the hydrogen hyperfine fre- 
quency with the resonance line for Mg^, they obtained 
a constraint on gp^J^ct'^^ (see Section V.D). Prom the 
comparison with the Mg"*" fine structure separations they 
constrained gpUct^^ , and the Mg"*" fine structure doublet 
splitting gave 

|Aa™/«EMl <3xl0-2. (98) 

Potekhin and Varshalovich (1994) extended this method 
based on the absorption lines of alkali-like atoms and 
compared the wavelengths of a catalog of transitions 
2si/2 — 2p3/2 and 2si/2 ~ for a set of five ele- 

ments. The advantages of such a method are that (1) 



it is based on the measurement of the difference of wave- 
lengths which can be measured much more accurately 
than (broader) emission lines and (2) these transitions 
correspond to transitions from a single level and are thus 
not affected by differences in the radial velocity distribu- 
tions of different ions. They used data on 1414 absorption 
doublets of C IV, N V, O VI, Mg II, Al III and Si IV and 
obtained 

Aa,„/a,„ = (2.1 ± 2.3) X 10-3 (99) 

at ~ 3.2 and |dlnaj,„/dz| < 5.6x lO"' * between z = 0.2 
and z = 3.7 at 2a level. In these measurements Si IV, 
the most widely spaced doublet, is the most sensitive to 
a change in a-^-^ . The use of a large number of systems 
allows to reduce the statistical errors and to obtain a red- 
shift dependence after averaging over the celestial sphere. 
Note however that averaging on shells of constant red- 
shift implies that we average over a priori non-causally 
connected regions in which the value of the fine structure 
constant may a priori be different. This result was fur- 
ther constrained by Varshalovitch and Potekhin (1994) 
who extended the catalog to 1487 pairs of lines and got 

I Aa,„/a,„ I < 1.5x10-3 (100) 

at z ~ 3.2. It was also shown that the fine structure split- 
ting was the same in eight causally disconnected regions 
at 2 = 2.2 at a 3(7 level. 

Cowie and Songaila (1995) improved the previous anal- 
ysis to get 

Aa,„/a,„ = (-0.3 ± 1.9) x 10"^ (101) 

for quasars between z = 2.785 and z = 3.191. Var- 
shalovich et al. (1996a) used the fine-structure doublet 
of Si IV to get 

Aa™/«EM = (2±7)xlO-^ (102) 

at 2(7 for quasars between z = 2.8 and z = 3.1 (see also 
Varshalovich et al., 1996b). 

Varshalovich et al. (2000a) studied the doublet lines of 
Si IV, C IV and Ng II and focused on the fine-structure 
doublet of Si IV to get 

Aa^M/ttEM = (-4.5 ± 4.3[stat] ± 1.4[syst]) x 10"^ 

(103) 

for 2; = 2 — 4. An update of this analysis (Ivanchik et 
al., 1999) with 20 absorption systems between z = 2 and 

z = 3.2 gave 

Aa^M/a™ = (-3-3 ± 6.5[stat] ± 8[syst]) x 10"^ (104) 

Murphy et al. (2001d) used the same method with 21 
Si IV absorption system toward 8 quasars with redshift 
2; ~ 2 — 3 to get 
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Aa^/^EM = (-0.5 ± 1.3) X 10-^ (105) 
hence improving the previous constraint by a factor 3. 

Recently Dzuba et al. (1999a, b) and Webb et al. 
(1999) introduced a new method referred to as the many 
multiplet method in which one correlates the shift of the 
absorption lines of a set of multiplets of different ions. It 
is based on the parametrization (76) of the computation 
of atomic spectra. One advantage is that the correlation 
between different lines allows to reduce the systematics. 
An improvement is that one can compare the transitions 
from different ground-states and using ions with very dif- 
ferent atomic mass also increases the sensitivity because 
the difference between groimd-states rclativistic correc- 
tions can be very large and even of opposite sign (see the 
example of Ni II by Dzuba et al., 2001). 

Webb et al. (1999) analyzed one transition of the Mg II 
doublet and five Fe II transitions from three multiplets. 
The limit of accuracy of the method is set by the fre- 
quency interval between Mg II 2796 and Fe II 2383 which 
induces a fractional change of Aa^^^/a^^j ~ 10~^. Using 
the simulations by Dzuba et al. (1999a,b) it can be de- 
duced that a change in ol-^^^ induces a large change in the 
spectrum of Fe II and a small one for Mg II (the mag- 
nitude of the effect being mainly related to the atomic 
charge). The method is then to measure the shift of the 
Fe II spectrum with respect to the one of Mg II. This com- 
parison increases the sensitivity compared with methods 
using only alkali doublets. Using 30 absorption systems 
toward 17 quasars they obtained 

Aa^M/^EM = (-0-17 ± 0.39) x 10"^ (106) 
Aa™/aBM = (-1-88 ± 0.53) x IQ-^ (107) 

respectively for 0.6 < z < 1 and 1 < z < 1.6. There is 
no signal of a variation of for redshift smaller than 
1 but a 3.5(7 deviation for redshifts larger than 1 and 
particularly in the range z ~ 0.9 — 1.2. The summary 
of these measurements are depicted on figure 7. A pos- 
sible explanation is a variation of the isotopic ratio but 
the change of ^^Mg/^'*Mg would need to be substantial 
to explain the result (Murphy et al., 2001b). Calibra- 
tion effects can also be important since Fe II and Mg II 
lines are situated in different order of magnitude of the 
spectra. 

Murphy et al. (2001a) extended this technique of fit- 
ting of the absorption lines to the species Mg I, Mg II, 
Al II, Al III, Si II, Cr II, Fe II, Ni II and Zn II for 
49 absorption systems towards 28 quasars with redhsift 
2; ~ 0.5 — 3.5 and got 

H^a^^Ja^^^ = (-0.2 ± 0.3) x 10"^ (108) 
AQem/Qem = (-1-2 ± 0.3) X 10"^ (109) 

respectively for 0.5 < z < 1 and 1 < z < 1.8 at 4.1cr. The 
low redshift part is a re-analysis of the data by Webb et 



al. (1999). Over the whole sample {z = 0.5 — 1.8) it gives 
the constraint 

Aa™/«EM = (-0.7 ± 0.23) x 10"^ (110) 

Webb et al. (2001) re-analyzed their initial sample 
and included new optical QSO data to have 28 absorp- 
tion systems with redshift 2; = 0.5 — 1.8 plus 18 damped 
Lyman-a absorption systems towards 13 QSO plus 21 
Si IV absorption systems toward 13 QSO . The analysis 
used mainly the multiplets of Ni II, Cr II and Zn II and 
Mg I, Mg II, Al II, Al III and Fe II were also included. 
One improvement compared with the analysis by Webb 
et al. (1999) is that the "g" coefficient of Ni II, Cr II and 
Zn II in Eq. (76) vary both in magnitude and sign so that 
lines shift in opposite directions. The data were reduced 
to get 72 individual estimates of Aa^^Ja^^^^ spanning a 
large range of redshift. From the Fe II and Mg II sample 
they obtained 

Aa™/«EM = (-0.7 ± 0.23) x 10"^ (m) 

for ^ = 0.5 - 1.8 and from the Ni II, Cr II and Zn II they 
got 

Ac^em/c^em = ("0.76 ± 0.28) X 10"^ (112) 

for z = 1.8 — 3.5 at a 4a level. The fine-structure of Si IV 
gave 

Aa^/^EM = (-0.5 ± 1.3) X 10-5 (113) 
for z = 2 - 3. 

This series of results is of great importance since all 

other constraints are just upper bounds. Note that they 
are incompatible with both Oklo {z ~ 0.1) and mete- 
orites data (z 0.45) if the variation is linear with 
time. Such a non-zero detection, if confirmed, will have 
tremendous implications concerning our understanding 
of physics. Among the first questions that arise, it is 
interesting to test whether this variation is compatible 
with other bounds (e.g. test of the universality of free 
fall), to study the level of detection needed by the other 
experiments knowing the level of variation by Webb et 
al. (2001), to sort out the amplitude of the variation of 
the other constants and to be ensure that no systematic 
effects has been forgotten. For instance, the fact that 
Mg II and Fe II are a priori not in the same region of 
the cloud was not modelled; this could increase the er- 
rors even if it is difficult to think that it can mimic the 
observed variation of a^.^ . If one forgets the two points 
arising from HI 21 cm and molecular absorption systems 
(hollow squares in figure 7), the best fit of the data of 
figure 7 does not seem to favor today's value of the fine 
structure constant. This could indicate an unknown sys- 
tematic effect. Besides, if the variation of a-^^^ is mono- 
tonic then these observations seem to be incompatible 
with the Oklo results. 
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C. Cosmological constraints 

1. Cosmic microwave background 

The Cosmic Microwave Background Radi- 
ation (CMBR) is composed of the photons emitted at 
the time of the recombination of hydrogen and hehum 
when the universe was about 300,000 years old [see e.g. 
Hu and Dodelson (2002) or Durrer (2002) for recent re- 
views on CMBR physics]. This radiation is observed to 
be a black body with a temperature T = 2.723 K with 
small anisotropies of order of the /iK. The temperature 
fluctuation in a direction (^, ^p) is usually decomposed on 
a basis of spherical harmonics as 



5T 



m=+e 



{^,^) = J2 J2 aemYem{'&,<p)- (114) 



e TO=- 



Thc angular power spectrum miltipolc Ci = {\aim\'^) is 
the coefficient of the decomposition of the angular corre- 
lation function on Legendre polynomials. Given a model 
of structure formation and a set of cosmological parame- 
ters, this angular power spectrum can be computed and 
compared to observational data in order to constraint 
this set of parameters. 

Prior to recombination, the photons are tightly coupled 
to the electrons, after recombination they can be consid- 
ered mainly as free particles. Changing the fine structure 
constant modifies the strength of the electromagnetic in- 
teraction and thus the only effect on CMB anisotropies 
arises from the change in the differential optical depth of 
photons due to the Thomson scattering 



(115) 



which enters in the collision term of the Boltzmann equa- 
tion describing the evolution of the photon distribution 
function and where is the ionization fraction (i.e. the 
number density of free electrons with respect to their 
total number density rie)- The first dependence of the 
optical depth on the fine structure constant arises from 
the Thomson scattering cross-section given by 



3 m^c^ 



(116) 



and the scattering by free protons can be neglected since 
irie/mp ~ 5 X 10^''. The second, and more subtle depen- 
dence, comes from the ionization fraction. Recombina- 
tion proceeds via 2-photon emission from the 2s level or 
via the Ly-a photons which are redshifted out of the res- 
onance line (Peebles, 1968) because recombination to the 
ground state can be neglected since it leads to immediate 
reionization of another hydrogen atom by the emission of 
a Ly-a photon . Following Ma and Bertschinger (1995) 



and Peebles (1968) and taking into account only the re- 
combination of hydrogen, the equation of evolution of the 
ionization fraction takes the form 



dt 



■ C 



/? (1 - Xe) exp - 



TZupX^ 



Bn = —Ej/n^ is the energy of the nth hydrogen atomic 
level, P is the ionization coefficient, TZ the recombination 
coefficient, C the correction constant due to the redshift 
of Ly-a photons and to 2-photon decay and Up = is 
the number of proton. (3 is related to TZ by the principle 
of detailed balance so that 
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The recombination rate to all other excited levels is 



(117) 
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where (T„; is the ionization cross section for the (n, I) ex- 
cited level of hydrogen. The star indicates that the sum 
needs to be regularized and the a^^-, me-dependence 
of the ionization cross section is complicated to ex- 
tract. It can however be shown to behave as (T„; oc 
a->-2/(WSi). 

Finally, the factor C is given by 



C = 



1+XA2,,(1-Xe) 
1 + K{l3 + A2s){l-Xe) 



(118) 



where is the rate of decay of the 2s excited level 
to the ground state via 2 photons; it scales as rriea^^. 
The constant K is given in terms of the Ly-a phot;on 
Act = 167i'?i/(3mea^^c) by if = npX^/{8TrH) and scales 
as m~^a~^, . 

6 EM 

Changing q:j,^j will thus have two effects: first it 
changes the temperature at which the last scattering hap- 
pens and secondly it changes the residual ionization after 
recombination. Both effects influence the CMB temper- 
ature anisotropies [see Kaplinghat et al. (1999) and Bat- 
tye et al. (2001) for discussions]. The last scattering can 
roughly be determined by the maximum of the visibility 
function g = f exp(— r) which measures the differential 
probability for a photon to be scattered at a given red- 
shift. Increasing a^^j^^ shifts g to higher redshift at which 
the expansion rate is faster so that the temperature and 
Xe decrease more rapidly, resulting in a narrower g. This 
induces a shift of the Ce spcctriim to higher multipoles 
and an increase of the values of the C^. The first ef- 
fect can be understood by the fact that pushing the last 
scattering surface to a higher redshift leads to a smaller 
sound horizon at decoupling. The second effect results 
from a smaller Silk damping. 

Hannestad (1999) and then Kaplinghat et al. (1999) 
implemented these equations in a Boltzmann code, tak- 
ing into account only the recombination of hydrogen and 
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neglecting the one of helium, and showed that com- 
ing satellite experiments such as MAP^ and Planck^ 
should provide a constraint on a^.^ at recombination 
with a precision |Aj,^j/aj,^j | < 7 x 10~^'^yr^^. which cor- 
responds to a sensitivity jAa^^/aj^^l ~ 10^^ — 10^"^ 
at a redshift of about z ^ 1,000. Avelino et al. 
(2000) studied the dependence of the position of the 
first acoustic peak on a^M- Hannestad (1999) chose the 
underlying ACDM model (f2, f2b, A, /i, n, TVj^, r, a^.^^) = 
(1, 0.08, 0, 0.5, 1, 3, 0, ^^'^ performed a 8 parameters 
fit to determine to which precision the parameters can be 
extracted. Kaplinghat et al. (1999) worked with the pa- 
rameters (/i, J7b, A, iVjy, Fp, q;jj,^j). They showed that the 
precision on Aa^^^/a,^^ varies from 10~^ if the maximum 
observed CMB multipole is of order 500-1000 to 10"^ if 
one observes multipoles higher than 1500. 

Avelino et al. (2000) claim that BOOMERanG and 
MAXIMA data favor a value of a^^^ smaller by a few 
percents in the past (see also Martins et al, 2002) and 
Battye et al. (2001) showed that the fit to current CMB 
data are improved by allowing Aa^^ ^ and pointed 
out that the evidence of a variation of the fine structure 
constant can be thought of as favoring a delayed recom- 
bination model (assuming 0=1 and n = 1). Avelino et 
al. (2001) then performed a joint analysis of nucleosyn- 
thesis and CMB data and did not find any evidence for a 
variation of a^^ at one-sigma level at either epoch. They 
consider Ob and Aa^^^ as independent and the marginal- 
ization over one of the two parameters lead to 



-0.09 < Aoj^^j < 0.02 



(119) 



at 68% confidence level. Martins et al. (2002) concluded 
that MAP and Planck will allow to set respectively a 
2.2% and 0.4% constraint at la if all other parameters 
are marginalized. Landau et al. (2001) concluded from 
the study of BOOMERanG, MAXIMA and COBE data 
in spatially flat models with adiabatic primordial fluctu- 
ations that, at 2cr level. 



-0.14 < Aa^„ < 0.03. 



(120) 



All these works assume that only q:j,^j is varying but, 
as can been seen from Eqs. (114-118), one has to assume 
the constancy of the electron mass. Battye et al. (2001) 
show that the change in the fine structure constant and 
in the mass of the electron are degenerate according to 
^'^EM ~ 0.39Ame but that this degeneracy was broken 
for multipoles higher than 1500. The variation of the 
gravitational constant can also have similar effects on 
the CMB (Riazuelo and Uzan, 2002). All the works also 
assume the a^.^ -dependence of TZ to be negligible and 



Battye et al. (2001) checked that the helium recombina- 
tion was negligible in the range of Aa^,^ considered. 

In conclusion, strong constraints on the variation of 
a^.^, can be obtained from the CMB only if the cosmolog- 
ical parameters are independently known. This method 
is thus non competitive unless one has strong bounds on 
Ob and h (and the result will always be conditional to the 
model of structure formation) and assumptions about the 
variation of other constants such as the electron mass, 
gravitational constant are made. 



2. Nucleosynthesis 

The amount of '*He produced during the big bang nu- 
cleosynthesis is mainly determined by the neutron to pro- 
ton ratio at the freeze-out of the weak interactions that 

interconvcrt neutrons and protons. The result of Big 
Bang nucleosynthesis (BBN) thus depends on G, a^, 
a^^j and respectively through the expansion rate, the 
neutron to proton ratio, the neutron- proton mass differ- 
ence and the nuclear reaction rates, besides the standard 
parameters such as e.g. the number of neutrino fami- 
lies. The standard BBN scenario (sec e.g. Malaney, 1993, 
Reeves, 1994) proceeds in three main steps: 

1. for T > 1 MeV, (t < 1 s) a first stage during 
which the neutrons, protons, electrons, positrons 
an neutrinos are kept in statistical equilibrium by 
the (rapid) weak interaction 



n + e~' 



3 +Ve, n + Ve 



P + e , 



(121) 



As long as statistical equilibrium holds, the neutron 
to proton ratio is 



(n/p) = e 



(122) 



where Q = {rn^ - mp)c^ = 1.29 MeV. The abun- 
dance of the other light elements is given by (Kolb 
and Turner, 1993) 



A-l 



2(3A-5)/2^5/2 
3(A-l)/2 



^A-lyZyA-Z^B^/k^T^ (123) 



^http: //map. gsfc .nasa.gov/ 

®http : / / astro . estec . esa . nl/SA-general/Pro jects/Planck/ 



where gA is the number of degrees of freedom of the 
nucleus ^X, ton is the nucleon mass, t] the baryon- 
photon ratio and Ba = {Zmp + {A — Z)m^ — mA)c? 
the binding energy. 

2. Around T ~ 0.8 MeV (t 2 s), the weak inter- 
actions freeze out at a temperature Tf determined 
by the competition between the weak interaction 
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rates and the expansion rate of the universe and 
thus determined by r„(rf) ~ H{Ti) that is 



(124) 



where Gp is the Fermi constant and N^, the num- 
ber of rclativistic degrees of freedom at Tf. Be- 
low Tf , the number of neutrons and protons change 
only from the neutron /3-decay between Tf to Tn ~ 
0.1 MeV when p + n reactions proceed faster than 
their inverse dissociation. Tn is determined by de- 
manding that the relative number of photons with 
energy larger that the deuteron binding energy, 
i^D, is smaller than one, i.e. so that n^/rip ~ 
ex.p{En/Tf^) ~ 1. 

3. For 0.05 MeV< T < 0.6 MeV {3s < t < 6min), 
the synthesis of light elements occurs only by two- 
body reactions. This requires the deuteron to be 
synthetized {p + n ^ D) and the photon density 
must be low enough for the photo-dissociation to 
be negligible. This happens roughly when 



^~r?2exp(-^D/TN) 



(125) 



with ?7 ^ 3 X 10 The abundance of ^He by mass, 
Yp, is then well estimated by 



(n/p)^ 



1 -I- (n/p)N 



with 



{n/p)N = {n/p)fexp{-tf^/Tn) 



(126) 



(127) 



with tN oc G-'^/^T^^ and t'^ = 1.636G2(1 + 
3g\)ml/{2TT^), with qa ^ 1-26 being the ax- 
ial/vector coupling of the nucleon. Assuming 
that Ed (x a^, this gives a dependence ti<s/Tp cx 
G-V2a2G2 (see Section V.B). 

The helium abundance depends thus mainly on Q. 
Tf and Tn (and hence mainly on the neutron life- 
time, Tn) and the abundances of the other elements 
depends also on the nuclear reaction rates. 

The light element abundances are thus sensible to the 

freeze-out temperature, which depends on Gp, G, on the 
proton- neutron mass difference Q, and on the values of 
the binding energies Ba so that they mainly depend a^^ , 
a^, ctg, and the mass of the quarks. An increase in 
G or N^, results in a higher expansion rate and thus to 
an earlier freeze-out, i.e. a higher Tf. A decrease in Gp, 
corresponding to a longer neutron lifetime, leads to a 
decrease of the weak interaction rates and also results in 
a higher Tf . It inplies, assuming uncorrelated variations. 



that |AG/G| < 0.25 (see Section IV) and |AGp/Gp| < 
6 X 10-2 (see Section V.A). 

First, the radiative and Coulomb corrections for the 

weak reactions (121) have been computed by Dicus et al. 
(1982) and shown to have a very small influence on the 
abundances. 

The constraints on the variation of these quantities 
were first studied by Kolb et al. (1986) who calculated 
the dependence of primordial *He on G, Gp and Q. They 
studied the influence of independent changes of the for- 
mer parameters and showed that the helium abundance 
was mostly sensitive in the change in Q. Other abun- 
dances are less sensitive to the value of Q, mainly be- 
cause "^He has a larger binding energy; its abundances 
is less sensitive to the weak reaction rate and more to 
the parameters fixing the value of (n/p). To extract the 
constraint on the fine structure constant, one needs a 
particular model for the -dependence of Q. Kolb et 
al. (1986) decomposed Q as 



(128) 



where the first part represents the electromagnetic con- 
tribution and the second part corresponds to all non- 
electromagnetic contributions. Assuming that Qa and 
Q/3 are constant and that the electromagnetic con- 
tribution is the dominant part of Q, they deduce 
that Q/Qo — OisM/ctEM ^'^d thus that (n/p) ~ 
{n/p)o [l — qoTfa^f^ / a'^^] . To consider the effect of the 
dependent variation of G, Gp and a^j, , the time vari- 
ation of these constants was related to the time varia- 
tion of the volume of an internal space of characteris- 
tic size R for a 10-dimcnsional superstring model and 
Kaluza-Klein models (see Section VI for details on these 
models)''. They concluded that 



lAapM/^BMl < 10" 



(129) 



and showed that if one requires that the abundances of 
and ^He remains unchanged it is impossible to com- 
pensate the change in Q!p„ by a change in the baryon-to- 
photon ratio. Indeed, the result depends strongly on the 
hypothesis of the functional dependence. Khare (1986) 
then showed that the eff'ect of the extra-dimensions can 
be cancelled if the primordial neutrinos are degenerate. 
This approach was generalized by Vayonakis (1988) who 
considered the 10-dimensional limit of superstring and by 
Coley (1990) for the case of 5-dimensional Kaluza-Klein 
theory. 

Campbell and Olive (1995) kept track of the changes 
in Tf and Q separately and deduced that 



Their hypothesis on the variation of the Fermi constant are 
questionable, see Section V.A for details. 
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They used this to study the constraints on (see Sec- 
tion V.A). 

Bergstrom et al. (1999) extended the original work by 
Kolb et al. (1986) by considering other nuclei. They 
assumed the dependence of Q on a,^^ 



Q ~ (1.29 - 0.76Aa^„/a^„) MeV 



(131) 



that relies on a change of quark masses due to strong and 
electromagnetic energy binding. Since the abundances of 
other nuclei depend mostly on the weak interaction rates, 
they studied the dependence of the thermonuclear rates 
on a^M- 111 the non-relativistic limit, it is obtained as 
the thermal average of the cross section times the relative 
velocity times the number densities. The key point is that 
for charged particles the cross section takes the form 



E 



(132) 



where r]{E) arises from the Coulomb barrier and is given 
in terms of the charges and the reduced mass fj, of the 
two particles as 



(133) 



The factor S{E) has to be extrapolated from experimen- 
tal nuclear data which allows Bergstrom et al. (1999) to 
determine the Q!jj,„ -dependence of all the relevant reac- 
tion rates. Let us note that the -dependence of the 
reduced mass ^ and of S{E) were neglected; the latter 
one is polynomial in a-^^ (Fowler et at, 1975). Keep- 
ing all other constants fixed, assuming no exotic effects 
and taking a lifetime of 886.7 s for the neutron, it was 
deduced that 



lAa^M/^EMl <2x 10- 



(134) 



In the low range of ~ 1.8 x 10~^° the ''Li abundance 

does not depend strongly on and the one of "^Hc has 
to be used to constrain a^,^ . But it has to be noted that 
the observational status of the abundance of ^He is still 
a matter of debate and that the theoretical prediction of 
its variation with depends on the model-dependent 
ansatz (131). For the high range of ?7 ^ 5 x 10"^'^, the 
variation of ''Li with a^^^^ is rapid, due to the exponential 
Coulomb barrier and limits the variation of a^M • 

Nollet and Lopez (2002) pointed out that Eq. (132) 
does not contain all the q:jj,^j -dependence. They argue 
that (i) the factor S depends linearly on a^,^ , (ii) when 
a reaction produces two charged particles there should 
be an extra contribution arising from the fact that 
the particles need to escape the Coulomb potential, (iii) 
the reaction energies depend on a,^^ and (iv) radiative 



captures matrix elements are proportional to a-^^ . The 
most secure constraint arising from D/H measurements 
and combining with CMB data to determine fls gives 



Aa™/aEM = (3±7)xlO- 



(135) 



at la level. 

Ichikawa and Kawasaki (2002) included the effect of 
the quark mass and by considering a joint variation of 
the different couplings as it appears from a dilaton. Q 
then takes the form 



Q = a"EM^QCD + KVd - yu)v 



(136) 



where a and b are two parameters and yd, Uu the Yukawa 
couplings. The neutron lifetime then behaves as 



{l/vyl)f-\Q/m, 



(137) 



with / is a known function. Assuming that all the cou- 
plings vary due to the effect of a dilaton, such that the 
Higgs vacuum expectation value v remains fixed, they 
constrained the variation of this dilaton and deduced 

Aa,M /"em = (-2.24 ± 3.75) x 10"^ (138) 

In all the studies, one either assumes all other con- 
stants fixed or a functional dependence between them, 
as inspired from string theory. The bounds are of the 
same order of magnitude that the ones obtained from the 
CMB; they have the advantage to be at higher redshift 
but suffer from the drawback to be model-dependent. 



3. Conclusion 

Even if cosmological observations allow to test larger 
time scales, it is difficult to extract tight constraints on 
the variation of the fine structure constant from them. 

The CMB seems clean at first glance since the effect 
of the fine structure constant is well decoupled from the 
effect of the weak and strong coupling constants. Still, 
it is entangled with assumption on G. Besides, it was 
shown that degeneracy between some parameters exists 
and mainly between the fine structure constant, the elec- 
tron to proton mass ratio, the baryonic density and the 
dark energy equation of state (Huey et al., 2001). 

Nucleosynthesis is degenerate in the four fundamental 
coupling constants. In some specific models where the 
variation of these constants are linked it allows to con- 
straint them and definitively the helium abundance alone 
cannot constraint the fine structure constant. 



D. Equivalence Principle 

The equivalence principle is closely related to the de- 
velopment of the theory of gravity from Newton's the- 
ory to general relativity (see Will, 1993 and Will, 2001 
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for reviews). Its first aspect is the weak equivalence 
principle stating that the weight of a body is propor- 
tional to its mass or equivalently that the trajectory of 
any freely faUing body does not depend on its internal 
structure, mass and composition. Einstein formulated 
a stronger equivalence principle usually referred to as 
Einstein equivalence principle stating that (1) the weak 
equivalence principle holds, (2) any non-gravitational ex- 
periment is independent of the velocity of the laboratory 
rest-frame (local Lorentz invariance) and (3) of when an 
where it is performed (local position invariance). 

If the Einstein equivalence principle is valid then grav- 
ity can be described as the consequence of a curved space- 
time and is a metric theory of gravity, an example of 
which are general relativity and the Brans-Dicke (1961) 
theory. This statement is not a "theorem" but there are 
a lot of indications to back it up (see Will, 1993, 2001). 
Note that superstring theory violates the Einstein equiv- 
alence principle since it introduces additional fields (e.g. 
dilaton, moduli...) that have gravitational-strength cou- 
plings which violates of the weak equivalence principle. A 
time variation of a fundamental constant is in contradic- 
tion with Einstein equivalence principle since it violates 
the local position invariance. Dicke (1957, 1964) was 
probably the first to try to use the result of Eotvos et al. 
(1922) experiment to argue that the strong interaction 
constant was approximatively position independent. All 
new interactions that appear in the extension of standard 
physics implies extra scalar or vector fields and thus an 
expected violation of the weak equivalence principle, the 
only exception being metric theories such as the class of 
tensor-scalar theories of gravitation in which the dilaton 
couples universally to all fields and in which one can have 
a time variation of gravitational constant without a vio- 
lation of the weak equivalence principle (see e.g. Damour 
and Esposito-Farese, 1992). 

The difference in acceleration between two bodies of 
different composition can be measured in Eotvos- type ex- 
periments (Eotvos et al., 1922) in which the acceleration 
of various pairs of material in the Earth gravitational 
field are compared. The results of this kind of laboratory 
experiments are presented as bounds on the parameter r] 



\ai - a2\ 
\ai + a2\ 



(139) 



The most accurate constraints on t] are t] = (— 1.9 ± 
2.5) X lO"'^^ between beryllium and copper (Su et al., 
1994) and \r]\ < 5.5 x 10~^^ between Earth-core-like 
and Moon-mantle-like materials (Baessler et al, 1999). 
The Lunar Laser Ranging (experiment) gives the bound 
ri = (3.2 ± 4.6) X lO'^^ (Williams et al., 1996) and 
r] = (3.6 ± 4) X 10-13 (Miiller and Nordtvedt, 1998; 
Miiller et al., 1999). Note however that, as pointed by 
Nordtvedt (1988, 2001a), the LLR measurement are am- 
biguous since the Earth and the Moon have (i) a different 



fraction of gravitational self-energy and (ii) a difference 
of composition (the core of the Earth having a larger 
Fe/Ni ratio than the Moon). This makes this test sensi- 
titive both to self-gravity and to non-gravitational forms 
of energy. The experiment by Baessler et al, (1999) lifts 
the degeneracy by considering miniature "Earth" and 
"Moon". 

As explained in Section II. C, if the self-energy depends 
on position, the conservation of energy implies the exis- 
tence of an anomalous acceleration. In the more general 
case where the long range force is mediated by a scalar 
field (j), one has to determine the dependence mi{(j)) of 
the different particles. If it is different for neutron and 
proton, then the force will be composition dependent. 
At the Newtonian approximation, the interaction poten- 
tial between two particles is of the form (Damour and 
Esposito-Farese, 1992) 



mim2 



V{r) = -G (l + aiae-'-/^) 
with ai2 = /1/2 and /j defined as 
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where = SttG/Uc is the four dimensional Planck 

mass. The coefficient ai2 is thus not a fundamental con- 
stant and depends a priori on the chemical composition 
of the two test masses. It follows that 



V12 = 



./cxt|./l-/2| 
l + /ext(/l + /2)/2 



M4/e: 



mi 
m2 



(142) 



To set any constraint, one has to determine the 
functions fi{(j>), which can only be made in a model- 
dependent approach [see e.g. Damour (1996) for a dis- 
cussion of the information that can be extracted in a 
model-independent way]. For instance, if couples to a 
charge Q the additional potential is expected to be of the 
form 



y(r-)=-/g^e-A 



(143) 



with /q being a fundamental constant (/g > for scalar 

exchange and /g < for vector exchange). It follows 
that ai2 depends explicitly of the composition of the two 
bodies as 



0:12 
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Ml At2 
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where fii = mi/m,^ and £_q = .fq/Gmfj^. Their relative 
acceleration in an external field (fcxt is 



Aa 
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For instance, in the case of a fifth force induced by a dila- 
ton or string moduli, Damour and Polyakov (1994a, b) 
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showed that there are three charges B = N + Z, D = 
N — Z and E = Z(Z — \)B^I^ representing respectively 
the baryon number, the neutron excess and a term pro- 
portional to the nuclear Coulomb energy. The test of the 
equivalence principle results in an exclusion plot in the 
plane (^q,A) (see Figure 6). 

To illustrate the link between the variation of the 
constants and the tests of relativity, let us considered 
the string-inspired model developed by Damour and 
Polyakov (1994a, 1994b), in which the fine structure 
constant is given in terms of a function of the four di- 
mensional dilaton as a„,, = Bp^{4>). The QCD mass 
scale can be expressed in terms of the string mass scale, 
Ms ~ 3 X 10^^ GeV [see Section VLB for details and 
Eq. (276)]. In the chiral limit, the (Einstein- frame) 
hadron mass is proportional to the QCD mass scale so 
that, 



/hadron 



In- 



1\ 91n< 



^hadron 



(146) 



With the expected form lni?j?((/>) = —«(</>— (^mf' jl (see 
Section VLB), the factor of the r.h.s. of the previous 
equation is of order 40k (0 — (\)m)- The exchange of the 
scalar field excitation induces a deviation from general 
relativity characterized, at post-Newtonian level, by the 
Eddington parameters 



l-7Eaa-2(40K)2(.^0-.^. 

- 1 ~ (40k)2(,^o - <^m)V2. 



(147) 
(148) 



Besides, the violation of the universality of free fall is 
given by 7712 = — ^2 with 



^1 = (l-7Edd) 



B_' 

1 



-1-0.943 X 10"^ I - 



cd\ — 
M/ 1 

E 



(149) 



obtained from the expression (27) for the mass. In 
this expression, the third term is expected to dominate. 
We sec on this example that the variation of the con- 
stants, the violation of the equivalence principle and post- 
Newtonian deviation from general relativity have to be 
considered together. 

Similarly, in an effective 4-dimensional theory, the only 
consistent approach to make a Lagrangian parameter 
time dependent is to consider it as a field. The Klein- 
Gordon equation for this field {(j) + 'i^H cj) + rn^ (j) + . . . = 0) 
implies that ^ is damped bs (f) (x if its mass is much 
smaller than the Hubble scale. Thus, in order to be vary- 
ing during the last Hubble time, </> has to be very light 
with typical mass m ^ /fo ~ 10"^^ eV. This is analo- 
gous to the case of quintessence models (see Section V.E 
for details). As a consequence, (j) has to be very weakly 
coupled to the standard model fields. To illustrate this. 



Dvali and Zaldarriaga (2002) [followed by a re-analysis 
by Chiba and Khori (2001), Wetterich (2002)] expanded 
around it value today as 



«em=«bm(0) + A— + 



(150) 



from which it follows, from Webb et al. (2001) measure, 
that AA(/)/Af4 ^ 10"'^ during the last Hubble time. The 
change of the mass of the proton and of the neutron due 
to electromagnetic effects was obtained from Eqs. (27-28) 
but with neglecting the last term. The extra-Lagrangian 
for the field 6 is thus 



5L = \— {Bppp + Bnun) . 

IVI4 



(151) 



A test body composed of rin neutrons and rip protons will 
be characterized by a sensitivity 



fi = — (t'p-Bp + I'n-Bn) 

rriN 



(152) 



where (resp. Vp) is the ratio of neutrons (resp. 

protons) and where it has been assumed that nin ~ 
nip ~ Ton- Assuming* that I'n^^^ ~ 1/2 and us- 
ing that the compactness of the Moon-Earth system 
51n(TOEarth/™Moon)/91nQ!j,^, ^ 10"^, ouc gcts r]i2 ~ 
10~^A^. Dvali and Zaldarriaga (2002) obtained the same 
result by considering that Aun,p 6 x 10~^ — 10~^. This 
implies that A < 10~^ which is compatible with the vari- 
ation of a^^ if A0/M4 > 10"^ during the last Hubble 
period. 

From cosmological investigations one can show that 
(A^/M4)^ ~ {P4, + Pff,) / Ptota.\- If dominates the matter 
content of the universe, Ptotaii then A(f) ^ M4 so that 
A ^ 10"'' whereas if it is sub-dominant Acp <C M4 and 
A ^ lO"*". In conclusion 



10"^ < A < 10"^ 



(153) 



This explicits the tuning on the parameter A. 

An underlying approximation is that the ^-dependence 
arises only from the electromagnetic self-energy. But, in 
general, one would expect that the dominant contribu- 
tion to the hadron mass, the QCD contributions, also in- 
duces a (?i-dependence (as in the Damour and Polyakov, 
1994a, b approach). 

In conclusion, the test of the equivalence principle 
offers a very precise test of the variation of constants 
(Damour, 2001). The LLR constrain rj < 10'^^, i.e 
[ttEarth — flMoon | ^ 10^^^ cm.s^'^, implies that on the size 



'For copper Vp = 0.456, for uranium i/p = 0.385 and for lead 
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of the Earth orbit |Vlna^^J < lO"'^^ _ iQ-:i2 ^^-i 
tending this measurement to the Hubble size leads to the 
estimate Aag^/a^M ^ 10~^ — 10~^. This indicates that 
if the claim by Webb et al. (2001) is correct then it should 
induce a detectable violation of the equivalence princi- 
ple by coming experiments such as MICROSCOPE^ and 
grpgpio ^jjj ^^g^ respectively at the level 77 ^ 10"^^ 
and rj ~ 10^ Indeed, this is a rough estimate in which 
d^M is assumed to be constant, but this is also the con- 
clusion indicated by the result by Dvali and Zaldarriaga 
(2002) and Bekenstein (1982). 

Let us also note that this constraint has been discarded 
by a some models (sec Section VI. C) and particularly 
while claiming that a variation of of 10~^ was re- 
alistic (Sandvik et al., 2002; Barrow et al., 2001) [see 
however the recent discussion by Magueijo et al. (2002)] . 

IV. GRAVITATIONAL CONSTANT 

As pointed by Dickc and Peebles (1965), the impor- 
tance of gravitation on large scales is due to the short 
range of the strong and weak forces and to the fact that 
the electromagnetic force becomes weak because of the 
global neutrality of the matter. As they provide tests of 
the law of gravitation (planetary motions, light deflec- 
tion,. . . ), space science and cosmology also ofl^er tests of 
the constancy of the gravitational constant. 

Contrary to most of the other fundamental constants, 
as the precision of the measurements increased, the dis- 
parity between the measured values of G also increased. 
This led the CODATA" in 1998 to raise the relative un- 
certainty for G from 0.013% to 0.15% (Gundlach and 
Merkowitz, 2000). 

A. Paleontological and geophysical arguments 

Dickc (1964) stressed that the Earth is such a complex 
system that it would be difficult to use it as a source 
of evidence for or against the existence of a time varia- 
tion of the gravitational constant. He noted that among 
the direct effects, a weakening of the gravitational con- 
stant induces a variation of the Earth surface tempera- 
ture, an expansion of the Earth radius and a variation of 
the length of the day (Jordan, 1955, and then Murphy 
and Dicke, 1964; Hoyle, 1972). 



®http : / / sci2 . esa. int/Microscope/ 
^°http : //einstein . Stanford . edu/STEP/ 

"The COD ATA is the COmmittee on Data for Science and 
Technology, see http://www.codata.org/. 



1. Earth surface temperature 

Teller (1948) first emphasized that Dirac hypothesis 

may be in conflict with paleontological evidence. His ar- 
gument is based on the estimation of the temperature at 
the center of the Sun Tq oc GMq/Rq using the virial 
theorem. The luminosity of the Sun is then propor- 
tional to the radiation energy gradient times the mean 
free path of a photon times the surface of the Sun, that 
is Lq cx TqRqMq'^ , hence concluding that Lq cx TgAfg. 
Computing the radius of the Earth orbit in Newtonian 
mechanics, assuming the conservation of angular momen- 
tum (so that GMg i?Earth IS coustaut) and stating that 
the Earth mean temperature is proportional to the fourth 
root of the energy received, he concluded that 

TEarth OC G^'^' M^'\ (154) 

If Mq is constant and G was 10% larger 300 million years 
ago, the Earth surface temperature should have been 20% 
higher, that is close to the boiling temperature. This 
was in contradiction with the existence of trilobites in 
the Cambrian. 

Teller (1948) used a too low value for the age of the 
universe. Gamow (1967a) actualized the numbers and 
showed that even if it was safe at the Cambrian era, there 
was still a contradiction with bacteria and alga estimated 
to have lived 4 x 10^ years ago. It follows that 

\AG/G\<0.1. (155) 

Eichendorf and Reinhardt (1977) re-actualized Teller's 
argument in light of a new estimate of the age of the uni- 
verse and new paleontological discoveries to get \G/G\ < 
2.0 X 10"" yr"^ (cited by Petley, 1985). 

When using such an argument, the heat balance of 
the atmosphere is affected by many factors (water vapor 
content, carbon dioxide content, circulatory patterns,. . . ) 
is completely neglected. This renders the extrapolation 
during several billion years very unreliable. For instance, 
the rise of the temperature implies that the atmosphere 
is at some stage mostly composed of water vapor so that 
its convective mechanism is expected to change in such 
a way to increase the Earth albedo and thus to decrease 
the temperature! 

2. Expanding Earth 

Egeyed (1961) flrst remarked that paleomagnetic data 
could be used to calculate the Earth paleoradius for dif- 
ferent geological epochs. Under the hypothesis that the 
area of continental material has remained constant while 
the bulk of the Earth has expanded, the determination 
of the difference in paleolatitudes between two sites of 
known separation give a measurement of the paleoradius. 
Creer (1965) showed that data older than 3 x 10^ years 
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form a coherent group in r^.^,^ and Wesson (1973) con- 
cluded from a compilation of data that the expansion was 
most probably of 0.66 mm per year during the last 3 x 10^ 
years. 

Dicke (1962c, 1964) related the variation of the Earth 
radius to a variation of the gravitational constant by 



Alnr„ 



-O.lAlnG. 



(156) 



McElhinny et al. (1978) re-estimated the paleoradius of 
the Earth and extended the analysis to the Moon, Mars 
and Mercury. Starting from the hydrostatic equilibrium 

equation 



dP 
dr 



-G 



p{r)M{r) 



(157) 



where M(r) is the mass within radius r, they generalized 
Dicke's result to get 



A In r„ 



-aAlnG 



(158) 



where a depends on the equation of state P{p), e.g. 
a = l/(3n — 4) for a polytropic gas, P = Cp^ . In the 
case of small planets, one can work in a small gravita- 
tional self-compression limit and set P = Ko{p/pQ — 1). 
Eq. (157)then gives a = {2/15){Ap/ po), Ap being the 
density difference between the center and surface. This 
approximation is poor for the Earth and more sophis- 
ticated model exist. They give a^arth = 0.085 ± 0.02, 
aM„. = 0.032, a,_^ = 0.02 ±^0.05 and = 
0.004±0.001. Using the observational fact that the Earth 
has not expanded by more than 0.8% over the past 4 x 10* 
years, the Moon of 0.06% over the past 4 x 10^ years and 
Mars of 0.6%, they concluded that 



-G/G<8x 10-1^ yr 



(159) 



Despite any real evidence in favor of an expanding Earth, 
the rate of expansion is also limited by another geophys- 
ical aspect, i.e. the deceleration of the Earth rotation. 

Dicke (1957) listed out some other possible conse- 
quences on the scenario of the formation of the Moon 
and on the geomagnetic field but none of them enable 
to put serious constraints. The palcontological data give 
only poor limits on the variation of the gravitational con- 
stant and even though the Earth kept a memory of the 
early gravitational conditions, this memory is crude and 
geological data are not easy to interpret. 



B. Planetary and stellar orbits 

Vinti (1974) studied the dynamics of two-body system 
in Dirac cosmology. He showed that the equation of mo- 
tion 



dV 
di2 



^ = -Go- 



k + t 



to r 



(160) 



where fc is a constant and Gq the gravitational constant 
today, can be integrated. For bounded orbits, the so- 
lution describes a growing ellipse with constant eccen- 
tricity, e, pericentcr argument, u and a linearly growing 
semi-latus rectum p(t) = {l"^ /Gom){k + t)/{k + to), where 
I is the constant angular momentum, of equation 



r = 



Pit) 



1 + e cos{9 — (jj) ' 



(161) 



Similarly, Lynden-Bell (1982) showed that the equations 
of motion of the TV-body problem can be transformed to 
the standard equation if G varies as t~^. 

It follows that in the Newtonian limit, the orbital pe- 
riod of a two-body system is 



P 



2ttI 



{Gmf (1 - e2)3/2 



Ghn_ 

r.212 



2 \ 1 



(162) 



in which the correction terms represent the post- 
Newtonian corrections to the Keplerian relationship. It 
is typically of order 10~^ and 10~^ respectively for So- 
lar system planetary orbits and for a binary pulsar. It 
follows that 



P 
P 



^G 



(163) 



Only for the orbits of bodies for which the gravitational 
self-energy can be neglected does the previous equation 
reduce to 



P G 



(164) 



This leads to two observable effects in the Solar system 
(Shapiro, 1964; Counselman and Shapiro, 1968). First, 
the scale of the Solar system changes and second, if G 
evolves adiabatically as G = Go + Go{t — to), there will be 
a quadratically growing increment in the mean longitude 
of each body. 

For a compact body, the mass depends on G as well 
as other post-Newtonian parameters. At first order in 
the post-Newtonian expansion, there is a negative con- 
tribution (29) to the mass arising from the gravitational 
binding energy and one cannot neglect m in Eq. (163). 
This is also the case if other constants are varying. 



1. Early works 

Early works mainly focus on the Earth-Moon system 

and try to relate a time variation of G to a variation of 
the frequency or mean motion (n = 2tt/P) of the Moon 
around the Earth. Arguments on an expanding Earth 
also raised interests in the determination of the Earth 
rotation rate. One of the greatest problem is to evaluate 
and subtract the contribution of the spin-down of the 
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Earth arising from the friction in the seas due to tides 
raised by the Moon [Van Flandern (1981) estimated that 
'^tidai = (— 28. 8±1. 5)" century"^] and a contribution from 
the Moon recession. 

The determination of ancient rotation rates can rely 
on paleontological data, ancient ecUpse observations as 
well as measurements of star declinations (Newton, 1970, 
1974). It can be concluded from these studies that there 
were about 400 days in a year during the Devonian. In- 
deed, this studies are entailed by a lot of uncertainties, for 
instance, Runcorn (1964) compared telescope observa- 
tion from the 17th century to the ancient eclipse records 
and found a discrepancy of a factor 2. As an example, 
MuUer (1978) studied eclipses from 1374BC to 1715AD 
to conclude that 

G/G = (2.6 ± 15) X 10-^Vr"^ (165) 

and Morrison (1973) used ephemeris from 1663 to 1972 
including 40,000 Lunar occupations from 1943 to 1972 
to deduce that 

|G/G| < 2 X 10-"yr-^ (166) 

Paleontological data such as the growth rhythm found 
in fossil bivalves and corals also enable to set constraint 
on the Earth rotational history and the Moon orbit (Van 
Diggelen, 1976) [for instance, in the study by Scrutton 
(1965) the fossils showed marking so fine that the phases 
of the Moon were mirrored in the coral growth]. Blake 
(1977b) related the variation of the number of sidereal 
days in a sidereal year, Y = riE/ns, and in a sidereal 
month, M = he/um, {nE,ns and um being respectively 
the orbital frequencies of the motion of the Earth, of the 
Moon around the Earth and of the Earth around the Sun) 
to the variation of the Newton constant and the Earth 
momentum of inertia / as 

^,AY AM A7 „AG 

with 7 = 1.9856 being a calculated constant. The fossil 
data represent the number of Solar days in a tropical year 
and in a synodic month which can be related to Y and 
M so that one obtains a constraint on AI/I + 2AG/G. 
Attributing the variation of / to the expansion of the 
Earth (Wesson, 1973), one can argue that AI/I repre- 
sents only 10-20% of the r.h.s of (167). Blake (1977b) 
concluded that 

G/G= (-0.5 ±2) X lO-^Vr'^ (168) 

Van Flandern (1971, 1975) studied the motion of the 
Moon from Lunar occultation observations from 1955 to 
1974 using atomic time which differs from the ephemeris 
time relying on the motion of the Earth around the Sun. 
He attributed the residual acceleration after correction of 



tidal effect to a variation of G, '^Moon/^'^Moon — ^ 
5) X 10^^ ccntury^^ to claim that 

G/G = (-8 ±5) X 10"" yr-^ (169) 

In a new analysis, Van Flandern (1981) concluded that 
»^Moon/»^Moon = (3.2±1.1) X 10"" yr"^ hence that G was 
increasing as 

G/G= (3.2 ±1.1) X 10-"yr-i (170) 

which has the opposite sign. In this comparison the time 
scale of the atomic time is 20 years and the one of the 
ephemeris 200 years but is less precise. It follows that 
the comparison is not obvious and that these results are 
far from being convincing. In this occultation method, 
one has to be sure that the proper motions of the stars 
are taken into account. One also has to assume that (1) 
the mass of the planets are not varying (see Eq. 163), 
which can happen if e.g. the strong and fine structure 
constants are varying (2) the fine structure constant is 
not varying while comparing with atomic time and (3) 
the effect of the changing radius of the Earth was not 
taken into account. 

2. Solar system 

Monitoring the separation of orbiting bodies offers a 
possibility to constrain the time variation of G. This 
accounts for comparing a gravitational time scale (set by 
the orbit) and an atomic time scale and it is thus assumed 
that the variation of atomic constants is negligible on the 
time of the experiment. 

Shapiro et al. (1971) compared radar-echo time de- 
lays between Earth, Venus and Mercury with a cesium 
atomic clock between 1964 and 1969. The data were fit- 
ted to the theoretical equation of motion for the bodies 
in a Schwarzschild spacetime, taking into account the 
perturbations from the Moon and other planets. They 
concluded that 

|G/G| < 4 X 10-^°yr"^ (171) 

The data concerning Venus cannot be used due to impre- 
cision in the determination of the portion of the planet 
reflecting the radar. This was improved to 

|G/G| < 1.5 X 10-1° yr"^ (172) 

by including Mariner 9 and Mars orbiter data (Reascn- 
berg and Shapiro, 1976, 1978). The analysis was further 
extended (Shapiro, 1990) to give 

G/G = (-2 ± 10) X 10-1^1-"^ (173) 

The combination of Mariner 10 an Mercury and Venus 
ranging data gives (Anderson et al., 1991) 
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G/G= (0.0±2.0) X 10"^^ yr 



(174) 



3. Pulsars 



The Lunar laser ranging (LLR) experiment has mea- 
sured the position of the Moon with an accuracy of about 
1 cm for thirty years. This was made possible by the 
American Appolo 11, 14 and 15 missions and Soviet- 
French Lunakhod 1 and 4 which landed retro-reflectors 
on the Moon that reflect laser pulse from the Earth 
(see Dickey et al. (1994) for a complete description). 
Williams et al. (1976) deduced from the six first years of 



LLR data that > 29 so that 



\G/G\ < 3 X 10-'' yv 



(175) 



Miillcr et al. (1991) used 20 years of data to improve this 

result to 



\G/G\ < 1.04 X 10-"yr-\ 



(176) 



the main error arising from the Lunar tidal acceleration. 
Dickey et al. (1994) improved this constraint to 



\G/G\ < 6 X 10-12 yr"^ 



(177) 



and Williams et al. (1996) with 24 years of data con- 
cluded that 



\G/G\ < 8 X 10-12 yr-i 



(178) 



Rcascnbcrg et al. (1979) considered the 14 months 
data obtained from the ranging of the Viking spacecraft 
and deduced that co^^ > 500 which implies 



-G/G < 10-12 yr-i. 



(179) 



Hellings et al. (1983) using all available astrometric data 
and in particular the ranging data from Viking landers 
on Mars deduced that 



\G/G\ = (2±4) X 10-12 yr"^ 



(180) 



The major contribution to the uncertainty is due to the 
modeling of the dynamics of the asteroids on the Earth- 
Mars range. Hellings et al. (1983) also tried to attribute 
their result to a time variation of the atomic constants. 
Using the same data but a different modeling of the as- 
teroids, Reasenberg (1983) got 

< 3 X 10-iiyr-i (181) 

which was then improved by Chandler et al. (1993) to 

|G/G| < 10-11 yr-i. (182) 

All these measurements allow to test more than just 
the time variation of the gravitational constant and offer 
a series of tests on the theory of gravitation and constrain 
PPN parameters, geodetic precession etc... (see Will, 
1993). 



Contrary to the Solar system case, the dependence 

of the gravitational binding energy cannot be neglected 
while computing the time variation of the period (Dicke, 
1969; Eardley, 1975; Haugan, 1979). Here two ap- 
proaches can be followed; either one sticks to a model 
(e.g. scalar-tensor gravity) and compute all the effects in 
this model or one has a more phenomenological approach 
and tries to put some model- independent boimds. 

Eardley (1975) followed the first route and discussed 
the effects of a time variation of the gravitational con- 
stant on binary pulsar in the framework of the Brans- 
Dicke theory. In that case, both a dipole gravitational 
radiation and the variation of G induce a periodic varia- 
tion in the pulse period. Nordtvedt (1990) showed that 
the orbital period changes as 



P 
P 



2 + 



2(mici -I- m2C2) -|- 3(miC2 + m2Ci) 

TOi + 1712 



G 
G 

(183) 



where Cj = i51nmi/(51nG. He concluded that for the 
pulsar PSR 1913-1-16 (mi ~ m2 and ci ~ C2) one gets 



| = -P + 5cl§, 



(184) 



the coefficient c being model dependent. As another 
application, he estimated that Cj,^_.^^ ~ — 5 x lO^i", 
Cmoo„ ~ —10-^ and Cg^„ ~ — 4 x 10-^ justifying the ap- 
proximation (164) for the Solar system. 

Damour et al. (1988) used the timing data of the bi- 
nary pulsar PSR 1913-1-16. They implemented the effect 
of the time variation of G by considering the effect on 
P/P and making use of the transformation suggested by 
Lynden-Bell (1982) to integrate the orbit. They showed, 
in a theory-independent way, that G/G = —0.5SP/P, 
where SP is the part of the orbital period derivative that 
is not explained otherwise (by gravitational waves radi- 
ation damping). It has to be contrasted with the result 
(184) by Nordtvedt (1990). They got 



G/G= (1.0 ±2.3) X 10-11 yr 



(185) 



Damour and Taylor (1991) reexamined the data of 
PSR 1913-1-16 and the upper bound 



G/G < (1.10 ± 1.07) X 10-11 yr-i. 



(186) 



Kaspi et al. (1994) used data from PSR B1913-I-16 and 
PSR B 1855+09 respectively to get 



G/G = {A ±5) X 10-12 yr- 



and 
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G/G = (-9 ±18) X 10-12 yr 



(187) 



(188) 
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the latter case being more "secure" since the orbiting 
companion is not a neutron star. 

All the previous results concern binary pulsar but iso- 
lated can also be used. Heintzmann and Hillcbrandt 
(1975) related the spin-down of the pulsar JP1953 to a 
time variation of G. The spin-down is a combined ef- 
fect of electromagnetic losses, emission of gravitational 
waves, possible spin-up due to matter accretion. As- 
suming that the angular momentum is conserved so that 
I/P = const ant, one deduces that 



dln/\ G 



Pg VdlnCy G 
The observational spin-down can be decomposed as 



(189) 



P 



P 



P P 

GW J Cr 



(190) 



Since P/P^^^ and P/P^^ are positive definite, it follows 
that P / P,,^.^ > P / Pg ho that a bound on G can be in- 
ferred if the main pulse period is the period of rotation. 

Heintzmann and Hillebrandt (1975) modeled the pul- 
sar by a polytropic (P oc p") white dwarf and deduced 
that dln//dlnG = 2 - 3n/2 so that 



\G/G\ < 10"^" yr 



(191) 



Mansfield (1976) assumed a relativistic degenerate, zero 
temperature polytropic star and got 



-G/G < 5.8t J X 10-11 



(192) 



at a 2(7 level. He also noted that a positive G induces 
a spin-up counteracting the electromagnetic spin-down 
which can provide another bound if an independent esti- 
mate of the pulsar magnetic field can be obtained. Gold- 
man (1990), following Eardley (1975), used the scaling 
relations N oc G'^/"^ and M oc G'^/"^ to deduce that 
2dln7/dlnG = + 3dln7/dln A/'. He used the data 
from the pulsar PSR 0655+64 to deduce 



-G/G < (2.2 - 5.5) X 10"" yr 



(193) 



C. Stellar constraints 

In early works, Pochoda and Schwarzschild (1964), 
Ezer and Cameron (1966) and then Gamow (1967c) stud- 
ied the Solar evolution in presence of a time varying grav- 
itational constant. They came to the conclusion that un- 
der Dirac hypothesis, the original nuclear resources of the 
Sun would have been burned by now. This results from 
the fact that an increase of the gravitational constant is 
equivalent to an increase of the star density (because of 
the Poisson equation). 



A side effect of the change of luminosity is a change in 
the depth of the convection zone. This induces a mod- 
ification of the vibration modes of the star and partic- 
ularly to the acoustic waves, i.e p-modes, (Demarque et 
al., 1994). Demarque et al. (1994) considered an ansatz 
in which G oc and showed that |/3| < 0.1 over the last 
4.5 X 10^ years, which corresponds to 



G/G 



< 2 X 10-11 yr" 



(194) 



Guenther et al. (1995) also showed that ^r-modes could 
provide even much tighter constraints but these modes 
are up to now very difficult to observe. Nevertheless, 
they concluded, using the claim of detection by Hill and 
Gu (1990), that 



G/G 



< 4.5 X 10-1^ yr 



(195) 



Guenther et al. (1998) compared the p-mode spectra 

predicted by different theories with varying gravitational 
constant to the observed spectrum obtained by a network 
of six telescopes and deduced that 



G/G < 1.6 X 10- 



■^Vr-i. 



(196) 



The standard Solar model depends on few parameters 
and G plays a important role since stellar evolution is 
dictated by the balance between gravitation and other 
interactions. Astronomical observations determines very 
accurately GMq and a variation of G with GMq fixed 
induces a change of the pressue (P = GMq/Rq) and 
density {p = Mq/R^). Ricci and Villante (2002) studied 
the effect of a variation of G on the density and pressure 
profile of the Sun and concluded that present data cannot 
constrain G better than 10—2%. 

The late stages of stellar evolution are governed by the 
Chandrasekhar mass {fi.c/G)'^^^m^^ mainly determined 
by the balance between the Fermi pressure of a degen- 
erate electron gas and gravity. Assuming that the mean 
neutron star mass is given by the Chandrasekhar mass, 
one expects that G/G = -2M„s/3M^s. Thorsett (1996) 
used the observations of five neutron star binaries for 
which five Keplerian parameters can be determined (the 
binary period Pf,, the projection of the orbital semi-major 
axis oi sini, the eccentricity e, the time and longitude of 
the periastron Tq and uj) as well as the relativistic ad- 
vance of the angle of the periastron uj. Assuming that the 
neutron star masses vary slowly as M„g = M^J — M^^ t^^ , 
that their age was determined by the rate at which Pf, 
is increasing (so that tNS — '^Pb/Pb) and that the mass 
follows a normal distribution, Thorsett (1996) deduced 
that, at 2a, 



G/G = (-0.6 ± 4.2) X 10-12 yr" 



(197) 



Analogously, the Chandrasekhar mass sets the charac- 
teristic of the light curves of supernovae (Riazuelo and 
Uzan, 2002). 
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Garcia-Bcrro et al. (1995) considered the effect of a 
variation of the gravitational constant on the cooUng of 
white dwarfs and on their luminosity function. As first 
pointed out by Vila (1976), the energy of white dwarfs 
is entirely of gravitational and thermal origin so that a 
variation of G will induce a modification of their energy 
balance. Restricting to cold white dwarfs with luminosity 
smaller than ten Solar luminosity, the luminosity can be 
related to the star binding energy B and gravitational 
energy, as 



(198) 



which simply results from the hydrostatic equilibrium. 
Again, the variation of the gravitational constant inter- 
venes via the Poisson equation and the gravitational po- 
tential. The cooling process is accelerated if G/G < 
which then induces a shift in the position of the cut-off 
in the luminosity function. Garcia-Berro et al. (1995) 
concluded that 



-G/G < 3tl 



X 10"" yr" 



(199) 



The result depends on the details of the cooling theory 
and on whether the C/0 white dwarf is stratified or not. 

A time variation of G also modifies the main sequence 
time of globular clusters (Dicke 1962a; Roeder, 1967). 
Del'Innocenti et al. (1996) calculated the evolution of 
low mass stars and deduced the age of the isochrones. 
The principal effect is a modification of the main se- 
quence evolutionary time scale while the appearance 
of the color-magnitude diagram remained undistorted 
within the observational resolution and theoretical un- 
certainties. Since the globular clusters must be younger 
than the universe, and assuming that their age was be- 
tween 8 and 20 Gyr, they concluded 



G/G = (-1.4 ± 2.1) X 10"^^ yr" 



(200) 



This analysis was also applied to clusters of galaxies by 
Dearborn and Schramm (1974). In that case a lower 
gravitational constant allows the partic;le to escape from 
the cluster since the gravitational binding energy also 
decreases. They deduced that the decrease of G that 
allows the existence of clusters at the present epoch is 



-G/G<4x 10-^Vr"^ 



(201) 



(1) The variation of G modifies the Friedmann equa- 
tion and therefore the age of the Universe (and, hence, 
the sound horizon). For instance, if G is larger at earlier 
time, the age of the Universe is smaller at recombina- 
tion, so that the peak structure is shifted towards higher 
angular scales. 

(2) The amplitude of the Silk damping is modified. At 
small scales, viscosity and heat conduction in the photon- 
baryon fluid produce a damping of the photon perturba- 
tions (Silk, 1968). The damping scale is determined by 
the photon diffusion length at recombination, and there- 
fore depends on the size of the horizon at this epoch, and 
hence, depends on any variation of the Newton constant 
throughout the history of the Universe. 

(3) The thickness of the last scattering surface is mod- 
ified. In the same vein, the duration of recombination is 
modified by a variation of the Newton constant as the 
expansion rate is different. It is well known that CMB 
anisotropics are affected on small scales because the last 
scattering "surface" has a finite thickness. The net ef- 
fect is to introduce an extra, roughly exponential, damp- 
ing term, with the cutoff length being determined by the 
thickness of the last scattering surface. When translating 
redshift into time (or length), one has to use the Fried- 
mann equations, which are affected by a variation of the 
Newton constant. The relevant quantity to consider is 
the visibility function g. In the limit of an infinitely thin 
last scattering surface, t goes from oo to at recombi- 
nation epoch. For standard cosmology, it drops from a 
large value to a much smaller one, and hence, the visibil- 
ity function still exhibits a peak, but is much broader. 

Chen and Kamionkowski (1999) studied the CMB 
spectrum in Brans-Dicke theory and showed that CMB 
experiments such as MAP will be able to constrain these 
theories for uj-^-^ < 100 if all parameters are to be de- 
termined by the same CMB experiment, < 500 if 
all parameters are fixed but the CMB normalization and 
Lu^j^ < 800 if one uses the polarization. For the Planck 
mission these numbers are respectively, 800, 2500 and 
3200. 

As far as we are aware, no complete study of the impact 
of the variation of the gravitational constant (e.g. in 
scalar-tensor theory) on the CMB has been performed 
yet. Note that, to compute the CMB anisotropics, one 
needs not only the value of G at the time of decoupling 
but also its complete time evolution up to now, since it 
will affect the integrated Sachs- Wolfe effect. 



D. Cosmological constraints 

1. CMB 

A time-dependent gravitational constant will have 
mainly three effects on the CMB angular power spectrum 
(Riazuelo and Uzan, 2002): 



2. Nucleosynthesis 

As explained in details in section III.C.2, changing the 
value of the gravitational constant affects the freeze-out 

temperature Tf. A larger value of G corresponds to a 
higher expansion rate. This rate is determined by the 
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combination Gp and in the standard case the Friedmann 
equations imply that Gpt^ is constant. The density p 
is determined by the number TV* of relativistic particles 
at the time of nucleosynthesis so that nucleosynthesis al- 
lows to put a bound on the number of neutrinos iV^. 
Equivalently, assuming the number of neutrinos to be 
three, leads to the conclusion that G has not varied from 
more than 20% since nucleosynthesis. But, allowing for 
a change both in G and allows for a wider range of 
variation. Contrary to the fine structure constant the 
role of G is less involved. 

Steigmann (1976) used nucleosynthesis to put con- 
straints on the Dirac theory. Barrow (1978) assumed 
that G oc and obtained from the helium abundances 
that -5.9 X 10~^ < n< 7 X 10~^ which implies that 



G/G < (2±9.3)/ix IQ-^^yv- 



(202) 



assuming a flat universe. This corresponds in terms of 
the Brans-Dicke parameter to > 25, which is a much 
smaller bounds that the ones obtained today. Yang et 
al. (1979) included the computation of the deuterium 
and lithium. They improved the result by Barrow (1978) 
to n < 5 X 10~^ which corresponds to lo^^ > 50 and 
also pointed out that the constraint is tighter if there are 
extra-neutrinos. It was further improved by Rothman 
and Matzner (1982) to |n| <3 x 10~^ implying 



G/G 



< 1.7 X 10"" yr 



(203) 



Accetta et al. (1990) studied the dependence of the abim- 
dances of D, "^He, ''He and ^Li upon the variation of G 
and concluded that 



-0.3 < AG/G < 0.4 



(204) 



which roughly corresponds to9xl0 "^<n<8xl0 
and to \G/G\ < 9 x lO-^^yr"!. 

All previous investigations assumed that the other con- 
stants are kept fixed and that physics is unchanged. Kolb 
et al. (1986) assumed a correlated variation of G, a^^ 
and Gp and got a bound on the variation of the radius 
of the extra-dimensions. 

The case of Brans-Dicke (1961) theory, in which only 
the gravitational constant varies, was well studied. Casas 
et al. (1992a, 1992b) concluded from the study of hehum 
and deuterium abundances that cjj^^ > 380 when N^, = 3 
(see also Damour and Gundlach, 1991, and Serna et al., 
1992) and w^d > 50 when = 2. 

Kim and Lee (1995) calculated the allowed value for 
the gravitational constant, electron chemical potential 
and entropy consistent with observations up to lithium- 
7 and argued that beryllium-9 and bore- 11 abundances 
are very sensitive to a change in G. Kim et al. (1998) 
further included neutrino degeneracy. The degeneracy 
of the electron-neutrino not only increases the radiation 



density but also influences the weak interaction rates so 
that it cannot be absorbed in a variation of G. It was 
shown that a higher gravitational constant can be bal- 
anced by a higher electron-neutrino degeneracy so that 
the range of (electron chemical potential, G) was wider. 

Damour and Pichon (1999) extended these investiga- 
tions by considering a two-parameter family of scalar- 
tensor theories of gravitation involving a non-linear scalar 
field-matter coupling function. They concluded that even 
in the cases where before BBN the scalar-tensor theory 
was far from general relativity, BBN enables to set quite 
small constraints on the observable deviations from gen- 
eral relativity today. 

Let us also note the work by Carroll and Kaplinghat 
(2001) in which they tried to constrain the expansion his- 
tory of our universe in a model-independent way during 
nucleosynthesis. They assumed changes in the gravita- 
tional dynamics and not in the particle physics processes. 
For that purpose the expansion rate at the time of nucle- 
osynthesis is approximated as H{T) = (T/1 MeV)"Hi in 
order to infer the constraints on {a, Hi). This a simple 
way to compare an alternative to cosmology with data. 



V. OTHER CONSTANTS 

Up to now, we have detailed the results concerning the 
two most studied constants, and a^.^^. But, as we em- 
phasized, if a^.^, is varying one also expects a variation of 
other constants such as and a^^, . There are many the- 
oretical reasons for that. First, in Kaluza-Klein or string 
inspired models, all constants are varying due either to 
the dilaton or the extra-dimensions (see Section VI for 
details). 

Another argument lies in the fact that if we believe in 

grand unifled theories, there exists an energy scale A^y.^, 
at which all the (non-gravitational) couplings unify. 



[(■'^gut) '-'^w (■'^gut ) ^s(-^gut) — '-'^Q 



(205) 



The value of the coupling constants at any energy scale 
smaller than A^^^j, is obtained from the renormalization 
group equations. It follows that a time variation of a^^ 



induces a time variation of a^^T thus of 



and ttg . 



In such a framework, the varying parameters would then 
be ttcuTi ^gut/-^4 ^iid Yukawa couplings. 

The strong coupling at an energy scale E is related to 
the QCD scale A.^^^ by 



27r 



(206) 



with /3o = — 11 + 2nf/3, rif being the number of quark 
flavors. It follows that 



AA, 



In 



E \ 



A J a 

QCD / "S 



(207) 
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The time variation of is thus not the same at ah en- 
ergy scales. In the chiral Umit, in which the quarks are 
massless, the proton mass is proportional to the QCD 

energy scale, rup oc A^^^,, so that a change in (or in 
o:^-^J,J,) induces a change in /i and we have 



Amp/r 



•■QCD- 



(208) 



The energy-scale evolution of the three coupling con- 
stants in a 1-loop approximation takes the form 



2lT 



In 



E 



(209) 



where the numerical coefficients depend on the choice 
of the considered gauge group. For instance 6j = 
(41/10,-19/16,-7) in the standard model (SM) and 
hi = (33/5,1,-3) in its minimal supersymmetric exten- 
sion. In the case of supersymmetric models (SUSY), 
Eq. (209) has to be replaced by 



27r 



In 



E 
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Using (209), one can work out the variation of all cou- 
plings once the grand unified group is chosen, assuming 
or not supersymmctry. 

In the string-inspired model by Damour and Polyakov 
(1994a,b), Eq. (146) [obtained from Eq. (276)] implies 
that 



Amhadron/mhadron ^ 'iQAa^^/u^ 



(210) 



as was first pointed out by Taylor and Vencziano (1988). 

Recently Calmet and Fritzsch (2001,2002), Dent and 
Fairbairn (2001) and Langacker et al. (2002) tried to 
work out these relationships in different models and 
confirmed the order of magnitude (210). Calmet and 
Fritzsch (2001) computed low energy effects of a time 
varying fine structure constant within a GUT-likc theory 
with a constraint of the form (205) and focus their anal- 
ysis on the nucleon mass. Nevertheless, they assumed 
that the mechanisms of clcctroweak and supersymmctry 
breaking as well as fermion mass generation were left un- 
changed and thus that quarks, leptons, W and Z masses 
do not vary. But, as seen from Eqs (211-212) below, 
one cannot vary with being fixed without vary- 
ing the Higgs vacuum expectation value which induces 
a variation of the mass of the fermions. On this basis 
they concluded that the result by Webb et al. (2001) 
on the fine structure constant implies that Amp/nip ~ 
38 Aa^j^j /oj.^ ~ —4 x 10~^ (keeping the Planck mass 
constant) and that Ay/y ~ — 121Aag^/ag^ ~ 3 x 10"'', 
which is above the current observational constraints (see 



Section V.D). Calmet and Fritzsch (2002) considered 
different scenarios: (i) A^jj^, is constant and a^jj^ time- 
dependent, (ii) only A^^^t is time-dependent and (iii) 
both are varying. They concluded that the most "in- 
teresting" situation, in view of the variation of a^,^ and 
H, is the second case. Langacker et al. (2002) pointed 
out that changes in the quark masses and in the Higgs 
vacuum expectation value were also expected and they 
parameterized the effects of the variation of qigut on 
the electroweak and Yukawa sector. They assumed that 
<^GUT the vacuum expectation value of a slowly vary- 
ing scalar field. They concluded that AA^^^/A^^^ 
34 Aq!jj,„ /q:j,^j (with a precision of about 20% on the nu- 
merical factor) and that a variation of the fine structure 
of the magnitude of the one observed by Webb et al. 
(2001) would imply Amp/mp ~ -2.5 x 10"''. They also 
argued that Ax/x ~ —32Aa^j^/a^^ ~ 8x 10~^, which is 
consistent with current bounds, if one assumes the vari- 
ation of the proton gyromagnetic factor to be negligible. 
Earlier, Sisterna and Vucetich (1990) tried to determine 
the compatibility of all the bounds by restricting their 
study to ( Aq^j3 , Q!jj,^j , Gp , G) and then included the u, d 
and s quark masses (Sisterna and Vucetich, 1991). 

Since we do not have the theories of electroweak and 
supersymmctry breaking as well as the ones for the gen- 
eration of fermion masses, the correlations between dif- 
ferent low-energy observables remain model-dependent. 
But, in this unification picture, one is abale to derive 
stronger constraints. For instance Olive et al. (2002) 
expressed the constraints from a-, /3-decays and Oklo in 



fonction of \AA^^^/A 



■■QCD 



Av/v\ ~ 50AaEM/aEM to 
give respectively the tighter constraints < 10"'', < 10~^ 
and < 10"^'^. The goal of this section is to discuss the 
constraints on some of these constantss which are of im- 
portance while checking for consistency. 



A. Weak interaction 

Most of the studies on the variation either of Gp or 
concern BEN, Oklo, CMB and geochemichal dating. 

The Fermi constant can be expressed in terms of 
and of the mass of the boson W, , as 



9 



8M2 ■ 



(211) 



In the standard model, is simply the product of 5^/4 
by the Higgs vacuum expectation value = {(j))'^ , so that 

Gp = l/2v^. (212) 

Thus, at tree level, Gp is actually independent of the 
SU (2) coupling and is a direct measurement of the mag- 
nitude of the electroweak symmetry breaking. Note that 
a change in v is related to a change in the Yukawa cou- 
plings. 
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Kolb et al. (1985) considered the effect of the vari- 
ations of the fundamental constants on nucleosynthesis. 
As detailed in Section III.C.2, they found the dependence 
of the helium abundance on G, Gp and Q, the variation 
of which were related to the variation of a^^ (then re- 
lated to the size of extra-dimensions). Kolb et al. (1985) 
did not considered changes in Gp due to the variation in 
and assumed that 5G^ oc 5g^. Since G, 0;^^, and 
Gp where related to the volume of the extra-space, this 
study gives no bound on the variation of Gp . 

Dixit and Sher (1988) pointed out that the relation be- 
tween Gp and in the work by Kolb et al. (1985) was 
ill-motivated and that the only way to vary Gp was to 
vary v. Changing v has four effects on BBN: it changes 
(1) all weak interaction rates, (2) me, (3) the quark 
masses and hence Q and (4) the pion mass which affects 
the strong nuclear force and the binding of the deuteron. 
Using results on the dependence of rUe and on 
and V (Gasser and Leutwyler, 1982) they got 

— ^ = 1.293 -0.9^^ -h 2.193— . (213) 
IMeV V ' 

Besides, a change of 1% of the quark masses changes the 
pion mass by 0.5% which implies that the deuteron bind- 
ing energy changes also by 0.5% (Davies, 1972). They 
concluded that the helium abundance was given by 



Fp = 0.240 - 0.31At;/t; -|- 0.38AQ!p,^,/Q!p 



(214) 



and deduced that Au/t; < 0.032 if is fixed and 
Aofp^/apj, < 0.026 if v is fixed. They also noted that 
the changes in Q) and Gp induced by v tend to cancel 
making the change in a^^, appearing only in Q, larger. 

Scherrer and Spergel (1993) followed the same way and 
focused on two cases: (1) that the Yukawa couplings are 
fixed so that both Gp and the fermion masses vary in 
parallel and (2) that the Yukawa couplings vary so that 
Gp changes while the fermion masses are kept constant. 
Considering the abundances of helium they deduced, as- 
suming Q!g^^ fixed, that 



-0.22 < AGp/Gp < 0.01 
in the first case and 

-0.01 < AGp/Gp < 0.09 



(215) 



(216) 



in the second case. 

To finish with cosmological constraints, a change in Gp 
induces a change in me which can be constrained by the 
CMB. The electron mass appears in the expression of the 
Thomson cross section (116) and on the binding energy of 
hydrogen (117) which induces a change in the ionization 
fraction. Kujat and Scherrer (2000) implemented these 
changes as in Section III.C and showed that the upper 
limit on Ame/me is of order 10~^ — 10~^ (keeping the 



Planck mass constant) for a maximum multipole of € ^ 
500 — 2500 if a^^ is assumed constant. The degeneracy 
with a^yi is broken at high multipole so that one can hope 
to detect a 1% variation with a maximum multipole of 
£ > 1500. 

Prom Oklo data, Shlyakhter (1976) argued that the 

weak interaction contribution to the total energy of the 
nucleus is of order 10~^(m7r/mp)^ so that Ag^/g^ ~ 
5 X W^Ag^/g^ to conclude that 



|Aa„/a„| < 4 X 10" 



(217) 



But in fact, the change in is much more difficult 
to model than the change in a^^. Damour and Dyson 
(1996) used the estimate by Haugan and Will (1976) for 
the weak interaction contribution to the nuclear ground 
state energy of samarium £^(^^°Sm) — £(^^^Sm) ~ 5.6 eV 
to conclude that, if no subtle cancellation appears. 



\Aa^/a^\ < 0.02. 



(218) 



Concerning geolochemical data (see Section V.A), 
Dyson (1972) pointed out that all /3-decay rates are pro- 
portional to so that all constraints are in fact de- 
pendent on the combination a'^^a'^, s being defined 
in Eq. (49). The degeneracy can be lifted by compar- 
ing different nuclei, e.g. 75 ''Re (sRe = —18000) and 
IgK (sk = --30). The constancy of the decay rates of 
these two nuclei have approximatively the same accuracy. 
From the constancy of the ratio 

. ARe , , ARe Aap^ 

A-— = (SRe - Sk) 

Ak Ak OpM 

within a few parts in 10^° per year, one can deduce that, 
independently of any assumption on a^, 



< 2 X 10"^ 



and thus that 



\Aa^/a^\ < 10 



(219) 



(220) 



during the last 10^ years. 

Wilkinson (1958) studied the variation of by using 
pleochroic haloes, that is spheres formed by a-ray tracks 
around specks of uranium-bearing mineral in mica. The 
intensities of the haloes of different radii give a picture of 
the natural radioactive series integrated over geological 
time from which one can deduce the proportion of differ- 
ent daughter-activities in the decay chain from uranium 
to lead. This series contain elements undergoing both a- 
and /3-decay. For instance Ac branches 1.2% by a-decay 
and the rest by /3-decay. From 10^ years old samples, 
Wilkinson (1958) deduced that 



|Aa^/a^| < 10. 



(221) 
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Let us also point out some works (Agrawal et at, 
1998a, b) in which the mass scale of the standard model 
and the scale of electroweak symmetry breaking are con- 
strained by mean of the anthropic principle. Passarino 
(2001) investigated the effects of the time variation of 
the Higgs vacuum expectation value and showed that 
the classical equation of motion for the Higgs field in 
the standard model accepts time dependent solution. 

B. Strong interaction 

There is a very small number of works addressing this 
issue. Due to the strong energy dependence of a^, it 
makes more sense to constraint the variation of A^^.^,. 
It has a lot of implications on the stability properties of 
nuclei and it follows that most of the constraints arise 
from nuclear considerations. Let us remind that in the 
chiral limit, all dimensional parameters are proportional 
to A^^^ so that all dimensionless ratios will be, in this 
limit, pure numbers and thus insensitive to a change of 
the strong interaction. But, quark masses will play an 
important role in the variation of dimensionless ratios 
and have to be taken into account. 

A change in the strong interaction affects the light el- 
ements and (1) the most weakly bound nucleus, namely 
the dcutcron, can be unbind if it is weaker(2) there may 
exist stable dineutron and diproton if it is stronger [and 
hydrogen would have been burned catastrophically at the 
beginning of the imiversc, (Dyson, 1971)], (3) the rate of 
the proton capture {p + n ^ D + ^) is altered and (4) the 
neutron lifetime changes. All these effects influence the 
nucleosynthesis (Barrow, 1987). It will also be a catas- 
trophe if the deuteron was not stable (by affecting the 
the hydrogen burning properties in stars). 

Most of the early studies considered these stability 
properties by modelling the nuclear force by a Yukawa 
approximation of the form V{r) ^ exp(— m,r'')- In the 
following of this section, the cited bounds refer to suh 
a definition of a^. Davies (1972) studied the stability 
of two-nucleon systems in terms of a^^j and assum- 
ing that remains fixed and concluded that the dipro- 
ton is not bounded if Aag/ttg — Aa^^/a^^ < 0.034. 
Rozental (1980) assumed that the depth of the poten- 
tial well in the deuteron is proportional to , to state 
that a decrease of of 10-15% would make it unsta- 
ble. An increase of would render the diproton stable 
so that lAttg/agl < lO"'^ at nucleosynthesis. A previ- 
ous and more detailed analysis by Davies (1972) yield 
|Aas/ag| < 4 X 10~^ and Pochet et al. (1991) concluded 
that jAag/agl < 4 x 10~^ for the deuteron to be stable 
and I Aag/ag | < 6 x 10~^ for the diproton to be unstable. 

Concerning high-Z nuclei, Broulik an Trefil (1971) 
used the liquid drop model of the nucleus and the ob- 
served half lives and abundances of transuramium ele- 
ments to constraint the variation of a^/a^^^. In this 



model, the stability of a nucleus can be discussed by 
comparing the Coulomb repulsion between protons to 
the strong interaction attraction modeled by a surface 
tension T proportional to . With increasing atomic 
weight, the individual nucleons become progressively 
more weakly bound as the Coulomb force dominates. A 
nucleus is stable against spontaneous fission if 

If ag/o^BM larger in the past some unstable nuclei 
would have been stable. The idea is thus to find unsta- 
ble nuclei with long half-life which do not occur naturally. 
The variation of a^/a^^ would make them stable in the 
past but this must have occured roughly more than about 
ten times their lifetime since otherwise they will be in de- 
tectable abundances. Assuming that a^^ is fixed, they 
concluded from data on 94^Pu that I Aag /ttg I < 1.7x10"^ 
on a time scale of about 7.6 x 10* yr. Unfortunately, four 
month later it was reported that 94^ Pu occurs naturally 
on Earth (Hoffmann et al., 1971) hence making the ar- 
gument invalid. Davies (1972) argued that the binding 
energy is expected to vary as (contrary to the ansatz 
by Broulik and Trefil, 1971) so that the previous bound 
becomes lAag/ofgl < 8.5 x 10~^. 

Barrow (1987) studied the effect of the change of 
on the BBN predictions in Kaluza-Klein and superstring 
theories in which all the couplings depends on the com- 
pactification radius. Assuming that the deuteron binding 
energy, probably the most sensitive parameter of BBN, 
scales as , he concluded that 

iN/Tn OC G-1/2q2g2 (223) 

which affects the helium abundances from Eqs. (126-127). 
Flambaum and Shuryak (2001) discussed the effects of 
the variation of and took the quarks masses into ac- 
count. They expressed their results in terms of the pa- 
rameter = Sln{m^/A^^,^) = Sln{^{niu + md)/A^^^) 
where the pion mass tt?,^ determines the range of the nu- 
clear force. They concluded that < 0.005 between 
BBN and today. 

As detailed in Section III.A.l, Shlyakhter (1976) ar- 
gued that the change in the energy of the resonance is 
related to a change in by 

Agjg, ~ AEr/Vo (224) 

and deduced that \Ag^/g^\ < 1.9 x 10"^. Clearly, this 
analysis is not very reliable. Flambaum and Shuryak 
(2001) estimated the variation of the resonance energy 
due to a variation of the pion mass and concluded that 
AEr/Er ~ 3 X 10^\S„\ so that \d^\ < 7 x IQ-^o. 

Flambaum and Shuryak (2001) also argued that in the 
worst case all strong interaction phenomena depend on 
Aqpq + ^''T^s where K is some universal constant and 
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the strange quark mass but a real study of the effect of 
jTig on all hadronic masses remains to be done. It also 
follows that the proton gyromagnetic factor can be time 
dependent and constrained by observations such as those 
presented in Section V.D. 



C. Electron to proton mass ratio 

An early limit on the variation of fi was derived by 
Yahil (1975) who compared the concordance of K-Ar and 
Rb-Sr gcochcmical ages and deduced that |A/i//i| < 1.2 
over the past 10^° yr. 

As first pointed out by Thomson (1975) molecular ab- 
sorption lines can provide a test of the variation of /j,. The 
energy difference between two adjacent rotational levels 
in a diatomic molecule is proportional to Mr~^, r being 
the bond length and M the reduced mass, and that the 
vibrational transition of the same molecule has, in first 
approximation, a \/M dependence. For molecular hydro- 
gen M — m-p/2 so that comparison of an observed vibro- 
rotational spectrum with its present analog will thus give 
information on the variation of mp and rrin- Compar- 
ing pure rotational transitions with electronic transitions 
gives a measurement of /z. 

Following Thompson (1975), the frequency of vibra- 
tion-rotation transitions is, in the Born-Oppenheimer ap- 
proximation, of the form 



|A/i//i| < 2 X 10" 



(228) 



(225) 



where c , . c and c . arc some numerical coefficients. 

clcc ' vib rot 

Comparing the ratio of wavelengths of various electronic- 
vibration-rotational lines in quasar spectrum and in the 
laboratory allow to trace the variation of n since, at low- 
est order, Eq. (225) implies 



A£,,(0) 



(226) 



Varshalovich and Levshakov (1993) used the observations 
of a damped Lyman-a system associated with the quasar 
PKS 0528-250 (which is believed to have molecular hy- 
drogen in its spectrum) of redshift 2; = 2.811 and deduced 
that 



|A/x//^| < 4 X 10" 



(227) 



A similar analysis was first tried by Foltz et al. (1988) but 
their analysis did not take into account the wavclcngth- 
to-mass sensitivity and their result hence seems not very 
reliable. Nevertheless, they concluded that 



^^In the literature refers either to m^/m-p or to its inverse. 
In the present work we choose the first definition and we har- 
monize the results of the different articles. 



at 2; = 2.811. Cowie and Songaila (1995) observed the 
same quasar and deduced that 



A/i//i= (0.75 ±6.25) X 10" 



(229) 



at 95% C.L. from the data on 19 absorption lines. Var- 
shalovich and Potekhin (1995) calculated the coefficient 
Kij with a better precision and deduced that 



|A/x//i| < 2 X 10" 



(230) 



at 95% C.L.. Lanzetta et al. (1995) and Varshalovich et 
al. (1996b) used 59 transitions for H2 rotational levels in 
PKS 0528-250 and got 



A/i//i = (-8.3;^-5) X 10" 
at l.Gcr level and 

A/i//i= (-1± 1.2) X 10" 



(231) 



(232) 



at 2(7 level. These results were confirmed by Potekhin et 
al. (1998) using 83 absorption lines to get 



A;u/;u = (-7.5 ± 9.5) X 10" 



(233) 



at a 2a level. 

More recently, Ivanchik et al. (2001) measured, with 
the VLT, the vibro-rotational lines of molecular hydrogen 
for two quasars with damped Lyman-a systems respec- 
tively at z = 2.3377 and 3.0249 and also argued for the 
detection of a time variation of fx. Their most conserva- 
tive result is (the observational data were compared to 
two experimental data sets) 



A;u/At= (-5.7±3.8) X 10" 



(234) 



at 1.5cr and the authors cautiously point out that addi- 
tional measurements arc necessary to ascertain this con- 
clusion. 1.5cr is not really significant and this may not 
survive further extended analysis. The result is also de- 
pendent on the laboratory dataset used for the compar- 
ison since it gave A/i//i = (—12.2 ± 7.3) x 10~^ with 
another dataset. 

As in the case of Webb et al. (1999, 2001), this mea- 
surement is very important in the sense that it is a 
non-zero detection that will have to be compared with 
other bounds. The measurements by Ivanchik et al. 
(2001) is indeed much larger than one would expect 
from the electromagnetic contributions. As seen above, 
the change in any unified theory, the changes in the 
masses are expected to be larger than the change in . 
Typically, we expect A/i//i ~ AAq^d/Aqcd - ^v/v ~ 
(30 — 40)Aq!e^5/q!jj,^5, so that it seems that the detection 
by Webb et al. (2001) is too large by a factor of order 10 
to be compatible with it. 
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Wiklind and Combes (1997) observed the quasar 
PKS 1413+135 with redshift z = 0.247 and used different 
transitions from the same molecule to constrain the varia- 
tion of /i. They compared different lines of HCO"'", HCN, 
CO and showed that the redshift difference are likely to 
be dominated by the velocity diflFerence between the two 
species which limits the precision of the measurements to 
A/i//x ~ 10"^ at 3(7 level. In one source (B3 1504+377) 
they observed a discrepancy of i^n/n ^ 10~^ 

Pagcl (1977, 1983) used another method to constrain jj, 
based on the measurement of the mass shift in the spec- 
tral lines of heavy elements. In that case the mass of 
the nucleus can be considered as infinite contrary to the 
case of hydrogen. A variation of /z will thus influence 
the redshift determined from hydrogen [see Eq. (68)] . He 
compared the redshifts obtained from spectrum of hydro- 
gen atom and metal lines for quasars of redshift ranging 
from 2.1 to 2.7. Since 

A^^^„-2„etal = (1 + ^)7^' (235) 

1 - Mo 

he obtained that 

|A/u//u| < 4 X IQ-i (236) 

at 3(7 level. This result is unfortunately not conclusive 
because usually heavy elements and hydrogen belong to 
different interstellar clouds with different radial velocity. 

D. Proton gyromagnetic factor 

As seen in Section III.B, the hyperfine structure in- 
duces a splitting dependent on gp^^o^^y^■ The ratio be- 
tween the frequency V21 of the hyperfine 21 cm absorption 
transition an optical resonance transition of frequency 
z/^pj mainly depends on 

!^2i/!^opt "L^pM = X. (237) 

By comparing the redshift of the same object determined 
from optical data and the 21 cm transition, one deduces 
that 

A-z = -^opt - -221 = (1 + z)/^x/x (238) 

SavedofT (1956) used the spectrum of Cygnus A and 
deduced that 

IS.X/X = (3 ± 7) X IQ-^ (239) 

at z ^ 0.057. Wolfe et al. (1976) discovered a BL Lac ob- 
ject (AO 0235+164) having the same redshift determined 
either by the 21 cm absorption line or by the ultraviolet 
doublet of Mg"*". Using that 

z^H/z^MgOt5pMaL(l-3M+---) (240) 



they concluded that 

Ax/a; = (5 ± 10) x 10"^ (241) 

at redshift of z = 0.5. They also got a constraint on the 

variation of gp/i by comparing the separation of Mg II 
doublet to hydrogen to get |Agp/x/s(p/i| < 6x 10"^. Wolfe 
and Davis (1979) used the 21 cm absorption lines of neu- 
tral hydrogen in front of the quasar QSO 1331+170 at 
a redshift z ~ 2.081. They determined that the cloud 
was at redshift z ~ 1.755. The agreement between the 
21 cm and optical redshifts is limited by the error in the 
determination of the optical redshift. They concluded 
that 

|Aa;/a;| < 2 x 10"^ (242) 

at a redshift z ~ 1.755 and from another absorber at 
redshift z ^ 0.524 around the quasar AO 0235+164 gives 

|Aa;/.T| < 2.8 x 10"*. (243) 

Tubbs and Wolfe (1980) used a set of four quasars 
among which MC3 1331+17 for which ^21 = 1.77642 ± 
2 X 10~^ is known with very high precision and deduced 
that 

|Aa;/a;| < 2 x 10"^. (244) 

Cowie and Songaila (1995) used the observations of C° 

absorption and fine structure to get the better optical 
redshift z , = 1.77644 + 2 x 10"^ which enables them to 

opt 

improve the constraint to 

Aa;/.T = (7 ± 11) x 10"^. (245) 

Besides the uncertainty in the determination of the opti- 
cal redshift, since the 21 cm optical depth depends sen- 
sitively on spin temperature while resonance-line opti- 
cal depths do not, the two regions of absorption need 
not coincide. This induces an uncertainty A.z = ±(1 + 
2;)(Au^p,/c) into Eq. (237) [see e.g. Wolfe and Davis 
(1979) for a discussion]. 

Drinkwater et al. (1998) compared the hydrogen hy- 
perfine structure to molecular absorption for three sys- 
tems at redshift z = 0.24,0.67 and 0.68 and used CO 
absorption lines. This allows to constrain y = g^a'^y^ 
and they got 

|Ay/y| < 5 X 10"^ (246) 

Assuming that the change in ^p and ^re not cor- 
related they deduced that |Ag(p/,gp| < 5 x 10^^ and 
\Aa^^/a^y^\ < 2.5 x 10~^. Varshalovich and Potekhin 
(1996) used the CO and hyperfine hydrogen redshift to- 
ward PKS 1413+135 {z = 0.247) to get 

Ay/y = (-4 ± 6) x 10"^ (247) 
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and PKS 1157+0.14 {z = 1.944) 

Ay/y = (7 ± 10) x 10"^ (248) 

Murphy et al. (2001c) improved the precision of this 

measurement by fitting Voigt profiles to the H 21 cm 
profile instead of using published redshifts and got 

Ay/y = (-0.2 ± 0.44) x 10"^ (249) 

at ^; = 0.25 and 

Ay/y = (-0.16 ± 0.54) x 10"^ (250) 

at z = 0.68. With the same systems Carrilli et al. (2001) 
found 

\Ay/y\ < 1.7 x 10"^ (251) 

both aXz = 0.25 and z = 0.68. Murphy et al. (2001c) ar- 
gued that one can estimate the velocity to 1.2 km • s~^ in- 
stead of the 10 km • s~^ assumed by Carrilli et al. (2001) 
so that their results in fact lead to Ay/y = (1 ± 0.03) x 
10-5 at z = 0.25 and Ay/y = (1.29 ± 0.08) x 10"^ at 
z = 0.68. 

E. The particular case of the cosmological constant 

The cosmological constant has also been loosing its 
status of constant. In this section, we briefly review the 
observations backing up this fact and then describe the 
theoretical models in favor of a time dependent cosmolog- 
ical constant and some links with the variation of other 
fundamental constants. 

The combination of recent astrophysical and cosmolog- 
ical observations [among which the luminosity distance- 
redshift relation up to z ~ 1 from type la supernovae 
(Riess et al., 1998; Perlmutter et al., 1998), the cos- 
mic microwave background temperature anisotropics (de 
Bernardis et al, 2000) and gravitational lensing (Mellier, 
1999)] seems to indicate that the universe is accelerating 
and that about 70% of the energy density of the universe 
is made of a matter with a negative pressure (i.e. having 
an equation of state w = P/p < 0). 

There are many different candidates to account for this 
exotic type of matter. The most simple solution would 
be a cosmological constant (for which w = —1) but one 
will then have to face the well known cosmological con- 
stant problem (Weinberg, 1989), i.e. the fact that the 
value of this cosmological constant inferred from the cos- 
mological observations is extremely small — about 120 
order of magnitude — compared with the energy scales of 
high energy physics (Planck, GUT and even electroweak 
scales). Another way is to argue that there exists a (yet 
unknown) mechanism which makes the cosmological con- 
stant strictly vanish and to find another matter candidate 



(referred to as "dark energy" ) able to explain the cosmo- 
logical observations. 

Among all the proposals (see e.g. Binetruy 2000 and 
Carroll. 2000 for a review) quintessence seems to be a 
promising mechanism. In these models, a scalar field is 
rolling down a runaway potential decreasing to zero at 
infinity hence acting as a fluid with an effective equation 
of state in the range — 1 < w < 1 if the field is mini- 
mally coupled. Runaway potentials such as exponential 
potential and inverse power law potentials 

V{(l)) = M'^+^/cp", (252) 

with a > and M a mass scale, arise in models where 

supersymmetry is dynamically broken (Binetruy, 1999) 
and in which flat directions are lifted by non-perturbative 
effects. 

One of the underlying motivation to replace the cosmo- 
logical constant by a scalar fleld comes from superstring 
models in which any dimensionful parameter is expressed 
in terms of the string mass scale and the vacuum expec- 
tation value of a scalar field. However, the requirement 
of slow roll (mandatory to have a negative pressure) and 
the fact that the quintessence field dominates today im- 
ply, if the minimum of the potential is zero, that (i) 
it is very light, roughly of order 10"^^ eV (Carroll, 
1998) and that (ii) the vacuum expectation value of the 
quintessence field today is of order of the Planck mass. 
It follows that coupling of this quintessence field leads to 
observable long-range forces and time dependence of the 
constant of nature. 

Carroll (1998) considered the effect of the coupling of 
this very light quintessence field to ordinary matter via 
an interaction of the form f3i{4>/M)Ci and to the elec- 
tromagnetic field as 4>F^"'F^^. Chiba and Kohri (2001) 
also argued that an ultra-light quintessence field induces 
a time variation of the coupling constant if it is coupled 
to ordinary matter and studied a coupling of the form 
(j)F'"'Ff,^. Dvafi and Zaldarriaga (2002) showed that it 
will be either detectable as a quintessence field or by tests 
of the equivalence principle, as also concluded by Wet- 
terich (2002). 

It was proposed that the quintessence field is also 
the dilaton (Uzan, 1999; Banerjee and Pavon, 2001; 
Esposito-Farcse and Polarski, 2001; Riazuelo and Uzan, 
2000; Gasperini et al, 2002). The same scalar field drives 
the time variation of the cosmological constant and of 
the gravitational constant and it has the property to also 
have tracking solutions (Uzan, 1999). 

Another motivation for considering the link between 
a dynamical cosmological constant and the time varia- 
tion of fundamental constants comes from the origin of 
the inverse power law potential. As shown by Binetruy 
(1999), it can arise from supersymmetry breaking by non 
perturbative effects such as gaugino condensation. The 
same kind of potential was also considered by Vayonakis 



37 



(1988) while discussing the variation of the fundamental 
couplings in the framework of 10-dimensional supergrav- 
ity. 

The variation of fundamental constants has also other 
implications on the measurement of the cosmological con- 
stant. Riazuelo and Uzan (2002) considered the effect of 
the variation of the gravitational constant on supernovae 
data. Besides changing the luminosity distance-redshift 
relation, the variation of G changes the standard picture, 
according to which type la supernovae are standard can- 
dles, in two ways. First the thermonuclear energy release 
proportional to the synthetized nickel mass is changing 
(and hence the maximum of the light curve); second the 
time scale of the supernovae explosion and thus the width 
of the light curve is also changed. Riazuelo and Uzan 
(2002) derived the modified magnitude-redshift relation 
to include the effect of the variation of G, using a one- 
zone analytical model for the supernovae and was con- 
firmed by numerical simulations (Gaztanaga et al., 2002). 

Barrow and Magueijo (2001) considered the effect of a 
time dependent fine structure constant on the interpreta- 
tion of the supernovae data. Their study was restricted 
to a class of varying speed of light theories (see Sec- 
tion VI. C) which have cosmological solutions very similar 
to quintessence. But, only the effect on the Hubble di- 
agram was studied and the influence of the change of 
the fine structure constant on the thermonuclear burst 
of the supernovae, and hence on its light curve, was not 
considered at all. 

Up to now there is no observational evidence of a time 
variation of the cosmological constant. The measurement 
of the equation of state of the dark energy can be hoped 
to be possible very soon, the best candidate method being 
the use of large-scale structure growth and weak gravi- 
tational lensing (Benabed and Bernardeau, 2001). But, 
it seems that the variation of constants and the dark en- 
ergy arc somehow related (Dvali and Zaldarriaga, 2002: 
Chiba and Khori, 2001; Banks et al, 2002; Wetterich, 
2002; Fujii, 2002), at least they share the properties to 
be very light and to appear in many models with a run- 
away potential. 



emitted by two quasars. Their idea is based on the re- 
mark that a prism is sensitive to the energy E of the pho- 
ton and a diffraction grating to its wavelength A so that 
any difference in the comparison of the wavelengths of a 
particular spectral line could be attributed to a change 
in %. Their study led to a null result in terms of experi- 
mental errors. 

Noerdlinger (1973) [and later Blake (1977a)] tried to 
measure E\. His argument was that the intensity of the 
Rayleigh- Jeans tail of the Planck spectrum of the CMB 
photon determines /c^ T whereas the turnover point of the 
spectrum determines hu/k^T. It follows that one can 
determine the value of he at the time of recombination, 
leading to the constraint |Aln/ic| < 0.3. 

Further works were performed by Solheim et al. (1973) 
and Baum and Florentin-Nielsen (1976) who compared 
the light of nearby and distant galaxies in order to test 
the constancy of EX. Bekenstein (1979) demonstrated 
that these experiments were meaningless since the con- 
stancy of EX was interpreted as the constancy of he but 
that this latter fact was implicitely assumed in the two 
experiments since the wave vector and momentum of the 
photon were both parallely propagated. This is only pos- 
sible if their proportionality factor he is constant, hence 
ensuring the null result of the two experiments. 

VI. THEORETICAL MOTIVATIONS 

One general feature of extra-dimensional theories, such 
as Kaluza-Klein and string theories, is that the "true" 
constants of nature are defined in the full higher di- 
mensional theory so that the effective 4-dimcnsional con- 
stants depends, among other things, on the structure and 
sizes of the extra-dimensions. Any evolution of these sizes 
either in time or space, would lead to a spacctime depen- 
dence of the effective 4-dimensional constants. 

We present in Sections VI. A and VI. B some results 
concerning Kaluza-Klcin theories and string theories. We 
end in Section VI. C by describing some phenomenological 
approaches initiated by Bekenstein (1982). 



F. Attempts to constrain the variation of dimenslonful 
constants 

As emphasized in Section I, considering the variation 

of dimenslonful constants is doubtful and seems meaning- 
less but such attempts have nevertheless been performed. 
We briefly review and comment them. These investiga- 
tions were mainly motivated by the construction of cos- 
mological models alternative to the big bang scenario and 
in which the redshift needs to have another interpreta- 
tion. 

Bahcall and Salpeter (1965) proposed to look for a time 
variation of the Planck constant by comparing the light 



A. Kaluza-Klein theories 

The aim of the early model by Kaluza (1921) and Klein 
(1926) to consider a 5-dimensional spacctime with one 
spatial extra-dimension (assumed to be of radius -R^k) 
was to unify electromagnetism and gravity (for a review 
see e.g. Overduin and Wesson, 1997). Starting from the 
5-dimensional Einstein- Hilbert Lagrangian 

S5 = IJ d^x^/^MfiJs, (253) 
we decompose the 5-dimensional metric as 
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dsi = 5,..da;''da;" + e^^ {A^dx'' + dy) 



(254) 



This form still allows 4-dimcnsional reparametrizations 
of the form y' = y + A (a;'') provided that = — 
dfj,\ so that gauge transformations arise from the higher 
dimensional coordinate transformations group. Any field 
(j) can be decomposed as 



(255) 



The 5-dimensional Klein-Gordon equation for a massless 
field becomes 



^2 



(256) 



so that has a mass m„ = n/R^^. At energies small 
with respect to m^^ = R^^, only y-independent fields 
remain and the physics is 4-dimcnsional. The effective 
action for the massless fields is obtained from the relation 
R5=Ri- 2e-^Ae^ - e'^^F^/A with F^, = d^A, - d,A^ 
so that 



^e rij_i,r 



(257) 



The field equations do not determine the compactifica- 

tion radius and only the invariant radius p = R^^^^ cxp{a) 
distinguishes non-equivalent solutions (one can set R^^ 
to unity without loss of generality). 



Sotting Af, 



the covariant derivative is 



9p + ipyAf^ = + inA^ so that the charges are inte- 
gers. The 4-dimensional Yang-Mills coupling, identified 
as the coefficient — 1/45'ym ^'^'^ gravitational con- 

stant are given by 



Ml = 2'KpMl, 



(258) 



2-Kp being the volume of the extra-space. Note that as 
long as one considers vacuum as in Eq. (254), there is 
a conformal undeterminacy that has to be lifted when 
adding matter fields. This generalizes to the case of D 
extra-dimensions (sec e.g. Cremmer and Scherk, f977 
and Forgacs and Horvath, 1979 for the case of two extra- 
dimensions) to 



G en p 



-D 



K,{D)Gp' 



(259) 



where the constants Ki depends only on the dimension 
and topology of the compact space (Weinberg, 1983b) 
so that the only fundamental constant of the theory is 
Mij^D- A theory on Ma x Md where Md is a D- 
dimensional compact space generates a low-energy quan- 
tum field theory of the Yang-Mills type related to the 
isometrics oi Md [for instance Witten (1981) showed 
that for £) = 7, it can accommodate the Yang-Mills group 



SU{3) X SU{2) X [/(!)]. Indeed the two main problems of 
these theories is that one cannot construct chiral fermions 
in four dimensions by compactification on a smooth man- 
ifold with such a procedure and that gauge theories in five 
dimensions or more are not renormalisable. 

The expression for the structure constants at lower en- 
ergy arc obtained by the renormalisation group (Mar- 
ciano, 1987; Wu and Wang, 1986) 



"i ^(mc^) = ^(m^KC^) - -^Ci- 



In- 



—0{m — rrij) In 



m 



(260) 



where the sum is over all Icptons, quarks, gluons... and 
the Cij are constants that depend on the spin and group 
representation (Georgi et al., 1974). Note however that 
this relation is obtained by considering the rcnormaliza- 
tion group in four dimensions and does not take into 
account the contribution of the Kaluza-Klein modes in 
loops. 

Chodos and Detweiler (1980) illustrated the effect of 
the fifth dimension by considering a 5-dimensional vac- 
uum solution of the Kasner form 



ds' 



-dt' 



(261) 



with J2pi = J2p1 = 1 ^'^d assuming compact spatial 
sections < Xi < L. In order to ensure local isotropy 
and homogeneity, they choose the solution pi = p2 = 
Ps = 1/2 and P4 = —1/2 so that the universe has four 
macroscopic spatial dimensions at the time to and looks 
spatially 3-dimensional at a time t ^ to with a small 
compact dimension of radius {To/tY^'^L. Considering 
as a small metric perturbation, they deduced that 



(262) 



hence offering a realization of Dirac large number hy- 
pothesis. Freund (1982) studied (4 -|- D) Kaluza-Klein 
cosmologies starting both in a (4 -|- £>)-dimensional Ein- 
stein gravity or a (4 -|- Z))-dimensional Brans-Dicke grav- 
ity. 

Using the expressions (259-260), Marciano (1984) re- 
lated the time dependence of the different couplings and 
restricted his discussion to the cases where Ki and rhj 
vanish. In the case where ai{mj^^) = (as studied in 
Chodos and Detweiler, 1980) one can relate the time vari- 
ation of the gravitational and fine structure constant as 



2tt 



2j 



G 
G' 
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In the case where cxiiiTij^j^) ^ (as studied in Freimd, 
1982), it was shown that the time variation of is en- 
hanced at low energy so that constraints on the time 
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variation of mc/m-p provide a sensitive test. It is also 
claimed that in the case of an oscillating the am- 
plitude of the oscillations will be damped by radiation 
in our 3-dimcnsional spacctimc due to oscillating charges 
and that experimental bounds can be circumvented. 

Kolb et al. (1985) used the variation (259) to constrain 
the time variation of the radius of the extra-dimensions 
during primordial nucleosynthesis (see section III.C.2) 
but their ansatz concerning the variation of Gp was ill- 
motivated. They deduced | Ai?j^j^/i?j^j^ | < 1%. Barrow 
(1987) took the effects of the variation of oc (see 
Section V.B) and deduced from the helium abundances 
that for \AR^^/R^J < 0.7% and \ AR^^/R^^\ < 1.1% 
respectively for £) = 2 and D = 7 Kaluza-Klein theory 
and that \ AR^^/R^^ | < 3.4 x 10"^° from the Oklo data. 

It follows that the radius of the extra-dimensions has 
to be stabilized but no satisfactory and complete mech- 
anism has yet been found. Li and Gott (1998) consid- 
ered a 5-dimensional Kaluza-Klcin inflationary scenario 
which is static in the internal dimension and expand- 
ing in the other dimensions and solve the 5-dimensional 
semi-classical Einstein equations including the Casimir 
effect. In particular, it was deduced that the effective 4- 
dimensional cosmological constant is related to the fine 
structure constant by GA^^ = {15gt,/204:8TT'^)a'^^. 



B. Superstring theories 

There exist five anomaly free, supersymmctric pc^rtur- 
bative string theories respectively known as type I, type 
IIA, type IIB, SO(32) heterotic and EgxEs heterotic the- 
ories (see e.g. in Polchinski, 1997). One of the definitive 
predictions of these theories is the existence of a scalar 
field, the dilaton, that couples directly to matter (Tay- 
lor and Vencziano, 1988) and whose vacuum expectation 
value determines the string coupling constant (Witten, 
1984). There are two other excitations that are common 
to all perturbative string theories, a rank two symmetric 
tensor (the graviton) g^^, and a rank two antisymmetric 
tensor 5^,^ . The field content then differs from one theory 
to another. It follows that the 4-dimensional couplings 
are determined in terms of a string scale and various dy- 
namical fields (dilaton, volume of compact space, ...). 
When the dilaton is massless, we expect three effects: (i) 
a scalar admixture of a scalar component inducing devia- 
tions from general relativity in gravitaional eff'ects, (ii) a 
variation of the couplings and (iii) a violation of the cak 
equivalence principle. Our purpose is to show how the 4- 
dimensional couplings are related to the string mass scale, 
to the dilaton and the structure of the extra-dimensions 
mainly on the example of heterotic theories. 

The two heterotic theories originate from the fact that 
left- and right-moving modes of a closed string are inde- 
pendent. This reduces the number of supersymmetry to 



A'^ = 1 and the quantization of the left-moving modes im- 
poses that the gauge group is either 5*0(32) or Es x Es 
depending on the fermionic boundary conditions. The 
effective tree-level action is 



Sh=J di°x^/=^e-2* [Ml {Rw + 4n$ - 4(V$)'} 

(264) 



FabF^'' + ... 



When compactified on a 6-dimensional Calabi-Yau space, 
the effective 4-dimensional action takes the form 




R4 + 




(265) 



where (j) = Vf,c couples identically to the Einstein and 
Yang-Mills terms. It follows that 



Mi = M 



Ml 



(266) 



at tree-level. Note that to reach this conclusion, one has 
to assume that the matter fields (in the 'dots' of Eq. (265) 
are minimally coupled to 34, see e.g. the discussion by 
Maeda, 1988). 

The strongly coupled SO (32) heterotic string theory 
is equivalent to the weakly coupled type I string the- 
ory. Type I superstring admits open strings, the bound- 
ary conditions of which divide the number of supersym- 
metries by two. It follows that the tree- level effective 
bosonic action isA^=l,_D=10 supergravity which 
takes the form, in the string frame. 



Si = j diOx,/^Mfe-* [e-*M2i?io 



■ + ... 



(267) 



where the dots contains terms describing the dynamics of 
the dilaton, fermions and other form fields. At variance 
with (264), the field $ couples differently to the gravi- 
tational and Yang-Mills terms because the graviton and 
Yang-Mills fields are respectively excitation of close and 
open strings. It follows that M/ can be lowered even 
to the weak scale by simply having exp $ small enough. 
Type I theories require D9-brancs for consistancy. When 
Ve is small, one can use T-duality (to render Ve large, 
which allows to use a quantum field theory approach) 
and turn the -D9-brane into a £>3-brane so that 



Si = j di°x,/=^e-2*Mf i?io 



(268) 
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where the second term describes the Yang-Mills fields 
localized on the Z)3-brane. It follows that 



-2$ 



3-2 



(269) 



at tree-level. If one compactifies the Z)9-brane on a 6- 
dimensional orbifold instead of a 6-torus, and if the brane 
is localized at an orbifold fixed point, then gauge fields 
couple to fields Mi living only at these orbifold fixed 
points with a (calculable) tree-level coupling Ci so that 



Ml = e-'^VeMf, 



+ CiMi. 



(270) 



The coupling to the field q is a priori non universal. At 
strong coupling, the 10-dimensional x Es heterotic 
theory becomes M-theory on x /Z2 (Hofava and 
Witten, 1996). The gravitational field propagates in the 
11-dimensional space while the gauge fields are localized 
on two 10-dimensional branes. 

At one-loop, one can derive the couplings by including 
Kaluza-Klein excitations to get (see e.g. Dudas, 2000) 



Ml 



ba 



{RM,f + ... 
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when the volume is large compared to the mass scale and 
in that case the coupling is no more universal. Otherwise, 
one would get a more complicated function. Obviously, 
the 4-dimensional effective gravitational and Yang-Mills 
couplings depend on the considered superstring theory, 
on the compactification scheme but in any case they de- 
pend on the dilaton. 

Wu and Wang (1986) studied the cosmological behav- 
ior of the theory (264) assuming a 10-dimensional metric 
of the form dia,g{-l, R3{t)'^gij{x), Re{t)'^gmn{y)) where 
Rs and Re are the scale factors of the external and in- 
ternal spaces. The rate of evolution of the size of the 
internal space was related to the time variation of the 
gravitational constant. The effect of a potential for the 
size of the internal space was also studied. 

Maeda (1988) considered the {N = 1,D = 10)- 
supergravity model derived from the heterotic super- 
string theory in the low energy limit and assumed that 
the 10-dimensional spacetime is compactified on a 6- 
torus of radius i?(x^) so that the effective 4-dimensional 
theory described by (265) is of the Brans-Dicke type 
with ui = —1. Assuming that <p has a mass /z, 
and couples to the matter fluid in the universe as 
^»atte, = /d'°xy=^exp(-2$)£„^,,^^(5io), the re- 
duced 4-dimensional matter action is 



-f 



.id). 



(272) 



The cosmological evolution of ^ and R can then be 
computed and Maeda (1988) deduced that aEM/o^EM — 
10^''(/x/l eV)-2 yr-i. In this approach, there is an ambi- 
guity in the way to introduce the matter fluid. 



Vayonakis (1988) considered the same model but 
assumed that supersymmetry is broken by non- 
perturbative effects such as gaugino condensation. In 
this model, and contrary to the work by Maeda (1988), 
(j) is stabilized and the variation of the constants arises 
mainly from the variation of i? in a runaway potential. 

Damour and Polyakov (1994a, 1994b) argued that the 
effective action for the massless modes taking into ac- 
count the full string loop expansion is of the form 




{Bg($)i? + 4B$($) 



(V$)2]} 
(273) 



in the string frame, Mg being the string mass scale. The 
functions Bi are not known but can be expanded as 



Bm 



-2<S> 
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in the limit $ —00, so that these functions can ex- 
hibit a local maximum. After a conformal transforma- 
tion {g^^i, = CBgg^^,ip = {CBg)-^/'^B]^'^ '(()), the action 
in Einstein frame takes the form 



d'^x 



167rG ^ 



R-2{yci>f-\Bp{cl^)F 



(275) 



from which it follows that the Yang-Mills coupling be- 
haves as (?-2 = kBp{(j)). This also implies that the QCD 
mass scale is given by 



A, 



M,{CBg) 



-l/2g-87r^/sSF/fe 



(276) 



where b depends on the matter content. It follows that 
the mass of any hadron, proportional to A^^^ in first ap- 
proximation, depends on the dilaton, mj^(Bg, Bp, . . .). 
With the anstaz (274), mA{4>) can exhibit a minimum 
(l)m that is an attractor of the cosmological evolution that 
drives the dilaton towards a regime where it decouples 
from matter. But, one needs to assume for this mecha- 
nism to apply, and particularly to avoid violation of the 
equivalence principle at an unacceptable level, that all 
the minima are the same, which can be implemented by 
setting Bi = B. Expanding In B around its maximum 
as Ini? oc —K{(t> — (t>„i)'^ /2, Damour and Polyakov (1994a, 
1994b) constrained the set of parameters {k, (/>o — (pm) us- 
ing the different observational bounds. This toy model 
allows to address the unsolved problem of the dilaton 
stabilization and to study all the experimental bounds 
together. 

Damour, Piazza and Veneziano (2002a, b) extended 
this model to a case where the coupling functions have 
a smooth finite limit for infinite value of the bare string 
coupling, so that Bi = Ci + 0{e~'^). The dilaton runs 
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away toward its attractor at infinity during a stage of in- 
flation. Tfie amplitude of residual dilaton interaction is 
related to the amplitude of the primordial density fluctu- 
ations and it can induce a variation of the fundamental 
constants, provided it couples to dark matter or dark en- 
ergy. It is concluded that, in this framework, the largest 
allowed variation of a.-^^^^ is of order 2 x 10~^, which is 
reached for a violation of the universality of free fall of 
of order 10~^^. 

Kolb et al. (1985) argued that in 10-dimensional su- 
perstring models, G oc and oc to deduce 
that |Ai?/i?| < 0.5%. This was revised by Barrow (1987) 
who included the effect of to deduce that helium 
abundances impose |Ai?/i?| < 0.2%. Recently Ichikawa 
and Kawasaki (2002) considered a model in which all 
the couplings vary due to the dilaton dynamics and con- 
strain the variation of the dilaton field from nucleosyn- 
thesis as -1.5 X 10-* < yi67fGA<?!. < 6.0 x lO""^. From 
the Oklo data. Barrow (1987) concluded that \AR/R\ < 
1.5 X 10-1°. 

To conclude, superstring theories offer a theoretical 

framework to discuss the value of the fundamental con- 
stants since they become expectation values of some 
fields. This is a first step towards their understanding 
but yet, no complete and satisfactory mechanism for the 
stabilization of the extra-dimension and dilaton is known. 



C. Other investigations 

Independently of string theory, Bekenstein (1982) for- 
mulated a framework to incorporate a varying fine struc- 
ture constant. Working in units in which Ti and c are 
constant, he adopted a classical description of the elec- 
tromagnetic field and made a set of assumptions to obtain 
a reasonable modification of Maxwell equations to take 
into account the effect of the variation of the elementary 
charge [for instance to take into account the problem of 
charge conservation which usually derived from Maxwell 
equations]. His eight postulates are that (1) for a con- 
stant a^^ electromagnetism is described by Maxwell the- 
ory and the coupling of the potential vector A^, to mat- 
ter is minimal, (2) the variation of a^^ results from dy- 
namics, (3) the dynamics of electromagnetism and a^^ 
can be obtained from an invariant action that is (4) lo- 
cally gauge invariant, (5) electromagnetism is causal and 
(6) its action is time reversal invariant, (7) the short- 
est length scale is the Planck length and (8) gravitation 
is described by a metric theory which satisfies Einstein 
equations. 

Assuming that the charges of all particles vary in the 
same way, one can set e = eoe(a;^) where e{x^) is a di- 
mensionless universal field (it should be invariant under 
e constant x e through a redefinition of cq). The elec- 
tromagnetic tensor generalizes to 



and the electromagnetic action is given by 



-1 
16^ 



(277) 



(278) 



The dynamics of e can be shown to derive from the action 

where t is length scale which needs to be small enough to 
be compatible with the observed scale invariance of elec- 
tromagnetism (i'pi < £ < 10" 1^ — 10" 1® cm around which 
electromagnetism merges with the weak interaction) . Fi- 
nally, the matter action for point particles of mass m 
takes the form Sm = Y. J[-mc^ + {e/c)uf'A^]j-^S^{x' - 
a;*(T))d^x where 7 is the Lorentz factor and r the proper 
time. 

Varying the total action gives the electromagnetic 

equation 



and the equation for the dynamics of e 



□ = — 

he 



'e—-—F F^'^ 



(280) 



(281) 



with a = X;wc27-i(53(a;* - x*(r))/y^. The Maxwell 
equation (280) is the same as electromagnetism in a 
material medium with dielectric constant and per- 
meability [this was the original description proposed 
by Fierz (1955) and Lichnerowicz (1955); see also Dicke 
(1964)]. 

On cosmological scales, it can be shown that the dy- 
namical equation for e can be cast under the form 



(a^e/e)' = -a^C^PmC^ 



(282) 



where Q = 0(10"^) is a dimensionless (and approxi- 

matively constant) measuring the fraction of mass in 
Coulomb energy for an average nucleon compared with 
the free proton mass and pm is the matter density. Since 
Pm OC a"^, Eq. (282) can be integrated to relate (e/e)o 
to £/£p, and the cosmological parameters. In order to 
integrate this equation, Bekenstein assumed that C was 
constant, which was a reasonable assumption at low red- 
shift. Livio and Stiavelli (1998) extended this analysis 
and got C = 1.2 X IQ-'^{X + 4/3F) where X and Y are 
the mass fraction of hydrogen and helium. 

Replacing the quantity in the brackets of the r.h.s. of 
Eq. (281) by CPmC^ with C = ^(lO-^), the static form 
Eq. (281) is analogous to the standard Poisson equation 
so that In e is proportional to the gravitational potential 



Ine 



c ^ 

47rc2 L 



(283) 



42 



from which it foUows that a test body of mass m and of 
electromagnetic energy E^^^ experiences an acceleration 
of a = - V$ - M-'^{dE^^/de)V€. 

From the confrontation of the results of the spatial and 
cosmological variation of e Bekenstein (1982) conduced, 
given his assumptions on the couplings, that "is a 
parameter, not a, dynamical variahle^^ . This problem was 
recently by passed by Olive and Pospelov (2001) who 
generalized the model to allow additional coupling of a 
scalar field = Bp((j)) to non-baryonic dark matter (as 
first proposed by Damour et al., 1990) and cosmological 
constant, arguing that in certain classes of dark matter 
models, and particularly in supersymmetric ones, it is 
natural to expect that (p would couple more strongly to 
dark matter than to baryon. For instance, supersym- 
metrizing Bekenstein model, (f) will get a coupling to the 
kinetic term of the gaugino of the form M^^(j)xdx so that, 
assuming that the gaugino is a large fraction of the sta- 
ble lightest supersymmetric particle, then the coupling 
to dark matter would be of order lO'^ — 10'' times larger. 
Such a factor could almost reconcile the constraint aris- 
ing from the test of the universality of free fall with the 
order of magnitude of the cosmological variation. This 
generalization of Bekenstein model relies on an action of 
the form 



+ 
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^d'^x (284) 



M|Ba(<^)A+-M^B^(</.)x^X 



where the sum is over proton []p= 7^(9^ — ieoA^)] and 
neutron [Ip = 7^9^]. The functions B can be expanded 
(since one focuses on small variations of the fine structure 
constant and thus of 0) as Bx = 1 + Cx<t> + (.xfj)'^ /2. It 
follows that q;em(0) ~ e^/A-KBp{(i)) so that /^a^^Ja^^^ = 
Cf4'+{^f — 2C|')'/>^/2. This framework extends the anal- 
ysis by Bekenstein (1982) to a 4-dimensional parameter 
space (M*, Cf, Cm, Ca)- It contains the Bekenstein model 
(Cf = —2, Ca = 0, Cm ~ IO^^'Cf), a Jordan-Brans-Dicke 
model (Cf = 0, Ca = -2v^2/2j+3, ^„ = -l/V4w + 6), 
a string-like model (Cf = — V^, Ca = V^, Cm — "\/2/2) 
so that A/a^^/a^^ = 3) and supersymmetrized Beken- 



stein model (Cf 



-2, C. 



2, Cn 



Cx 



so that 



Aq!jj,j^j/Qjj,^j ^ 5/^). lu aU the models, the imiversal- 
ity of free fall sots a strong constraint on (p/y/uj (with 
CO = M*/2Af|) and the authors showed that a small set 
of models was compatible with the cosmological variation 
and the equivalence principle tests. 

The constraint arising from the universality of free fall 
can be fulfilled if one sets by hand Bp — 1 oc [0 — (l>{0)]'^ 



where ct){0) is the value of the field today. It then follows 
that the cosmological evolution will drive the system to- 
ward a state in which (f) is almost stabilized today but 
allowing for cosmological variation of the constants of 
nature. In their two-parameter extension, Livio and Sti- 
avelli (1998) found that only variations of Aa^^/a^^^ of 
8 X 10~^ and 9 x 10"'' respectively for z < 5 and z < 1.6 
were compatible with Solar system experiments. 

The formalism developed by Bekenstein (1982) was 
also applied to the the strong interaction (Chamoun et 
al, 2000, 2001) by simply adding a term fabcA''^Al to 
describe the gluon tensor field G^^, fabc being the struc- 
ture constants of the non-Abelian group. It was also im- 
plemented in the braneworld context (e.g. Youm, 2001) 
and Magueijo et al. (2001) studied the effect of a varying 
fine structure constant on a complex scalar field undergo- 
ing an electromagnetic U{1) symmetry breaking in this 
framework. Armendariz-Picon (2002) derived the most 
general low energy action including a real scalar field that 
is local, invariant under space inversion and time rever- 
sal, diffeomormism invariant and with a U(l) gauge in- 
variance. This form includes the previous form (284) of 
Bekenstein's theory as well as scalar-tensor theories and 
long wavelength limit of bimctric theories. 

Recently Sandvik et al. (2001) claimed to have gener- 
alized Bekenstein model by simply redefining = eA^, 
fnv = d^i^tti,] and V' = In e so that the covariant derivative 
becomes = + It follows that the total action 

including the Einstein-Hilbert action for gravity the ac- 
tions (278) and (279) for the modified electromagnetism 
and normal matter takes the form 



5* = 



J ^/^d^x (£g,av + £mat + + -C^^e"''^) (285) 



with £^ = —{u)/2)dij,4'd'^4' so that the Einstein equa- 
tion are the "standard" Einstein equations with an ad- 
ditive stress-energy tensor for the scalar field ip. Indeed, 
Bekenstein (1982) did not take into account the effect of 
e (or 0) in the Friedmann equation and studied only the 
time variation of c in a matter dominated universe. In 
that sense Sandvik et al. (2002) extended the analysis 
by Bekenstein (1982) by solving the coupled system of 
Friedmann and Klein-Gordon equations. They studied 
numerically in function of C/'^SBm (with Wsbm = ^pj^) 
and showed that cosmological and astrophysical data can 
be explained with o^sbm = 1 if C ranges between 0.02% 
and 0.1% (that is about one order of magnitude smaller 
than Bekenstein's value based on the argument that dark 
matter has to be taken into account). An extension of the 
discussion of the cosmological scenarios was performed 
in Barrow et al. (2002c) and it was shown that a^^ 
is constant during the radiation era, then evolves loga- 
rithmically with the cosmic time during the matter era 
and then tends toward a constant during a curvature or 
cosmological constant era. The scalar-tensor case with 
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both varying G and a^,^ was considered by Barrow et al. 
(2002a,b). 

Sandvik et al. (2002), following Barrow and O 'Toole 

(2001) , estimated the spatial variations to be of order 
Alna^„ ~ 4.8 x lO-^GM/cV (Magueijo, 2001) to con- 
clude that on cosmological scale Alna^^ ~ 10~^ if 
GM/c^r ~ lO""*^, as expected on cosmological scales 
if {5T/T)^„^ ~ GM/c^r. On the Earth orbit scale, 
this leads to the rough estimate jVlnaj,^! 10~^^ — 
10^^^ cm^^ which is about ten orders of magnitude 
higher than the constraint arising from the test of the 
universality of free fall. Nevertheless, Magueijo et al. 

(2002) re-analyzed the violation of the universality of free 
fall and claimed that the theory is still compatible with 
equivalence principle tests provided that Cm ^ 1 for dark 
matter. This arises probably from the fact that only 

is varying while other constants are fixed so that the dom- 
inant factor in Eq. (146) is absent. 

Wetterich (2002) considered the effect of the scalar 
field responsible for the acceleration of the universe (the 

"cosmon") on the couplings arising from the coupling 
of the cosmon to the kinetic term of the gauge field as 
Zp{^)F'^ /4:. Focusing on grand unified theory, so that 
one gets a coupling of the form C = Zp{(j))TT{F'^)/A + 
iZ^{(l))%l)Ipil) and assuming a runaway expotenial poten- 
tial, he related the variation of a^,^, ttg, the nucleon 
masses to the arbitrary function Zp and to the 4>- 
dependent electroweak scale so that the different bounds 
can be discussed in the same framework. 

Chacko et al. (2002) proposed that the variation 
of the fine structure constant could be explained by a 
late second order phase transition at z ~ 1 — 3 (that 
is around T ^ 10"*^ eV) inducing a change in the 
vacuum expectation value of a scalar field. This can 
be implemented for instance in supersymmetric theo- 
ries with low energy symmetry breaking scale. This will 
induce a variation of the masses of electrically charge 
particle and. From the renormalization group equa- 
tion Q'jT^ij = a.^l^{K) + 5];,(&j+i/27r) ln(mi+i/mj) and as- 
suming that ftj^Jj (A) was fixed, one would require that 
Y^Smi/mi = 0(10^^) to explain the observations by 
Webb et al. (2001), so that the masses have to increase. 
Note that it will induce a time variation of the Fermi con- 
stant. Such models can occur in a large class of super- 
symmetric theories. Unfortunately, it is yet incomplete 
and its viability depends on the exitence of an adjustment 
mechanism for the cosmological constant. But, it offers 
new way of thinking the variation of the constants at odd 
with the previous analysis involving a rolling scalar field. 

Motivated by resolving the standard cosmological puz- 
zles (horizon, flatness, cosmological constant, entropy, 
homogeneity problems) without inflation, Albrecht and 
Magueijo (1998) introduced a cosmological model in 



which the speed of light is varying. Earlier related at- 
tempts were investigated by Moffat (1993a, 1993b). Al- 
brecht and Magueijo (1998) postulated that the Pried- 
mann equations are kept unchanged from which it fol- 
lows that the matter conservation has to be changed and 
get a term proportional to c/c. The flatness and hori- 
zon problems are not solved by a period of accelerated 
expansion so that, contrary to inflation, it does not offer 
any explanation for the initial perturbations (see however 
Harko and Mak, 1999). Albrecht and Magueijo (1998) 
considered an abrupt change in the velocity of light as 
may happen during a phase transition. It was extended 
to scenarios in which both c and G were proportional to 
some power of the scale factor by Barrow (1999) (see also 
Barrow and Magueigo, 1999a, b). The link between this 
theory and Bekenstein theory was investigated by Bar- 
row and Magueijo (1998). Magueijo et al, (2002) inves- 
tigated the test of universality of free fall. A Lagrangian 
formulation would probably requires the introduction of 
an "ether" vector field to break local Lorentz invariance 
as was used in e.g. Lubo et al. (2002). 

Clayton and Moffat (1999) implemented a varying 
speed of light model by considering a bimetric the- 
ory of gravitation in which one metric g^^ describe the 
standard gravitational vacuum whereas a second metric 
Qnv + P^fi^iy^ P being a dimensionless constant and ip^^ a 
dynamical vector field, describes the geometry in which 
matter is propagating (see also Bekenstein, 1993). When 
choosing ^'/^ = 9^'^ this reduces to the models developed 
by Clayton and Moffat (2000, 2001). Some cosmological 
implications were discussed by Moffat (2001, 2002) but 
no study of the constraints arising from Solar system ex- 
periments have been taken into account. Note that Dirac 
(1979) also proposed that a varying G can be reconciled 
with Einstein theory of gravity if the space metric was 
different from the "atomic" metric. Landau and Vucetich 
(2000) investigated the constraints arising from the viola- 
tion of the charge conservation. Other realizations arise 
from the brane world picture in which our universe is 
a 3-dimensional brane embedded in a higher dimensional 
spacetime. Kiritsis (1999) showed that when a test brane 
is moving in a black hole bulk spacetime (Kehagias and 
Kiritsis, 1999) the velocity of light is varying as the dis- 
tance between the brane and the black hole. Alexander 
(2000) generalizes this model (see also Steer and Parry, 
2002) by including rotation and expansion of the bulk so 
that the speed of light gets stabilized at late time. Carter 
et al. (2001) nevertheless showed that even if a Newton- 
like force is recovered on small scales such models are very 
constrained at the post-Newtonian level. Brane models 
allowing for the scalar field in the bulk naturally predicts 
a time variable gravitational constant (see e.g. Brax and 
Davis, 2001). 



44 



D. A new cosmological constant problem? 

The question of the compatibihty between an ob- 
served variation of the fine structure constant and par- 
ticle physics models was put forward by Banks et al. 
(2002). As seen above, in the low energy limit, the change 
of the fine structure constant can be implemented by cou- 
pling a scalar field to the photon kinetic term F^'^F^j^, 
but this implies that the vacuum energy computed in this 
low energy limit must depend on . Estimating that 

AA_~A*Aa,„, (286) 

leads to a variation of order AA^.^^. ^ 10^^ ("^V)^ for a vari- 
ation Attj,,, - 0(10"") and foTA = A^^^, - 100 McV. 
Indeed, this contrasts with the average energy density 
of the universe of about 10** (eV)" during the matter era 
so that the universe was dominated by the cosmological 
constant at z ~ 3, which is at odds with observations. It 
was thus concluded that this imposes that 

|Aa™/«KMl <10-''. (287) 

Contrary to the standard cosmological constant problem, 
the vacuum zero-point energy to be removed is time- 
dependent and one can only remove it for a fixed value 
of a^^ . Whereas the cosmological constant problem in- 
volves the fine tuning of a parameter, this now implies 
the fine tuning of a function! 

It follows that a varying a^^ cannot be naturally ex- 
plained in a field theory approach. A possible way out 
would be to consider that the field is in fact an axion 
(see Carroll, 1998; Choi, 2000; Banks and Dine, 2001). 
Some possible links with Heisenbcag relations and quan- 
tum mechanics were also investigated by Raiiada (2002). 
Besides, the resolution of the cosmological constant prob- 
lem may also provide the missing elements to understand 
the variation of the constants. Both preoblems can be 
hoped to be solved by string theory. 

VII. CONCLUSIONS 

The experimental and observational constraints on the 
variation of the fine structure, gravitational constants, of 
the electron to proton mass ration and different combi- 
nations of the proton gyromagnetic factor and the two 
previous constants, as well as bounds on are summa- 
rized in tables II, IV and V. 

The developments of high energy physics theories such 
as multi-dimensional and string theories provide new mo- 
tivations to consider the time variation of the fundamen- 
tal constants. The observation of the variability of these 
constants constitutes one of the very few hope to test di- 
rectly the existence of extra-dimensions and to test these 
high energy-physics models. In the long run, it may help 
to discriminate between diflferent effective potentials for 



the dilaton and/or the dynamics of the internal space. 
But indeed, independently of these motivations, the un- 
derstanding of the value of the fundamental constants of 
nature and the discussion of their status of constant re- 
mains a central question of physics in general: question- 
ing the free parameters of a theory accounts to question 
the theory itself. It is a basic and direct test of the law 
of gravity. 

As we have shown, proving that a fundamental con- 
stant has changed is not an obvious task mainly because 
observations usually entangle a set of constants and be- 
cause the bounds presented in the literature often as- 
sume the constancy of a set of parameters. But, in GUT, 
Kaluza-Klein and string inspired models, one expects all 
the couplings to vary simultaneously. Better analysis 
of the degeneracies as started by Sisterna and Vucetich 
(1990, 1991) (see also Landau and Vucetich, 2002) are 
really needed before drawing definitive conclusions but 
such analysis are also dependent in the progresses in our 
understanding of the fundamental interactions and par- 
ticularly of the QCD theory and on the generation of the 
fermion masses. 

Other progresses require (model-dependent) investiga- 
tions of the compatibility of the different bounds. It 
has also to be remembered that arguing about the non- 
existence of something to set constraints (e.g. Broulik 
and Trefil, 1971) is very dangerous. On an observational 
point of view, one needs to further study the systemat- 
ics (and remember some erroneous claims such as those 
by Van Flandern, 1975) and to propose new experiments 
(sec e.g. Karshenboim, 2000, 2001 who proposed experi- 
ments based on the hyperfine structure of deuterium an 
ytterbium-171 as well as atoms with magnetic moment; 
Braxmaier et al., 2001; Torgerson, 2000 who proposed to 
compare optical frequency references; Sortais et al, 2001 
who improved the sensitivity of frequency standards, the 
coming satellite experiments ACES, MICROSCOPE and 
STEP...). On local scales, the test of the imiversal- 
ity of free fall sets drastic constraints and one can hope 
to use similar methods on cosmological scales from the 
measurements of weak gravitational lensing (Uzan and 
Bernardeau, 2001) or from structure formation (Mar- 
tins et al, 2002). The complementarity between local 
experiments and geo-astrophysical observations is neces- 
sary since these methods test different time-scales and 
are mainly sensitive either to rapid oscillations or a slow 
drift of the constants. 

The recent astrophysical observations of quasars tend 
to show that both the fine structure constant and the 
electron to proton mass ratio have evolved. These two 
measurements are non-zero detections and thus very 
different in consequences compared with other bounds. 
They draw the questions of their compatibility with the 
bounds obtained from other physical systems such as e.g. 
the test of the universality of free fall and Oklo but also 
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on a more theoretical aspect of the understanding of such 
a late time variation which does not seem to be natural 
from a field theory point of view. Theoretically, one ex- 
pects all constants to vary and the level of their variation 
is also worth investigating. One would need to study the 
implication of these measurements for the other experi- 
ments and try to determine their expected level of detec- 
tion. Both results arise from the observation of quasar 
absorption spectra; it is of importance to ensure that all 
systematics are taken into account and are confirmed by 
independent teams, using e.g. the VLT which offers a 
better signal to noise and spectral resolution. 

The step from the standard model-|-general relativity 
to string theory allows for dynamical constants and thus 
starts to address the question of why the constants have 
the value they have. Unfortunately, no complete and 
satisfactory stabilization mechanism is known yet and 
we have to understand why, if confirmed, the constants 
are still varying and whether such a variation induces a 
new cosmological constant problem. 

The study of the variation of the constants offers a new 
link between astrophysics, cosmology and high-energy 
physics complementary to primordial cosmology. It is 
deeply related to the test of the law of gravitation, both 
of the deviations from general relativity and the violation 
of the weak equivalence principle. But yet much work is 
needed both to disentangle the observations and to relate 
them to theoretical models. 
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FIG. 1. The cube of physical theories as presented by Okun 
(1991). At the origin stands the part of Newtonian mechan- 
ics (NM) that does not take gravity into account. NG, QM 
and SR then stand for Newtonian gravity, quantum mechanics 
and special relativity which respectively introduce the effect of 
one of the constants. Special relativity 'merges' respectively 
with quantum mechanics and Newtonian gravity to give quan- 
tum field theory (QFT) and general relativity (GR). Bringing 
quantum mechanics and Newtonian gravity together leads to 
non-relativistic quantum gravity and all theories together give 
the theory of everything (TOE). [From Okun (1991)]. 
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ing to the determination of the electron mass. Note that, as 
expected, the determination of a^.^ requires no dimensional 
input. [From Mohr and Taylor (2001)]. 
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FIG. 3. The standard orbital parameters, a and b is the semi- 
major and semi-minor axis, c = ae the focal distance, p the 
semi-latus rectum, 9 the true anomaly. F is the focus, A 
and B the periastron and apoastron (see e.g. Murray and 
Dermott, 2000). It is easy to check that — {I ~ e^) and 
that p — a(l — e^) and one defines the frequency or mean 
motion as n — 2-k/P where P is the period. 
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FIG. 4. Hyperfine structure of the n = 1 level of the hydro- 
gen atom. The fine structure Hamiltonian induces a shift of 
— of the level Is. J can only take the value -1-1/2. 
The hyperfine Hamiltonian (74) induces a splitting of the level 
lsi/2 into the two hyperfine levels F = and F = -1-1. The 
transition between these two levels corresponds to the 21 cm 
ray with Ah^ = 1, 420, 405, 751.768 ± 0.001 Hz and is of first 
importance in astronomy. 
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FIG. 6. Constraints on the coupling (upper panel) and ^/ 
(lower panel) respectively to N -\- Z and N — Z as a. function 
of the length scale A. The shaded regions are excluded at 2a. 
[From Fischbach and Talmadge (1997)]. 



Redshift 




10 20 30 40 50 60 70 80 90 
ATOMIC NUMBER Z 



FIG. 5. The correction function [From Prestage et al. 

(1995)]. 
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FIG. 7. AQfj,„/aj,„ as a function of the look-back time com- 
puted with the cosmological parameters (f2m, ^a) = (0.3, 0.7) 
and h = 0.68. Squares refer to the data by (Murphy et al. 
(2001c) assuming g-p constant; triangles refer to Si IV systems 
by Murphy et al. (2001d); dots correspond to Mg II and Fe II 
systems for redshidts smaller than 1.6 (Webb et al., 2001) 
and higher redshifts come from Murphy et al. (2001a). [From 
Webb et al. (2001)]. 
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TABLE I. Comparison of absorption lines and the combinations of the fundamental constants that can be constrained. 
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v^ow iV, oi IV,... j 


(Oowic (irid SongcLils., 1995) 


(-0.3 ± 1.9) 


X 10-* 


2.785 - 3.191 






\^jr 1 oaid^t; c-t III. , iyyo i 


(0± 1.42) 


X 10-" 





140 days 


iT-tOIiiiC CiUivb 


^^j_7ciiiiuiJ.i diiiu. j-zyouii, ±yyu J 


(0.15 ± 1.05) 


X 10"'^ 




1.8 




(Varshalovich al.^ 1996a) 


(2 ± 7) 


X 10"^ 


2.8-3.1 




\^D\J [pi iV J 


(TKaTrvQVrr\Yr\ oi nl 1 QQQ^ 

1 J3ci ^oLioiii a (iL.j lyyy j 


(0±2) 


X 10-2 


10* 






(\\cr)D et ai., iyyyj 


(-0.17 ±0.39) 


X 10--^ 


0.6-1 




i^oU (^iVlg ii, rc iij 


l^vveDD et (it., lyyy^ 


(-1.88 ±0.53) 


X 10-^ 


1 - 1.6 




i^ivig ii, re iij 


1^1 vdiiiciiiK. et tt(. , ±yyyy 


(-3.3 ±6.5 ±8) 


X 10-^ 


2 — 3.5 




090 TV^ 


(Fujii et al, 2000) 


(-0.36 ± 1.44) 


X 10-'^ 




1.8 


Oklo 


(Varshalovicli et at, 2000a) (- 


-4.5 ± 4.3 ± 1.4) 


X lO-'"' 


2-4 




QSO (Si IV) 


(Avelino et al, 2001) 


(-3.5 ±5.5) 


X 10-2 


10^^ 




CMB 


(Landau et al, 2001) 


(-5.5 ±8.5) 


X 10-2 


10^ 




CMB 


(Webb et al., 2001) 


(-0.7 ±0.23) 


X 10-^ 


0.5 - 1.8 




QSO (Fe II, Mg II) 


(Webb et al, 2001) 


(-0.76 ±0.28) 


X 10-^ 


1.8-3.5 




QSO (Ni II, Cr II, Zn II) 


(Webb et al, 2001) 


(-0.5 ±1.3) 


X 10-^ 


2-3 




QSO (Si IV) 


(Murphy et al., 2001a) 


(-0.2 ±0.3) 


X 10-^ 


0.5 - 1 




QSO (Mgl, Mg 11,...) 


(Murphy et al., 2001a) 


(-1.2 ±0.3) 


X 10-^ 


1 - 1.8 




QSO (Mg I, Mg II,. . . ) 


(Murphy e< a/., 2001a) 


(-0.7 ±0.23) 


X 10-® 


0.5 - 1.8 




QSO (Mg I, Mg II,. . . ) 


(Murphy et al, 2001d) 


(-0.5 ±1.3) 


X 10-® 


2-3 




QSO (Si IV) 


(Sortais et al., 2001) 


(8.4 ± 13.8) 


X 10-^® 




24 months 


Atomic clock 


(NoUet and Lopez, 2002) 


(3 ±7) 


X 10-2 


10* 




BBN 


(Ichikawa and Kawasaki, 2002) 


(-2.24 ±3.75) 


X 10-* 


10* 




BBN 


(Olive et ai, 2002) 


(0± 1) 


X 10-'^ 




1.8 


Oklo 


(Olive et ai, 2002) 


(0±3) 


X lO-'' 


~ 0.45 


4.6 


Re/Os 


(Olive et al, 2002) 


(0± 1) 


X 10-® 






a-decay 
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TABLE III. The different atomic clock experiments. We recall the transitions which axe compared and the constraint on 
the time variation obtained. SCO refers to superconductor cavity oscillator and the reference to (Breakiron, 1993) is cited in 
Prestage et al. (1995). fs and hfs refer respectively to fine structure and hyperfine structure. 



Reference 


Experiment 


Constant 


Duration 


Limit (yr ^) 


(Turneature and Stein, 1974) 


hfs of Cs vs SCO 




12 days 


< 1.5 X 10"'^ 


(Godone et al, 1983) 


hfs of Cs vs fs of Mg 




1 year 


< 2.5 X 10"" 


(Dcmidov, 1992) 


hfs of Cs vs hfs of H 


"emAp/A^ 


1 year 


< 5.5 X 10"" 


(Breakiron, 1993) 


hfs of Cs vs hfs of H 


ObmSp/S/ 




< 5 X 10"" 


(Prestage et al., 1995) 


hfs of HG+ vs hfs of H 


"EMfp/fl'/ 


140 days 


< 2.7 X 10"" 


(Sortais et al, 2001) 


hfs of Cs vs hfs of Rb 




24 months 


(4.2 ±6.9) X 10"^^ 



TABLE IV. Summary of the constraints on the time variation of the Newton constant G. The constraints labelled by * refer 
to bounds on the rate of decrease of G (that is —G/G <...). 





aiiit lyi I 


ivieLnuQ 


(^ieller, iy4oj 


{0 ± 2.0) X lU 


Earth temperature 


l^oiiapiro at cll., iy/ij 


(U ± 4j X lU 


Planetary ranging 


(^Morison, lijio) 


[U ± Z) X iU 


Lunar occult ations 


(Dearborn anci SchrcLinm, 1974) 


* ^ /I V 1 


Clusters of galaxies 


(van r lanaern, iy/oj 


/ Q _l_ c;^ 1 n~ll 
\—o ± Oj X iU 


Lunar occultations 


(Heintzmann and Hillebrandt, 1975) 


(U ^ IJ X lU 


Pulsar spin-down 


(Rcasenberg and Shapiro, 1976) 


(U =t l.Oj X lU 


Planetary ranging 


(^iViansneiQ, iy/oj 


< ( — O.ozhlj X lU 


Pulsar spin-down 


(^vviiuams et ctt. , lyyoj 


(U ± o j X lU 


Planetary ranging 


(oiaKe, ly/ lO) 


( — U.O ^ Zj X iU 


Earth radius 


l^iviuiier, lyio) 


yz.o It i.o^ X iu 


Solar eclipses 


{^Mcrjinmny et at., ly^oj 


< X iU 


Planetary radii 


(^jjarrow, ly/oj 


(z it y.o jAi X lu 




(Rcasenberg et al., 1979) 


< iU 


Viking ranging 


(van r lanaern, iysi J 


(o.Z it i.lj X lU 


Lunar occultation 


(Rothman and Matzner, 1981) 


(0± 1.7) X 10"" 


BBN 


(Helhngs et al., 1983) 


(2 ±4) X 10"^^ 


Viking ranging 


(Reasenberg, 1983) 


(0±3) X 10"" 


Viking ranging 


(Damour et al., 1988) 


(1.0 ±2.3) X 10"" 


PSR 1913+16 


(Shapiro, 1990) 


(-2± 10) X 10"^^ 


Planetary ranging 


(Goldman, 1990) 


* < (3.85 ±1.65) X 10"" 


PSR 0655+64 


(Accctta et al, 1990) 


(0 ± 9) X 10"" 


BBN 


(Miiller et al, 1991) 


(0± 1.04) X 10"" 


Lunar laser ranging 


(Anderson et al., 1991) 


(0.0 ±2.0) X 10"^^ 


Planetary ranging 


(Damour and Taylor, 1991) 


(1.10 ± 1.07) X 10"" 


PSR 1913+16 


(Chandler, 1993) 


(0± 1) X 10"" 


Viking ranging 


(Dickey et al, 1994) 


(0±6) X 10"^^ 


Lunar laser ranging 


(Kaspi et al., 1994) 


(4 ±5) X 10"i2 


PSR B1913+16 


(Kaspi et al., 1994) 


(-9 ±18) X 10"^^ 


PSR B1855+09 


(Demarque et al, 1994) 


(0 ± 2) X 10"" 


Heliosismology 


(Gucnthcr et al, 1995) 


(0±4.5) X 10"^^ 


Heliosismology 


(Garcia-Bcrro et al, 1995) 


* < {3tl) X 10~" 


White dwarf 


(WiUiams et al, 1996) 


(0 ± 8) X 10"^^ 


Lunar laser ranging 


(Thorsett, 1996) 


(-0.6 ±4.2) X 10"^^ 


Pulsar statistics 


(Dcl'Innoccnti et at, 1996) 


(-1.4 ±2.1) X 10"" 


Globular clusters 


(Guenther et al, 1998) 


(0± 1.6) X 10"^^ 


Hehosismology 
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TABLE V. Summary of the constraints on the variation of the constant k. We use the notation = me/rrip, x = a^^Qpfi and 



Reference 


Constant 


Constraint 


redshift 


Time (lO'' yr) 


Method 


(Yahil, 1975) 


M 


(0± 1.2) 




10 


Rb-Sr, K-Ar 


(Pagel, 1977) 




(0±4) X 10-1 


2.1-2.7 




QSO 


(Foltz et al, 1988) 




(0 ± 2) X IQ-* 


2.811 




QSO 


(Varshalovich and Lcvshakov 1993) 




(0±4) X 10"^ 


2.811 




QSO 


(Cowic and Songaila, 1995) 




(0.75 ±6.25) X 10""* 


2.811 




QSO 


(Varshalovich and Potekhin, 1995) 


M 


(0±2) X 10-* 


2.811 




QSO 


(Varshalovich et al., 1996a) 


M 


(0±2) X 10"" 


2.811 




QSO 


(Varshalovich et al., 1996b) 




(-1 ± 1.2) X 10"* 


2.811 




QSO 


(Potekhin et al, 1988) 




(-7.5 ±9.5) X 10"^ 


2.811 




QSO 


(Ivanchik et al, 2001) 




(-5.7 ±3.8) X 10"^ 


2.3-3 




QSO 


(Savedoff, 1956) 


X 


(3 ± 7) X 10"* 


0.057 




Cygnus A 


(Wolfe et al, 1976) 


X 


(5 ± 10) X 10"^ 


~ 0.5 




QSO (Mg I) 


(Wolfe and Davis, 1979) 


X 


(0 ± 2) X 10"* 


1.755 




QSO 


(Wolfe and Davis, 1979) 


X 


(0±2.8) X 10"* 


0.524 




QSO 


(Tubbs and Wolfe. 1980) 


X 


(0 ± 1) X 10"* 


1.776 




QSO 


(C 'owie and Scjugaila. 1995) 


,r 


(7± li) X 10 


1.776 




QSO 


(Varshalovich and Potekhin, 1996) 


y 


(-4 ±6) X 10"" 


0.247 




QSO 


(Varshalovich and Potekhin, 1996) 


y 


(-7± 10) X 10"^ 


1.94 




QSO 


(Drinkwater et al, 1998) 


y 


(0±5) X lO"*^ 


0.25, 0.68 




QSO 


(Carrilh et al, 2001) 


y 


(0 ± 3.4) X 10"^ 


0.25, 0.68 




QSO 


(Murphy et al, 2001c) 


y 


(-0.2 ±0.44) X 10"^ 


0.25 




QSO 


(Murphy et al, 2001c) 


y 


(-0.16 ±0.54) X 10"^ 


0.68 




QSO 


(Wolfe et al, 1976) 




(0 ± 0.68) X 10 


0.524 




QSO 


(Turneature and Stein, 1994) 




(0±9.3) X 10""' 




12 days 


Atomic clocks 


(Godone et al, 1993) 




(0±5.4) X 10"" 




1 year 


Atomic clocks 


(Wilkinson, 1958) 


"w 


(0± 1) X 10' 




1 


Fission 


(Dyson, 1972) 


"w 


(0± 1) X 10"' 




1 


/3-decay 


(Shlyakhter, 1976) 


"w 


(0±4) X 10"^ 




1.8 


Oklo 


(Damour and Dyson, 1996) 


"w 


(0 ± 2) X 10-2 




1.8 


Oklo 
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